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FOREWORD 


This  document  consists  of  the  papers  presented  at  the 
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28  January  1975. 
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CALL  TO  ORDER  AND  OPENING  REMARKS 


V.  J.  Corcoran 

I want  to  welcome  all  of  you  to  the  Workshop  on  Atmospheric  Transmission 
Modeling. 

The  purpose  of  this  workshop  Is  to  bring  together  those  people  who  are 
actively  Involved  In  computer  modeling  of  the  atmosphere  for  optical  propa- 
gation. This  Includes  people  who  have  generated  models,  those  who  have  pro- 
vided Inputs  to  the  models,  and  various  users  of  the  models.  Also  Included 
are  persons  who  have  been  critics  of  models  and  have  expressed  dissatisfaction 
with  models  as  they  now  stand.  It  Is  hoped  that  the  meeting  will  help  to  show 
what  the  similarities  and  differences  of  the  models  are,  and  that  the  criticisms 
that  have  been  expressed  can  be  aired  by  the  group  and  especially  by  those  who 
have  generated  the  models. 

As  a result  of  the  discussions  means  for  Improving  current  models  and 
desired  extensions  of  the  models  should  become  apparent.  This  future  work 
should  eventually  result  In  a model,  or  models,  that  can  be  used  to  do  system 
calculations  which  will  not  be  In  doubt  because  of  the  transmission  model 
used. 

We  hope  that  this  meeting  Is  different  from  other  meetings  concerned  with 
atmospheric  propagation.  If  It  Is  not  It  will  be  a waste  of  time.  The  best 
thing  that  could  possibly  happen  Is  for  us  to  conclude  that  no  further  meetings 
are  necessary  because  of  the  progress  that  has  been  made  toward  having 
satisfactory  models. 

As  you  can  see  from  the  agenda  we  shall  start  In  the  morning  with  three 
users  papers  followed  by  a variety  of  papers  on  modeling  and  finally  a comparison 
of  two  models.  The  afternoon  is  concerned  with  two  workshops,  one  on  the  physics 
and  engineering  aspects  end  the  other  on  the  computational  aspects  of  modeling 
the  atmosphere.  The  discussions  In  the  afternoon  are  expected  to  be  the  most 
Important  aspect  of  this  meeting.  After  the  two  workshops  we  shall  reconvene 
here  In  the  IDA  Board  Room  for  the  summary  reports  and  closing  remarks. 
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ABSTRACT 

This  paper  examines  the  detection  probability  of  tanks  versus 
range.  It  employs  real  weather  data,  a model  developed  by  Sendall 
and  Biberman  growing  out  of  the  Rotell  perceived  slgnal-to-nolse 
concept,  and  the  LOWTRAN  III  atmospheric  model. 

The  results  show  large  variability  In  probability  of  detection 
caused  by  the  large  and  frequent  variations  In  weather  examined 
hourly  for  a full  year. 


SUMMARY 

The  traditional  way  of  thinking  about  the  operational  utility  of  Infrared  sensors  such  os  FLIRs 
may  very  well  be  both  Inappropriate  and  Incorrect. 

The  calculations  based  upon  long-time  averages  of  weather  do  not  help  one  understand  the 
probable  utility  at  soma  particular  hour. 

Weather  does  change  rapidly,  and,  more  to  the  point,  for  Infrared  modeling  It  changes  drasti- 
cally within  a single  day  often  In  a cyclic  manner  from  day  to  day  to  day. 

Further,  aerosols  dominate  the  FLIR  performance  far  more  than  the  absorption  effects  when  one 
calculates  the  variations  In  FLIR  performance  that  might  be  expected. 
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1 , INTRODUCTION  AND  ACKNOWLEDGEMENTS 


The  Aerospace  Applications  Studies  Committee  of  NATO  AGARD  sponsored  a study  on  the  use  of 
sensors  for  low-flying  fast  aircraft  at  night.  In  the  course  of  this  study,  It  became  dear  that 
although  adequate  methodology  was  developed  to  carry  out  the  study,  the  weather  picture  assumed 
appeared  far  too  uniform,  and  Indeed  It  was  a highly  smoothed  ten-year  average.  Curiosity  about 
the  variations  of  unsmoothed  data  led  to  a proposal  by  L.  M.  Blberman  of  IDA  and  M.  Caller  of  RAE, 
Farnborough  to  the  AASC  that  this  problem  be  Investigated  In  tume  detail  to  study  the  effects  of 
terrain  masking  for  a number  of  real  locations,  the  effect  of  cloud  obscuration,  and  the  effect  of 
hour  by  hour  weather  variations. 

This  preliminary  report  outlines  the  early  findings  from  calculations  so  far  carried  out  for 
two  sample  months,  January  and  August.  These  results  are  the  dally  and  seasonal  effects  of  real 
-weathor  In  the  vicinity  of  Hannover,  Germany  that  tend  to  control  the  applicability  of  the  electro- 
optical  sensors,  The  question  Is  how  much  and  how  often.  It  is  our  hope  to  answer  those  questions. 

This  reporting  of  the  work  by  Blberman  and  du  Mala  Is  actually  based  upon  a larger  cooperative 
effort  with  a long  history.  The  final  reporting  will  be  by  a committee.  This  preliminary  report- 
ing must  however  acknowledge  the  foundations  upon  which  this  work  rest*. 

The  model  used  In  the  study  was  created  by  the  cooperation  of  a number  of  Individuals  supported 
by  their  parent  organizations,  The  basic  FUR  model  was  developed  by  Robert  L.  Sendall  of  the 
Xerox  Corporation's  Electro-Optical  Systems  Division  (now  of  Hughes  Aircraft  Company)  and 
Luden  M.  Blberman  of  IDA. 

The  weather  data  was  selected  and  assembled  by  ETAC  through  the  efforts  of  Major  Thomas  E, 
Stanton  and  Clarence  B.  Elam,  This  data  was  further  processed  oy  the  Westlnghouse  Aerospace  1 
Electronic  Systems  Division  of  Baltimore  under  the  supervision  of  Kenneth  C.  Leonard,  Jr. 

The  need  for  cloud-free  line  of  sight  data  was  long  unforeseen,  end  a cooperative  effort  by 
Colonel  John  T.  McCabe  of  ETAC  and  Blberman  of  IDA  led  to  the  work  "Estimating  Mean  Cloud  and 
Cl Imataloglca!  Probability  of  Cloud-Free  Lines  of  Sight"  In  1965  (Ref.  1). 

In  an  effort  to  extend  tMs  work,  Norman  Slssenwlne,  then  of  AFCRL,  proposed  a date-collection 
program  on  a global  scale  to  accumulate  such  data  through  routine  observation  by  craws  of  commercial 
aircraft.  The  method  was  Implemented,  and  Ivar  Lund  of  AFCRL  published  three  papers  on  the  subject 
(Refs.  2,  3,  A). 


Ref,  1.  USAF/ETAC  Technical  Report  186, 

Ref,  2.  I. A.  Lund,  1 9651,  Estimating  the  Probability  of  Clear  Llnc-of-Slght  from  Sunshine  and 
Cloud  Cover  Observations,  J,  Appl . Meteor.,  4,  p,  714-722. 

Hef.  3.  I, A.  Lund,  1966:  Methods  for  Estimating  the  Probability  of  Clear  Llnes-of  Sight,  or 
Sunshine,  through  the  Atmosphere,  J.  Appl.  Meteor.,  5,  p.  769-777, 

Ref.  4,  I. A.  Lund,  and  M.D,  Shanklln.  1972i  Photogrammetlcal 1y  Determined  Cloud-Free  Llnes-of 
Sight  Through  the  Atmosphere,  J,  Appl.  Meteor.,  TJ_,  P-  773-782. 
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Th«  cloud-free  1 Ine-of-slght  data  to  be  used  with  thl*  modal  was  prepared  by  ETAC  on  the 
basis  of  tha  AFCRL  modal. 

Tha  atmospheric  transmission  model  Initially  used  was  the  LOWTRAN  Program  of  tha  AFCRL  ereatad 
by  John  Selby,  This  program  was  redone  for  Blberman  by  Selby  to  make  batter  and  fuller  use  of  a 
broader  base  of  scattering  medals  developed  at  AFCRL  under  Robert  Fann  and  to  allow  greater  flexi- 
bility to  accept  a variety  of  Inputs  and  data  formats  eonmonly  obtained  In  meteorological  records. 
This  new  program  will  be  known  as  tha  LOWTRAN  ui. 

Overall  program  design  was  accomplished  by  deorga  du  Mels  of  IDA,  who  Is  primarily  responsible 
for  the  final  numerical  analysis  reported  herein. 

The  overall  base  of  the  model  grew  from  a long-range  program  of  studies  of  system  analysis 
and  system  synthesis  by  Rosell  of  Westlnghouse-,  Sendall,  then  of  Xerox,  under  the  USAF  AFAL  Program 
698  DF ; and  the  work  of  Blberman  and  Schnltxler  of  IDA  under  a continuing  study  of  Infrared  and 
.night  vision  technology  supported  by  the  Office  of  Research  and  Technology,  Director  of  Defense 
Research  and  Engineering,  Department  of  Dafanse. 

The  work  reported  here  Is  supported  as  part  of  an  IDA-sponsored  study  under  Its  Internal 
central  research  program. 


2.  THE  FOUNDATIONS  OF  THE  BASIC  METHOD 

In  past  work,  Rosell  (Ref.  S)  has  davaloped  and  with  Willson  hat  shown  (Refs,  6,  7,  B)  that 
the  ability  of  an  observer  to  accomplish  detection,  recognition,  or  Identification  of  a targat 
object  In  a displayed  Image  Is  best  predicted  by  the  signal -to-nolse  ratio  of  the  target  Image. 
Further,  Rosell  develops  the  concept  of  equivalent  bar  pattern. 

This  concept  goes  back  to  the  work  of  John  Johnson  (Ref.  9),  who  showed  that  an  observar  who 
could  detect,  reengntxe,  or  Identify  an  object  such  as  a man  or  a vehicle  was  also  able  to  just 
resolve  one,  four,  or  slx-and-one-half  pairs  of  black  and  white  bars  (of  a typical  photographic 
resolution  chart)  when  one,  four,  or  slx-and-one  half  pairs  of  sueh  bar  widths  fit  Inside  the 
critical  (usually  minimum)  dimension.  (See  Tables  I and  II.) 


Ref.  S. 
Ref.  6. 
Ref.  7. 

Ref.  8. 

Ref.  9. 


L.  M.  Blberman  and  S.  Nudelman  (Eds,),  Photoelectrontc  Imaging  Devices,  Planum  Press, 
Naw  York  (1971),  Vol . 1,  Physical  Processes  and  Methods  of  Analysis,  chapter  If. The 
Limiting  Resolution  of  Low-Llght-Level  Imaging  Sensors  by  F.  A.  Rosell,  pp.  307-329. 

L.  M.  Blberman  (Ed,),  Perception  of  Displayed  Information,  Plenum  Press,  New  York 
(1973),  Chapter  5,  Recent  Psychophysical  Experiments  and  the  Display  Slgnal-to-Nolsa 
Ratio  Concept,  by  F.  A.  Rosell  and  R.  H.  Willson,  pp.  17S-232. 

F.  A.  Rosell  and  R,  H.  Willson,  "Recent  Psychophysical  Experiments  and  the  Display 
Slgnal-to-Nolsa  Ratio  Concept,"  (September  7,  1972),  Report  No,  ADTM  No.  110, 
Westlnghouse  Electric  Corporation,  Defense  A Electronic  System  Center,  Systems 
Development  Division.  P.0.  Box  746,  Baltimore,  Md,  21203. 

F.  A.  Rosell  and  R,  H.  Willson,  'Performance  Synthesis  of  Electro-Optical  Sensors," 
Report  No,  AFAL-TR- 74-1 04 , Contract  No,  F3361 5-73-C-41 32 , (April  1974),  Westlnghouse 
Electronic  Corporation,  Defense  & Electronic  Systems  Center,  Systems  Development 
Olvlslon,  P.O.  Box  746,  Baltimore,  Md,  21203 

John  Johnson  (1958),  Image  Intenslfler  Symposium,  Fort  Belvolr,  Va.,  October  6-7,  1968, 
AD  220160. 
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TABLE  I.  JOHNSON'S  LEVELS  OF  OBJECT  DISCRIMINATION 


. Classification  of 
discrimination  level 

Meaning 

Detection 

An  object  Is  present 

Orientation 

The  object  is  approximately  symmetric  or  asymmetric 
and  its  orientation  may  be  discerned 

Recognition 

The  class  to  which  the  object  belongs  may  be  discerned 
(e.g,,  house,  truck,  man,  etc,) 

Identification 

The  target  can  he  described  to  the  limit  of  the  observer's 
knowledge  (e.g.,  motel,  pickup  truck,  policeman,  etc.) 

TABLE  II,  JOHNSON'S  CRITERIA  FOR  THE  RESOLUTION  REQUIRED  PER  MINOR 
OBJECT  DIMENSION  VERSUS  DISCRIMINATION  LEVEL 


Discrimination  ’aval 

Resolution  par  minor  object 
dimension,  TV  Unas 

Detection 

,♦1.0 
-0.  B 

Orientation 

2 a+0'8 
z,8-o.a 

Recognl tlon 

B O*1'6 
"-0,4 

Identl f Icatlon 

12  8+3l? 

Llttlt  mnra  was  said  about  tha  conr.apt  In  a quantitative  sense  until  RotalVa  work  rafaranead 
abova,  Hart  It  becama  claar  that  tha  work  dona  by  Johnson  ralatad  tha  number  of  Una  pairs  of  a 
bar  chart  corraspondlns  to  the  Minor  dimension  of  an  object  with  tha  ability  of  an  observer  to 
llmlnally  detect,  recognise,  or  Identify  that  object,  l.e.,  perform  that  function  at  tha  level  of 
SO  percent  probability. 

Rosell  went  further  to  establish  the  concept  of  tha  oquivalent  bar  pattern  Illustrated  In 
Fig,  1.  This  concept  creates  a bar  pattern  equivalent  In  dimensions  to  the  object  to  be  viewed. 
For  the  various  Increasingly  difficult  visual  tasks,  the  bar  pattern  Is  composed  of  Increasingly 
narrower  bars,  according  to  the  demands  of  Johnson's  criteria,  If  tha  device  In  question  produces 
an  Image  of  tha  bar  chart  that  permits  an  observer  to  "break  out"  the  Individual  bars  In  tha 
pattern  75  parcant  of  the  time,  that  observer  should  be  able  to,  for  example,  recognise  the 
reel  object  75  percent  of  the  time  If  the  equivalent  bar  pattern  waa  composed  of  four  pairs 
of  black  and  white  bars. 
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TABLE  III*.  BEST  ESTIMATE  OF  THRESHOLD  SNRn  FOR  DETECTION, 
[ RECOGNITION,  AND  IDENTIFICATION  OF  IMAGES 


1 


kj , . 

ru  14...  Threshold  SNRn  for  a Single  Bar 
O.i  ' of  Spatial  Frequency 

Discrimination  Minimum  On,  Mnes/p1cture , height>,qya.1.  to 


Level 

Backoround 

Dimension 

M 

300 

500 

Zoo 

Detection 

Uniform^ 

1 

2.8 

2.8 

2.8 

2.8 

Detection 

Clutter 

2 

4.8 

2.9 

2.5 

2.5 

Recognition 

Uniform 

a 

4.8 

2,9 

2.5 

2.5 

Recognition 

Clutter 

8 

6.4 

3.9 

3.4 

3,4 

Identification 

Uniform 

13 

5.8 

3.6 

3.0 

3.0 

®For  a vlewlng-dlstance-to-plcture-helght  ratio  of  3.5. 
^Treated  as  an  aperiodic  object. 


TABLE  Illb,  SPATIAL  FREQUENCY  AND  LINES  PER  PICTURE  HEIGHT  FOR  EQUIVALENT 
BAR  PATTERN,  TANK  DETECTION,  AS  A FUNCTION  OF  RANGE 


Dlitinci, 

kn 

Tercet 

Angle, 

mrad 

Bar 

Angle, 

JD£4i- 

Lines  Per  Picture  Height  . 

- .for  Field. of  view  lndHttrii 

2,6“  6,n»  7,8*  10.0° 

0.6 

7.8 

3,9 

9 

18 

27 

36 

1.0 

3.9 

1.96 

18 

36 

53 

71 

1.5 

2.6 

1.3 

27 

53 

80 

107 

2.0 

1.96 

0.976 

36 

71 

107 

142 

2.6 

1.66 

0.78 

46 

89 

134 

178 

3.0 

1.30 

0,65 

63 

107 

160 

214 

3.6 

1.11 

0.56 

63 

126 

189 

252 

4.0 

0.976 

0,48 

72 

146 

215 

289 

4.6 

0.867 

0,44 

79 

156 

236 

316 

6.0 

0,780 

0,39 

»B 

178 

266 

356 

6.6 

0.709 

0.35 

96 

193 

2B9 

386 

6.0 

0.650 

0,33 

106 

210 

315 

420 

6.6 

0.600 

0,30 

118 

231 

34? 

463 

7.0 

O.S67 

0.28 

124 

248 

371 

435 

7.6 

0,620 

0.26 

134 

267 

400 

533 

8.0 

0,487 

0,24 

146 

269 

433 

578 

8,6 

0.469 

0.23 

161 

302 

453 

603 

5.0 

0.433 

0.22 

158 

316 

473 

630 

9.6 

0,411 

0.21 

163 

334 

496 

661 

10.0 

0,390 

0.20 

174 

347 

521 

694 

1 
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*The  fields  of  view  show  ere  in  the  horl»ont«l  dimension,  Usually,  the 
verticil  dimension  Is  three-fourths  of  the  horliontei,  The  lines  per 
picture  height  ere  therefore  taken  across  three-fourths  of  the  horliontei 
dimension, 

For  a dlstanco  Indicated  In  Column  1,  the  front  aspect  of  the  tank  target  subtends  an  angle  In 
mllllradlans  shown  In  Column  2.  The  equivalent  bar  width  for  the  detection  model  Is  shown  In 
Column  3.  The  bar  width  In  lines  per  picture  height  for  four  fields  of  view  Is  shown  In  Columns 

!j  4,  5,  6,  and  7. 
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For  my  probability  other  than  SO  percent,  the  value  of  the  slgnal-to-nolse  ratio  need*  to  be 
adjusted  according  to  the  gausslan  form  shown  In  Fig.  2.  Here  the  value  of  50  percent  probability 
Is  related  to  a normalized  value  of  1.0,  and  other  probabilities  require  greater  or  lesser  values, 
l.e,,  90  percent  probability  corresponds  to  1.5,  and  10  percent  probability  corresponds  to  0.5. 
These  normalization  constants  would  be  applied  to  the  slgnal-to-nolte  value  for  50  percent  proba- 
bility to  achieve  the  90  percent  or  the  10  pereent  value.  If  a signal -to-nolae  ratio  of  2.8  it 
required  for  50  percent  detection,  then  1.5  x 2,8  Is  required  for  90  percent  probability  of 
detection. 


■i 


!;> 


is 


NORMALIZED  SNR. 

s-is-m.i  0 

FIGURE  2.  NORMALIZED  DISPLAY  SIGNAL-TO-NOISE  RATIO 


For  a range  of  object  sizes,  and  therefore  a range  of  spatial  frequencies,  a range  of  slgnal- 
to-nolte  ratios,  rather  than  a single  value,  Is  Indicated,  This  affect  It  due  to  the  decreasing 
efficiency  of  the  eye  at  an  Integrator  for  Images  above  a given  angular  size, 

The  details  of  this  topic  will  not  be  repeated  here  but  can  be  found  primarily  In  Chapter  5 
and  the  related  Chapter  4 of  Perception  of  Displayed  Information  (Ref.  10),  and  In  "Low  Light 
Level  Devices:  A Designers'  Manual"  (Ref,  11).  These  two  texts  give  detailed  equations,  etc. 
for  analyzing  the  performance  of  television-type  devices. 


1 


Ref.  10, 


L,  M. 

(1973),  Chapter  5,  Recent 


. Information,  Plenum  Press,  New  York, 
xpariments  and  the  Display  Signal -to- 

‘ 2.  and  Chapter  4, 
Analysis  of  Noise -Required  Contrast  and  Modulation  In  Image-Detecting  and  Display 
Systems,  pp.  119-166. 


Blberman  (Ed,),  Perception  of  Displayed 
ReeenFPsychbphysIcal  Expe 

Noise  Ratio  Concept  by  F.  A,  Roseli  and  R.  H.  Willson,  pp.  175-232,  and 


Ref.  11. 


L.  M.  Blberman  at  al,  Low-Light-Level  Devices:  A Designers'  Msnual,  Institute 
for  Defense  Analyses  Report  R-169  (August  1971), 


LtiH'diti’LW' ..i1. i v .. ■ i^i ][,l  lliU 


.L.ImilM  •'  l (i  l t'A.K. . 
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3,  CALCULATIONS  FOR  FLIRS 

For  Image- forming  Infrared  devices  such  at  forward-looking  Infrared  (FL1R)  devices,  a new  con- 
cept has  grown  up  that  Is  closely  related  to  the  material  above,  This  concept  Is  Minimum  Resolvable 
Temperature  .(MRT).  This  MRT  is  the  Incremental  temperature  between  a four-bar  pattern  and  Its 
background  necessary  to  produce  a Just-discernible  pattern  to  a viewer  on  a given  display.  Thus, 

MRT  Is  related  to  the  device  and  the  slie  and  quality  of  the  display  device  producing  the  output 
Image  (tee  Fig.  3).  This  topic  Is  covered  In  detail  In  an  Appendix  of  the'Flnal  Report  of  Study 
Group  No.  S of  the  NATO/AGARD  Aerospace  Applications  Studies  Committee  (Ref.  12),  and  In  the 
Sendai 1-Rosell  report  "E/0  Sansor  Performance  Analysis  and  Synthesis  (TV/1R  Comparison  Study)," 
fief.  13, 


k,  oyslw/mrad 

fSAStO  ON  A FIXED  CHOICE  OF  OPTICS  AND  DETECTOR  SIZE  AND  NUMIER. 


FIGURE  3.  TYPICAL  MRT  CURVES 


Ref.  12.  NATO/AGARD  Aerospace  Applications  Studies  Committee  No.  5,  "Night  Vision  Devices  for 
Fast  Combat  Aircraft,  Final  Report,"  AGARD  AR-73,  December  1975,  NATO  SECRET. 

Ref.  13.  Robert  L.  Sendai  1 , Xerox  Corporation,  F.  A.  Resell,  Meetinghouse  Electric  Corporation, 
Defense  A Electronic  Systems  Division,  "E/0  Sensor  Performance  Analysis  end  Synthesis 
(TV/IR  Comparison  Study),"  AFAL-TR-72-37A,  Final  Report  (April  1973). 
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Since  MRT  Is  that  value  of  Incremental  temperature  that  producei  a value  of  signal -to-nolse 
ratio  sufficient  to  allow  an  observer  to  break  out  the  HUT  test  bar  pattern  at  SO  percent  probe* 
blllty,  then  corresponds  to  the  normalization  factor  relating  signal -to-nolse  to  probability. 


4.  FIRST  ORDER  CORRECTIONS  TO  BAR  PATTERN  DATA 

Though  It  Is  not  the  purpose  of  this  paper  to  establish  the  theoretical  basis  of  perception  of 
displayed  Imagery,  It  Is  quite  Important  to  understand  the  means  and  the  rationale  for  using  MRT 
data  to  establish  the  probability  of  detecting,  recognizing,  or  Identifying  real  objects  rather  than 
the  standard  four-bar  MRT  test  patterns, 

In  1960,  Coltman  and  Anderson  reported  (Ref.  14)  that  the  detectability  of  a group  of  bars  was 
a function  of  the  number  of  such  bars  In  the  bar  pattern. 

To  show  the  Impact  of  reducing  the  number  of  bars  available  to  the  observer,  Coltman  and 
Anderson  devised  the  experiment  shown  In  Fig.  4.  The  displayed  psttern  was  left  fixed,  and  a series 
of  cardboard  apertures  were  employed  to  vary  the  number  of  lines  teen  by  the  observer  (Ref.  1*, 
p.  662).  The  mask  was  of  square  aspect  ratio.  The  results  ts  shown  In  Fig,  4 "show  that  the 
observer  probably  uses  no  more  than  seven  line  pairs  In  making  an  Identification.  As  the  number 
which  he  Is  permitted  to  sac  Is  dicreasad,  the  signal  required  rises  rapidly,  being  greater  by  a 
factor  of  four  when  only  one  line  pair  It  presented."  (Ref.  14,  p,  862.) 


NUMBER  OF  LINE  PAIRS  SEEN  THROUGH  MASK 


io-e-rs-4 

FIGURE  4.  NUMBER  OF  LINE  PAIRS  SEEN  THROUGH  MASK 
(ADAPTED  FROM  REF.  14) 

Schade  (Ref.  15)  also  notes  that  "the  sampling  aperture  of  the  eye  for  lines  and  edges  Is  Its 
line  Image,  limited  In  length  to  fourteen  equivalent  point  Image  diameters."  These  two  observations 

Ref.  14.  U,  W.  Coltman  and  A.  E.  Anderson  (I960).  Noise  Limitations  to  Resolving  Power  In 
Electronic  Imaging,  Proceedings  IRE,  48(5),  pp.  868-865. 

Ref.  16.  0.  H.  Schade,  Sr.  (1956),  Optical  and  Photoelectric  Analog  of  the  Eye,  Journal  Optical 
Society  of  America,  46(9),  p.  731. 
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glvt  4 possible  explanation  for  the  use  of  the  elemental  Image  of  site  1/NTV  • 1/NTv>  However, 

(f  this  Is  to  hold  over  • wide  range  of  spatial  frequencies,  It  Is  necessary  to  conclude  that,  at 
the  pattern  spacing  changes,  the  eye's  ability  to  Integrate  along  the  line  changes  In  direct  pro- 
portion, or  else  It  reaches  some  limit.  This  Is  at  considerable  variance  with  the  results  obtained 
In  the  reetangJn  experiment  of  Rosell  and  Willson,  loc  clt.  In  Rosell  and  Willson  It  was  shown 
that  the  eye  could  Integrate  a line  of  length-to-wldth  aspect  ratio  from  1:1  to  at  least  48:1  and 
perhaps  aver,  more,  since  no  and  point  was  determined. 

The  point  of  the  above  Is  that  Coltman  and  Anderson  reduced  the  number  and  lenoth  of  the  bars 
In  their  pattern  and  ascribed  the  variation  In  results  to  the  number  of  bars,  Ignoring  the  effect 
of  length,  The  work  of  Hotel  1 and  Willson  predicts  the  same  results  based  upon  variation  In  bar 
iflhPTjcg  the  number  of  bars,  and  further  show  similar  results  based  upon  as  little  as  one 
bar  of  varying  aspect  ratio. 

The  aspect  ratio  of  standard  bar  patterns  Is  either  8:1  for  the  USAF  three-bar  pattern  or  7:1 
for  the  standard  MRT  test  pattern,  (See  Fig,  1b,  le.) 

If  we  choose  a more  or  lets  square  cross  section  of  a tank  at  a target  and  place  one  line  pair 
across  that  target,  the  aspect  one  of  those  two  bars  Is  approximately  2tl  at  opposed  to  the  7:1  of 
the  MRT  bar  pattern. 

Since  the  slgnal-to-nolse  ratio  In  the  Image  of  a bar  Is  proportional  to  the  square  root  of  the 
aspect  ratio,  the  slgnal-to-nois*  ratio  In  the  2:1  aspect  bar  will  be  the  square  root  of  2/7  timet 
that  of  the  7:1  aspect  ratio  or  about  0.83.  Thus,  the  signal -to-nolse  ratio  from  a 1:2  aspect  bar 
would  require  a value  of  MRT  greater  than  that  measured  by  tha  ratio  of  1//77T. 

One  can  correct  for  any  value  of  aspect  ratio,  c,  when  MRT  data  Is  used  by  correcting  the 
value  of  MRT  at  any  spatial  frequency  by  the  factor  1/v't7T, 

Thus,  the  Important  quantity  In  predicting  the  performance  of  devices  like  FLIRs  Is  the 
value  of 


AT 

MRT//E77 


Tables  iVa  and  IVb  show  the  values  of  MRT  versus  spatial  frequency  for  the  standard  aspect  bar 
target  and  for  an  equivalent  bar  pattern  representing  the  front  aspect  of  a tank  In  which  the 
individual  bars  exhibit  an  aspect  of  approximately  2:1.  These  two  sets  of  MRT  values  are  used  to 
show  the  uncorrected  data  at  well  at  the  bar  aspect  corrected  data  In  Table  Vb, 


One  may  thus  apply  the  AT/(MRT/v'c7T)  to  establish  the  signal  to  noise  In  an  Image  and/or  tha 
probability  that  the  equivalent  bar  pattern  or  Its  counterpart  object  will  be,  say,  recognised. 

One  further  Important  approximation  Is  used  In  this  work  and  It  concerns  the  computation  of 
AT  at  the  sensor.  It  can  be  shown,  and  Indeed  Sandal!  does  demonstrate  (Reft.  12  and  13),  that 
If  a small  Incremental  temperature  AT  between  an  object  and  Its  background  It  seen  by  a tensor 
through  an  attenuating  atmosphere  with  a transmission  at  a given  range  of  tR  than  the  value  of  AT 
tensed  at  the  tensor  location  Is  an  apparent  AT  or  AT#p  equal  to  ATt  . 

The  material  above  described  the  perceptual  factors  associated  with  tha  aspect  ratio  of  a bar 
In  a bar  pattern,  the  variation  of  MRT  with  that  aspect  ratio,  and  tha  method  of  correctly  applying 
standard  MRT  data  to  nonstandard  bar  patterns  and  thus  to  the  “equivalent  bar  pattern"  stressed 
In  this  report. 
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If  we  now  refer  to  the  various  source!  of  material  on  the  value  of  signal  to  noise  required 
for  different  ylsuel  task*  such  ts  detection,  recognition,  etc.  listed  In  Table  Illl,  It  becomes 
clear  that  the  signal  to  noise  required  Is  different  for  the  various  tasks  and  further  Is  different 
for  different  sited  objects  of  Interest. 

Table  Ilia. lists  the  perceived  slgnal-to-nolse  ratio  thresholds  at  SO  percent  probability  for  a 
number  of  visual  tasks,  This  material  Is  explained  In  a number  of  publications  such  at  Reft.  6,  7, 
8,  11,  12,  and  13.  For  the  reasons  explained  In  those  references,  the  perceived  signal -to-nolse 
ratio,  sometimes  referred  to  as  the  signal -to-nolse  ratio  at  the  display,  SNRq,  or  the  slgnal-to- 
noise  ratio  In  the  displayed  Image  must  be  larger  for  largo  Images  than  for  moderately  small  ones, 

In  the  cate  of  a FlIR,  the  measured  value  of  HAT  Is  the  measured  response  of  a human  observer 
and  thus  the  effects  of  the  Integration  of  the  human  eye  are  Included  In  the  measured  response. 

In  the  cate  of  computed  values  of  MRT,  this  factor  mutt  be  Introduced  Into  the  calculation. 

In  the  cate  of  predictions  of  observer  performance  computed  from  description  of  the  scene  and 
FLIR  parameters,  such  corrections  are  also  necessary. 


8,  DETAILS  OF  THE  CALCULATION  PROCEDURE 


Bated  upon  the  factors  discussed  above  we  constructed  an  algorithm  for  the  successive  compute* 
tlon  of  atmospheric  transmission,  signal  to  noise  at  a series  of  ranges  and  thus  at  a series  of 
appropriate  spatial  frequencies. and,  finally,  the  probability  of  detection  of  a tank  for  each  range 
of  Interest.  The  successive  operations  are  listed  below. 

1.  Determine  from  tables,  records,  etc.' the  value  of  AT  to  be  expected  for  a glvun  target. 

2.  Determine  the  dimensions  of  the  target, 

3.  Decide  on  the  various  ranges  of  Interest. 

A,  Determine  temperature,  dew  point  and  visibility  for  time  and  place  of  Interest.  See 
Figs.  6a-d  for  data  over  full  year. 

5.  Using  the  LOWTRAN  III  model,  calculate  what  the  atmospheric  transmission  Is  for  the  chosen 
spectral  band  and  range. 

6.  Calculate  the  angular  subtense  of  the  minimum  dimension  of  the  object  at  each  range. 

7.  Calculate  the  "equivalent  bar  pattern"  frequency  for  the  target  at  each  range. 

8.  Calculate  the  aspect  ratio,  c,  of  the  bar  In  the  equivalent  bar  pattern. 

9.  Obtain  or  calculate  the  MRT  as  a function  of  spatial  frequency  for  the  sensor  or  sensors 
of  Interest. 


10.  For  the  ranges  of  Interest  and  therefore  for  the  "equivalent  bar  pattern"  frequency,, 
determine  the  values  of  HRT/^e7T. 

11.  Calculate  the  value  of  ATAp  by  multiplying  the  AT  of  the  object  at  iero  range  by  the 
atmospheric  transmission  at  the  range  of  Interest. 

12.  Calculate  ATAp/(MRT//c77). 


13.  Enter  the  normal  lend  probability  curva,  (See  Fig.  8.) 

14,  Note  the  probability  of  achieving  the  given  function. 
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6.  FIRST  RESULTS 

The  application  of  this  algorithm  to  real  data  results  In  the  data  tabulated  In  Table  Va  and 
Vb.  Table  Va  shows  the  probabilities  of  detection  for  a frontal  aspect  tank  at  one  hour  In  January 
and  one  hour  In  August,  uncorrected  as  well  as  corrected  for  aspect. 

Although  Table  V Is  quits  carefully  calculated,  It  should  not  be  taken  very  seriously.  The 
reasons  for  Its  lack  of  credibility  lie  In  the  wide  variations  In  that  kind  of  data  as  weather  con- 
d 1 t Ions  vary,  and  they  do  vary  widely  and  rapidly.  Although  a glance  at  the  sample  weather  data  In 
Table  VI  doesn't  teem  to  show  much  change,  examine  column  6.  It  it  clear  that  although  the  humidity 
or  temperature  doesn't  vary  vary  widely,  tie  visibility  does,  and  that  affect  alone  can  at  times 
put  FLIRs  out  of  business.  Even  though  our  study  considers  only  the  8.6-11  micrometer  and  not 
visible  band  radiation,  the  same  particulates  and  aerosols  that  cause  trouble  In  visibility  also 
cause  trouble  In  a.S-11,  though  not  with  a one-to-one  correspondence. 

Figures  6a  and  6b  show  the  hourly  visibility,  air  temperature,  dew  point,  and  the  corresponding 
computed  atmospheric  transmittance  for  January  and  August  1970.  It  Is  clear  that  If  system  perfor- 
mance Is  weather-dependent,  the  performance  will  vary  widely  and  frequently.  NOTE:  DAY  OF  MONTH 
IS  ALWAYS  ABSCISSA. 

7.  DISCUSSION  OF  RESULTS 

from  the  data  on  MRT  weather  and  the  equivalent  bar  patterns  for  given  targets,  we  have  com- 
puted and  plotted  a series  of  curves  showing  and  comparing  dats  as  a function  of  the  more  Important 
parameters  of  weather  and  sensor. 

Figures  6a  and  6b  show  the  weather  data  and  the  computed  transmittance  data  at  8.5-11  urn  for 
distances  (ranges)  of  5 and  B km.  The  choice  of  ranges  was  made  on  the  basis  of  the  exponential 
effect  of  range  upon  transmission.  The  next  of  this  series,  Figs.  7 and  8,  show  the  range  at  SO 
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TABLE  VI*. 


SAMPLE , ILLUSTRATIVE  HEATHER  DATA  FOR  JANUARY  B 
WEST  OCRMANY 


1170,  HANNOVER, 


TtV 


Air 


"WWI' 


BSQSiL 

JiL 

-mm 

Point 

Temp . 
•0 

VlilDUlty 

X* 

caver 

Lumt 

Factor 

Velocity 

(Knotei 

Direction  HUilkere 
(Decreet)  tovort  m in 

1 

oioo 

-5. 

•3, 

1.6 

a. 

5. 

5. 

60. 

1 

) 

0300 

-5. 

• 3. 

1.6 

a, 

5. 

3. 

80. 

1 

1 

0)00 

•3 « 

3. 

a. 

5. 

5. 

ao. 

1037. 

1 

1 

0300 

-5. 

-5. 

1.6 

a. 

5. 

7'. 

70. 

1 

1 

0500 

•Jl 

-5. 

1.3 

a, 

3. 

a. 

100. 

1 

1 

otoo 

•a. 

-7. 

0.3 

a. 

5. 

a. 

no. 

1017. 

1 

a 

oroo 

-a. 

•7. 

0.3 

a. 

5. 

a. 

ao. 

1 

a 

oaso 

•a. 

-7. 

0.) 

a. 

a 

a. 

ao. 

1 

a 

oaoo 

•a. 

•7, 

0.3 

a. 

5. 

9. 

90. 

iota. 

1 

a 

1000 

-a. 

•a. 

0.6 

a. 

5. 

a. 

90. 

1 

a 

1100 

-a. 

-a, 

0.6 

a. 

5, 

7. 

100. 

1 

a 

1300 

-a. 

•a. 

1.) 

a. 

5. 

10. 

ao. 

iota. 

1 

a 

1)00 

-a. 

•8 1 

1,4 

a. 

5. 

a. 

ao, 

1 

a 

1*00 

-a. 

•8 1 

3. 

7, 

3. 

a. 

ao 

1 

a 

1500 

•a. 

-7. 

i.a 

7. 

3. 

i). 

90, 

1016. 

l 

a 

1600 

•a, 

.a  * 

1.6 

7, 

). 

13. 

ao. 

1 

a 

1700 

•a, 

-7. 

i.a 

7. 

). 

13. 

100. 

1 

a 

laoo 

-a. 

-7, 

3.6 

6. 

3. 

13, 

100. 

iota. 

1 

a 

1900 

>8. 

-7. 

3.0 

a. 

3. 

15, 

no. 

1 

a 

3000 

-a, 

-7. 

1.5 

). 

3. 

15. 

no. 

1 

a 

3100 

-a. 

-7. 

).o 

6. 

3. 

17. 

90. 

1035. 

l 

a 

3300 

-a, 

-8  * 

).0 

7, 

3. 

1), 

90. 

X 

a 

3)00 

•a. 

-7. 

3.0 

7. 

), 

17. 

100. 

1 

a 

3300 

•a. 

•6 , 

3.0 

7. 

3, 

17. 

90.  , 

1033. 

TABLE  Vlb.  SAMPLE,  ILLUSTRATIVE  WEATHER  DATA  FOR  AUGUST  B,  1970,  HANNOVER, 
WEST  GERMANY 


Month 

Day 

Hour 

DtK 

Point 

*0 

"Air'  ' 
Temp, 

•C 

vuieiuty 

Km 

Diou4 

Cover 

(elchthe) 

lumb 

-Fetter 

■vra — 

Velocity 

(Knot!) 

Ulna 

Direction 

(Decree! > 

etronetrie  rreeeure 

NUllbirt 

(every  3rd  hour) 

a 

a 

0100 

16. 

17. 

6.00 

a. 

3, 

0, 

0. 

a 

8 

0?00 

15. 

16. 

6.00 

a. 

4 i 

3. 

90. 

a 

a 

0300 

16. 

16. 

5.00 

a. 

3. 

2. 

90, 

1013. 

a 

a 

0300 

15. 

16. 

5.00 

a. 

4, 

3. 

40. 

a 

a 

0500 

15. 

16. 

5.00 

7. 

4, 

2. 

90. 

a 

a 

0600 

17. 

17. 

3.00 

7. 

4, 

5. 

50. 

1013. 

a 

a 

0700 

17. 

19. 

3.00 

7. 

4. 

4, 

70, 

a 

a 

oaoo 

19. 

30, 

3.00 

6. 

3. 

6. 

60. 

a 

a 

0900 

19. 

33. 

5.00 

7. 

4. 

7. 

90. 

1013. 

a 

a 

1000 

19. 

3). 

7.00 

6, 

4, 

9. 

100c 

a 

a 

1100 

16. 

33. 

9.00 

6, 

a, 

10. 

90. 

a 

8 

\?oo 

16. 

33. 

7.00 

8, 

4. 

10, 

100. 

1011. 

a 

e 

1)00 

17. 

19. 

3.00 

a. 

6. 

a. 

330. 

a 

a 

1300 

17. 

17. 

7.00 

a. 

6. 

4. 

90. 

a 

a 

1500 

17. 

17. 

a, oo 

a. 

5. 

1. 

90. 

1013. 

a 

a 

1600 

17. 

ia. 

a,  oo 

6. 

4, 

3. 

90. 

a 

a 

1700 

17, 

la. 

6.00 

a. 

5. 

4, 

10, 

a 

a 

1800 

IV 

it, 

7.00 

7. 

3. 

3 

360. 

1013. 

a 

a 

1900 

16. 

us. 

).eo 

5. 

3. 

4. 

350, 

a 

a 

3000 

16. 

16. 

3.00 

a 

3. 

a. 

370, 

a 

a 

3100 

15, 

It. 

s.oo 

a. 

3. 

t. 

310. 

1013. 

a 

a 

3300 

16. 

It. 

3.00 

a. 

3. 

3. 

310. 

a 

3)00 

16. 

It. 

a. oo 

a. 

3. 

3. 

350 
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JANUARY  1970 


FIGURE  ba.  HOURLY  VISIBILITY.  AIR  TEWERATURE  (— ).  DEM  POIHT  (o),  AND  CORRESPONDING  COMPUTED 
ATMOSPHERIC  TRAKSMITTANCES  AT  8.5-11  im  FOR  RANGES  OF  5 AND  8 k»;  HANNOVER,  GERMANY, 
JANUARY  1970.  {VISIBILITY  CURVE  IS  TRUNCATED  AT  10  ta-  EACH  TICK  ON  ABSCISSA  MARKS 
END  OF  A DAY.) 


(VISIBILITY  CURVE  IS  TOMCATED  AT  10 


percent  probability  of  front  atpact  detection  of  a tank  In  January  1970  for  four  dlffarant  fields 
of  view  from  2-1/2  dagraas  to  10  dagraas  with  displays  of  7 and  14  Inchat.  Tha  ranga  It  graatatt 
for  tha  smallest  f laid  of  view  and  tha  largar  display  and  laatt  for  tha  largatt  flald  of  view  and 
tha  tmallar  display. 

Although  we  started  by  showing  the  data  as  ranga  at  50  pareant  probability,  wa  believe  that 
showing  probability  at  specified  ranges  Is  more  meaningful  and  tha  rest  of  the  series  of  first 
results,  l.a.,  Figs.  9 through  13,  show  tha  probability  at  various  ranges  and  conditions  with 
sufficient  comparison  data  presented  simultaneously  for  tha  major  affects  to  be  clearly  obvious 
to  tha  reader. 

Figure  9 shows  tha  probability  of  datactlon  for  a tank  seen  In  frontal  aspect  with  a 2-1/2 
degree  FOV  FL1R  at  a ranga  of  6 km,  and  It  also  shows  similar  data  for  a 7-1/2  degree  FOV  FL1R 
at  a ranga  of  4 km. 

When  the  SNRQ  or  value  of  AT/MRT,  corrected  and  uneorreeted,  Is  significantly  greater  than 
about  2.0,  the  correction  matters  little,  and  tha  probability  of  detecting  appears  to  be  100  percent 
for  any  case. 

When,  however,  OT/MRT  falls  off  appreciably  below  that  point,  the  slopes  get  Increasingly  steep 
and  the  difference  In  corrected  and  uneorreeted  data  can  be  a factor  of  five  and  more,  For  example, 
the  values  that  previously  Indicated  a probability  of  detection  of,  say,  0.63  now  become  0.12  when 
corrected. 

8.  CONCLUSIONS 

It  Is  Important  to  note  that  there  are  two  major  different  but  Important  factors  affecting  the 
detection  of  a target  by  Infrared  through  a real  atmosphere.  These  factors  apply  to  any  spectral 
band  but  here  we  consider  only  the  8.5-11  micrometer  window, 

The  first  factor  Is  gaseous  absorption  such  as  that  caused  by  atmospheric  water  or  carbon 
dioxide,  etc,  The  second  Is  that  due  to  scatter  caused  by  aerosols  and  particulate  material. 

We  all  learned  that  water  was  the  real  vllllan  In  Infrared  system  performance r and  thus  the 
wintertime  with  cold  weather  and  little  water  should  be  much  better  than  summer  with  much  humidity. 

Out  data  shows  that  the  facts  do  not  support  these  contentions,  end  that  the  base,  mists,  fogs, 
and  smokes  of  winter  are  far  more  serious,  see  Figs.  11-14.  These  figures  show  the  probability  of 
detections  1)  limited  by  the  total  atmospheric  attenuation,  2)  limited  just  by  scatter,  and 
3)  limited  Just  by  absorption.  Note  the  effects  of  scatter  tend  to  cause  the  most  severe  degrada- 
tion, l.e.,  the  calculated  transmissions  of  Figs.  6e  and  6f  bob  up  and  down  rot  with  humidity  but 
with  computations  scaled  to  visibility. 

One  may  question  the  wisdom  of  scaling  from  data  taken  at  0.55  micrometers  to  10  micrometers. 

We  worried  about  that,  Therefore  we  asked  Selby  and  Fenn  of  the  Air  Force  Cambridge  Research 
Laboratories  to  examine  the  question  In  more  detail.  The  response  follows, 

John  Selby  has  compared  the  new  extinction  coefficient  provided  by  Bob  Fenn  for  rural  and  urban 
aerosols  with  those  currently  In  LOWTRAN  III  (average  continental  aerosol).  He  found*  that  although 
the  three  models  exhibit  small  differences  In  their  spectral  characteristics  In  the  7-14  micrometer 

Personal  correspondence,  24  February  1975, 
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region,  the  average  broadband  transmittance  appears  to  be  the  same.  Each  of  these  models  In  Tables 
Vlla-c  In  fact  contains  a different  size  distribution  and  particle  number  density:  urban  - 
2.3  x 10s  parti  cl es/cm3i  rural  - 2.7  x TO6;  average  continental  • 2.8  x 103, 

The  differences  between  these  models  as  a function  of  wavelength  result  from  the  different 
particle  size  distributions  and  refractive  Indices.  For  example,  the  maritime  model  has  more  large 
particles  than  the  others,  tn  order  to  use  these  curves,  the  visual  range  must  be  given,  and  the 
data  all  compared  at  0,5  micrometers. 

Although  the  three  models  considered  above  give  the  same  average  extinction  coefficient  In  the 
8.5-11  micrometer  region,  these  models  differ  more  significantly  at  other  wavelengths.  These  will 
be  examined  (especially  the  3. 4-4, 2 micrometer  band)  In  another  paper. 

It  should  be  clear  from  that  data  that  the  model  Is  not  very  sensitive  to  variations  In  particle 
size  distribution,  and  thus  we  believe  quite  strongly  that  our  results,  though  possibly  not  at  com- 
plete as  might  be  hoped  for,  do  show  the  effects  we  wished  to  demonstrate. 

Later  on  In  1976-77,  NATO  Project  OPAQUE  will  provide  more  experimental  data  against  which 
we  can  test  our  results,  In  the  Interim,  we  present  these  findings  above  to  show  as  best  we 
can  the  effects  of  weather. 

NOTE  IN  PROOF 

Since  this  paper  was  first  presented,  detailed  studies  by  Robert  Roberts  and  Luclen  Bibarman 
(Ref.  16)  have  lead  to  the  discovery  that  the  LOWTRAN  II  and  III  models  were  erroneous  In  the 
values  used  for  the  8-14  micron  region  of  the  water  vapor  continuum.  These  corrections  have  since 
been  applied  to  the  figures  presented  herein.  Tha  results  for  August,  when  humidity  Is  significant, 
are  considerably  Improved  over  prior  calculations  using  the  LOWTRAN  II  and  III  models.  Figure  14 
shows  a comparison  of  atmospheric  transmission  by  LOWTRAN  II  and  III  and  our  corrected  calculations. 


Ref.  16.  R.E.  Roberts,  L.M.  Blberman  (IDA)  and  J.E.A.  Selby  (AFCRL) , "Infrared  Continuum 

Absorption  by  Atmospheric  Water  Vapor  In  the  8-12  pm  Window,"  IDA  P-1184,  In  press. 


30 


- n\  * :V‘  


"■  -"'.  'il niu.i.t.fU  1a\i  ■■  *.il  |»(|  j,; ... ,, . , h -,n 


i'.iA -.'-v.'-c.'  t'.t.'iW,' 


1 





;i 

;•  i 


f! 


COPY  AVAILABLE  TO  DOC  DOES  NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 


TABLE  V tin,  b,  C. 


TRANSMITTANCE  AS  A FUNCTION  OF  VISIBILITY  AND 
RANGE  (KM), CONTINENTAL  ( LOUTRAN  III) 


7l<  7'MANJf  ■ .I  "■  “,5 


.5 

.41475 

,5)909 

1.0 

.45644 

, 74733 

2.6 

.47502 

,54114 

1.0 

.16327 

,92524 

6.0 

.44176 

« 44261 

5.6 

.16077 

,45321 

6.0 

•46202 

,966)4 

7.0 

.46266 

.96343 

0.0 

•46345 

.964  33 

9.0 

.16342 

.97237 

10.0 

•64421 

.97471 

11.0 

.46447 

,97673 

ii.o 

.16461 

,97344 

11.0 

.16465 

.97435 

u.o 

.64644 

.92106 

14,0 

.46513 

,4621 1 

* 0,0 

.61576 

.44371 

25.  o 

.16527 

,4633ft 

30.0 

•»6i94 

,44943 

35. 0 

.16546 

.44051 

4U®AL 

.? 

• H 336 

,41264 

1.0 

.13573 

.79 (42 

2.0 

.iT7b7 

.92461 

1.0 

.125Q5 

, 4 4'  4 1 1 7 

4.0 

• 12064 

,94171 

1.0 

.16043 

.4:250 

6.0 

.66191 

,95i76 

T.n 

• 11260 

,96444 

o.o 

.49341 

,9hrt9l 

0.0 

.49333 

.47(94 

10,0 

. 44414 

,67445 

11.0 

. 46443 

.97447 

12.0 

.44446 

,47.9  15 

11,0 

. 44413 

.47433 

U.O 

.44497 

.46031 

10,0 

• 44564 

,94(34 

20.0 

.14575 

.43541 

21.0 

.14547 

,93  71ft 

10.0 

. l434i 

, 9ai  1 7 

31.0 

. 46542 

,99044 

049*4 

. 5 

. 61507 

,4  >237 

1.0 

. liftb* 

,74759 

2.0 

•1761* 

,36(23 

1.0 

. 62515 

,9*349 

4.0 

• 11336 

,44222 

3.0 

• 19034 

.95.1)2 

6.0 

. 14206 

,96049 

T.O 

.44291 

,94641 

6.0 

.19344 

,96447 

9.0 

.49392 

,97*47 

10.0 

.49420 

.97439 

11.0 

.14.,  50 

,67697 

12.0 

.19471 

.9795* 

11,0 

.4445 7 

.9749* 

U.O 

.44501 

. *1111 

15,0 

.69512 

.96*17 

20.0 

.44377 

.9H591 

21,0  , 

. 1 1455.3  . 

_ _.43N|lj, 

30.0 

. 14995 

.466,1 

35.0 

.44599 

.44956 

r.T" 

• Hi  1»4 
•6344s 
.79.  li 
..85  (ji 
.39*1  ( 

• 41  <>*4 

.4*345 

.lUSl 

i<Mll 

,i44|3 

• l3  (90 
.46774 

• fe  i 0» 

»i6  17* 

• 4AM? 
.46<2) 

.1?lUT 
.47)11 
.19*73 
< 1648ft 


9.0 

10.0 

15,0 

20.0 

a 0 Hi  6 

.00061 

ToTooi 

;o6oo~ 

.1(663 

.01652 

,00205 

,0(70*0 

.3)965 

.11676 

.0*»10 

.01090 

u’ojh. 

-LL’io 

.-.JJtSl 

.0017* 

• 4 J > 97 

.13756 

.19077 

712072“ 

.6. *69 

.4)660 

,27510 

.10*20 

.66(96 

.45363 

.340)5 

.2*072 

. 7*6*5 

.53059 

,39005 

.29*1* 

.7586) 

.36024 

.**4*1 

.3*207 

.79*59 
. 9 . *1  5 

- jM467 
.4*751 

-:5tW~ 

.6(351 

.67)97 

, 53590 

.609)7 

.29*56 

.69667 

.554*7 

.4903* 

,6.441 

.71669 

,00440 

.5100* 

.35491 

.73461 

.03136 

,56*07 

. . . 7601  5,. 

.507)0 

,9i675 

.40769 

.72793 

>000)5 

,4.706 

. 34461 

.7773? 

.71040 

.4  (979 

,36966 

.513(5 

.700)0 

,4.074 

.99456 

.939J7 

.79305 

• 4ft  (2( 

, S ( '( 39 

>004 1 1 

. .00000 

.00000 

* 5 ) ( 4 / 

.l;.6*4 

.01179 

.00141 

.00016 

. 7 9u  7 ( 

.)<  32 

,106b* 

,01593  ' 

.012*1. 

.9372* 

.6, '.>43 

.226*1 

.10770 

. .05*06 

».1.4.,')i 

.3?1?S 

. . .12665 

oau?  . 

• »Mi6  . 

.11  >5.9 

.6)446 

.60310 

,25*45 

.16920 

.4*474 

, ftc  64 3 

.47431 

,32940 

.22777 

, 4 1 3 J 

.5*316 

.,39557 

.20192 

, 44 • 44 

.75(64 

.57253 

.4)497 

.3)075 

,4414* 

,7(663 

.6097) 

.47730 

.37469 

, 41  (2  1 

,74176 

.641)9 

..  UU07 

,46 71b 

,31547 

.06615 

.5*792 

,66939 

, 9614ft 

.9*964 

.69150 

,57559 

,48114 

.14)2/ 

.9  * 1 91 

.71169 

.60*31 

.50971 

.46364 

.11*42 

.72936 

,5*5(9 

.53570 

,14777 

.95171 

.74*30 

,64571 

.559*6 

a 4 7^1  4 

.ftHQft 

.604.4* 

7*1,16,  . 

...  950J  i 

.47134 

.4(574 

.34111 

,77405 

.71161 

• 46^ jft 

,96471 

«6477( 

,30946 

. 75614 

.49 .74 

.9)931 

.3967* 

, 4356C. 

• 70957 

. •*  l > 3 7 ' 

•0.(41 

. 0001  6 

,00000 

.00000 

,n3i2i 

,l'i407 

•0|?5) 

,001*4 

,170020 

, 7 9l 3ft 

, V.’  4 4 1 

.11694 

,03785 

.01305 

.36iii 

,47743 

.*3067 

, t ’ *07 

.05476 

,19 <44 

,57(61 

. *1?9* 

,1415* 

.11*44 

.91-57 

,6-*2?0 

.-16*5 

. 254(7 

.17491 

.4*1(9 

.6  <146 

.46129 

,3)543 

.23416 

. 4 1 >.  0 1 

.72431 

.63444 

,34*76 

,28457 

, 1 4 15* 

.7*991 

.5739* 

,44*04 

. 1)794 

.949  17 

• 79e J5 

. hi *74 

.69475 

• 16186 

.45-07 

.30*17 

.64466 

,52197 

.4*122 

U5J  »5_  . 

-♦4U.3P 

-»?S6Jt 

0)54,4  5 

.46119 

• 3 J* 7 4 

•69660 

,36330 

• 46906 

.46(7) 

,N’*6T1 

.7J674 

.60839 

.51537 

.454)0 

,4*514 

.73446 

.63105 

.54232 

.46335 

.*5427 

.75019 

.6*1)7 

. - .55*50 

.47657 

.64705 

.60793 

,7*794 

. 69579 

_*VLH2.  . 

-.J.665JA.— 

.,>780* 

. 7(6*5 

• 7)l6» 

.9)|(11 

.37*14 

,91  (411 

.76039 

, 4««, 96 

. 4313  > 

.9)446 

.743*3 

31 


OUT  IS 


FIGURE  14.  TRANSMISSION  III  8.5-11  m BAND  ACCORDING  TO  LOUT RAN  II  AND  III  (NARKED  ORIGINAL) 
AND  TO  ROBERTS  (NARKED  REVISED)  VERSUS  DAT  OF  MONTH. 
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SUMMARY  OF  PRESENTATION  AT  WORKSHOP  ON 
ATMOSPHERIC  TRANSMISSION  MODELING 

The  presentation  concerns  the  transmission  of  aircraft  or 
rocket  plume  radiation  through  the  atmosphere.  Cambridge 
Research  Laboratories  have  developed  line  by  line  transmission 
calculations  that  are  very  useful.  In  principle  these  cal- 
culations are  applicable  for  a continuum  emitter  or  a spectral 
line  emitter  provided  every  single  emitting  rotational  line  is 
known.  This  follows  from  the  fact  that  emitting  and  absorbing 
rotational  lines  are  not  always  identical  in  wavelength  position 
or  width  of  the  lines.  We  will  now  discuss  the  limitations  of 
the  knowledge  of  the  detailed  emission  spectrum  of  plumes  and 
restrict  ourselves  mainly  to  the  4.3m  C0g  band. 

At  ambient  temperature  and  small  optical  depth  the 
4.3m  C02  bancl  consists  essentially  of  the  001-‘  000  (or  CRL 
notation  00011-* 00001)  transition.  This  band  has  a R-branch 
that  extends  to  shorter  wavelengths  and  converges  towards  a 
hand  head.  The  P-branch  extends  towards  longer  wavelength  and 
the  spacing  between  the  lines  inorease  slowly.  Higher 
rotational  lines  are  not  sufficiently  included  in  the  CRL 
computer  tape  and  the  position  of  the  presently  listed  higher 
rotational  lines  is  not  sufficiently  accurate  to  decide  on  the 
coincidence  of  lines  (such  as  emission  and  absorption  lines). 

An  example  of  this  can  be  seen  in  fig.  1 which  is  taken  from 
recent  results  by  0.  Lindquist  of  ERIM  under  an  ARPA  oontraot. 
One  sees  the  missing  R-branch  lines  near  the  band  head  and  the 
mismatch  between  calculated  and  experimental  line  positions  due 
to  Inaccuracies  in  the  rotational  constants. 

The  CRL  line  by  line  calculation  contains  about  25  transi- 
tions between  upper  vibrational  levelB  (the  so-called  hot  bands) 
with  the  upper  levels  going  up  to  about  6000  om”^.  These  hot 
bands  become  important  under  two  circumstances.  The  first  is 


connected  with  transmission  over  large  absorbing  path  lengths. 
This  causes  the  fundamental  001  -*•  000  band  to  become  strongly 
saturated  whereas  the  weaker  upper  transition  lines  grow  with 
path  length  until  even  very  weak  lines  absorb  strongly.  The 
second  arises  from  elevated  temperatures.  At  plume  temperatures 
high  vibrational  levels  are  excited  and  all  transitions  leading 
to  a reduction  of  one  vibrational  level  in  the  levels  cause 
emission  at  or  near  4.3m*  Differences  in  vibrational  levels 
decrease  in  magnitude  for  higher  levels  as  the  levels  converge 
towards  the  dissociation  energy  of  the  molecule.  Thus  all 
(or  nearly  all)  upper  state  transitions  have  their  origin  at 
longer  wavelengths  than  the  001-*000  transition.  This  in  turn 
is  the  reason  why  the  band  head  of  the  R-branch  is  unperturbed 
by  the  hot  bands  (fig.  1). 

The  P-branoh  region  at  wavelength  > 4.3m  behaves  differently 
and  all  the  various  upper  Btate  transitions  (hot  bands)  overlap 
and  mask  the  P-branch  of  the  fundamental.  The  full  half  width 
at  sealevel  of  a C02  line  is  about  0.14  cm-1.  As  lines  are  on 
the  average  2 cm~^  apart,  it  takes  about  14  lines  on  the  average 
to  fill  in  the  spaces  between  the  lines.  Thus  the  P-branoh  of 
the  4.3m  C02  band  Bhould  become  quite  continuous  at  elevated 
temperatures.  Figure  2 shows  a line  by  line  calculation  by 
J.  Selby,  CRL  of  part  of  the  P-branch  for  a 1000°K  emitter  for 
thin  optical  conditions.  For  optioal  thicknesses  as  represented 
in  plumes  the  spectrum  will  become  less  spiky  aB  weaker  lines 
increase  whereas  stronger  lines  become  limited  by  the  black-body 
condition.  One  sees  that  for  all  practical  conditions  the  region 
has  become  a pseudo-continuous  emitter  (although  of  rather  rugged 
appearance).  For  relatively  short  absorbing  path  length  (say 
of  the  order  of  a kilometer)  the  atmospheric  C02  will  absorb 
mainly  through  its  000  -*001  and  010-011  bands  (including  the 
isotop  ^COg.  Thus  the  pseudo  continuum  emitter  is  absorbed 
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by  relatively  few  lines  and  transmission  losses  are  less  than 
if  a line  by  line  correlation  between  emitting  and  absorbing 
lines  would  pertain.  Thus  transmission  depends  on  minor  spectral 
details  of  the  emitting  plume  and  depends  on  the  temperature  as 
well  as  optical  depth  of  the  emitter.  There  is  little  value  in 
measuring  atmospheric  transmission  experimentally  unless  one  oan 
simulate  the  structure  of  the  emitter  i.e.,  a near  oontinuum  in 
the  P-branch  and  a line  structure  in  the  R-bramoh. 

A further  example  will  be  shown  for  the  CO  moleoule. 

Figure  3,  taken  from  preliminary  data  by  Aerodyne,  Inc.  shows 
the  emission  of  the  fundamental  band  centered  at  4.7m  for  a 
2000°K  hot  emitter  of  2.8  meter  total  thickness  with  5 mole  % 

CO  present.  Many  of  the  1-0  and  2-1  lines  become  optioally 
thlok.  As  the  atmosphere  contains  a small  amount  of  CO  one 
will  expect  that  1-0  lines  will  become  partially  absorbed  after 
an  absorbing  path  length  of  a few  kilometer.  It  will  be  the 
2-1,  3-2,  4-3,  5-4  and  even  6-5  transition  lines  that  will  be 
transmitted  and  although  in  this  case  the  line  density  is  not 
enough  to  render  the  spectrum  pseudo-continuous  nevertheless 
it  will  be  of  sufficient  complexity  to  render  the  transmission 
problem  quite  complex. 
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ABSORPTANCE 


WAVE  NUMBER  (cm-1) 

FIGURE  la.  Absorptance  as  Measured  Near  the  CO?  R-Branch  Band  Head 
In  a CO-Og  Flame  at  0.3  cm  Optical  Depth  of  2500°  K 
Temperature.  The  lines  can  be  seen  to  be  well  separated 
ami  undisturbed  by  hot  bands,  (from  G.  Lindquist,  ERIM) 

1 IGURE  lb.  CO2  R Branch  as  Calculated  by  CRL  Line-By-Line  Code.  Lines 
near  the  band  head  are  seen  to  be  missing  and  the  last  line 
near  2397.2  cm*'  coincides  with  a minimum  of  the  measured 
spectrum.  Thus  Its  position  Is  not  sufficiently  correct  for 
detailed  emission-absorption  calculations,  (from  G.  Lindquist, 
ERIM) 
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Vincent  J.  Corcoran 
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I would  like  to  point  out  a couple  of  problem  areas  associated 
with  using  the  high  resolution  codes  to  calculate  the  atmospheric 
transmittance  for  laser  systems. 

The  first  point  is  illustrated  in  Fig.  1 which  shows  a table 
taken  from  the  December  issue  of  Applied  Optics  in  which  a comparison 
between  measured  values  of  the  absorption  coefficient  for  DF  laser 
lines  is  compared  to  the  values  calculated  by  the  group  at  AFCRL. 

The  two  columns  on  the  right  show  that  there  is  a difference  in  the 
calculated  and  measured  values  by  as  much  as  a factor  of  two.  There 
are  two  DF  laser  lines,  the  P1(8)  and  P2(7),  that  are  only  7 wave- 
numbers  apart.  Now  if  a person  is  off  by  a factor  of  two  in  the  one 
way  transmission  and  the  receiver  is  located  at  the  same  point  as  the 
transmitter,  then  the  prediction  of  the  system  capability  will  be  off 
by  a factor  of  4.  That  may  not  sound  3.ike  much  to  someone  who  only 
wants  to  get  a ballpark  feeling  for  the  system  performance  but  for 
a soldier  who  expects  to  carry  a 10  lb  pack  and  ends  up  with  a 40  lb 
pack  because  of  a factor  of  4 uncertainty  the  difference  can  be  sig- 
nificant. 

The  other  point  I want  to  make  is  illustrated  in  Fig.  2,  which 
is  a plot  of  the  transmission  vs  visibility  over  12  km  horizontal 
paths  at  sea  level  for  10.6  micron  radiation  using  the  line  by  line 
program.  It  shows  that  for  a subarctic  winter  where  there  is  little 
water  vapor  the  transmission  can  approach  one  for  a clear  day;  however 
for  a tropical  atmosphere  where  there  is  a large  water  vapor  content, 
the  transmission  can  be  down  by  a factor  of  LOO  for  the  same  visibility. 
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Incidentally  if  the  curve  is  plotted  for  a 100%  humidity  condition, 
you  obtain  virtually  the  same  curve  as  for  the  tropical  atmosphere. 
You  know  that  if  you  do  a calculation  that  is  off  by  a factor  of  300 
you  are  not  going  to  be  believed  for  very  long. 

The  third  figure  shows  a photograph  of  a scene}  under  clear 
conditions  and  foggy  conditions.  It  also  indicates  that  the  in- 
creased propagation  loss  for  an  11  mile  path  was  only  4 db.  This 
points  out  an  apparent  discrepancy  between  the  calculations  made  by 
the  AFCRL  program  and  the  experimental  measurements  since  the  fog, 
in  addition  to  being  a high  humidity  condition,  has  water  droplets 
that  would  further  increase  the  attenuation.  I know  that  there  is 
a partial  explanation  of  this  discrepancy  based  on  particle  size, 
and  I would  like  to  see  this  problem  addressed  during  the  workshop. 

The  fourth  figure  summarizes  the  two  points  that  I wanted  to 
make.  First,  the  calculations  of  any  model  must  provide  accurate 
transmission  vs  wavelength,  and  second,  the  transmission  at  any 
wavelength  must  be  accurately  predictable  as  the  meteorological 
conditions  are  changed. 
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EXPERIMENTAL  ABSORPTION  COEFFICIENTS  FOR  SEVENTEEN  DF  LASER  LINES 

LXP  compared  with  theoretical  predictions4 
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Atmospheric  Transmission  of  10. 6y  Radiation 
Over  12  km  Horizontal  Path  at  Sea  Level 
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TRANSMITTANCE  MODELING 


A general  overview  of  the  AFCRL  transmittance  modeling  activity  will  be 
given.  The  first  slide  describes  the  tools  that  have  been  developed  for 
dealing  with  this  problem: 

A)  LINE  LIST  - A compilation  of  Atmospheric  Absorption  Line  Parameters 
ha«  been  developed  (McClatchey  et  al.,  1973)  for  dealing  with  the  problem 
of  atmospheric  transmittance  modeling.  In  order  to  apply  this  line  list, 

It  Is  first  necessary  to  define  the  atmosphere  In  terms  of  molecular 
abundances,  pressure  and  temperature.  In  practice  there  are  many  applica- 
tions where  uncertainties  In  atmospheric  models  are  larger  than  the  un- 
certainties In  the  molecular  parameters  In  terms  of  transmittance  calcula- 
tions. There  are  exceptions,  particularly  as  concerns  laser  (monochromatic) 
transmittance  calculations  In  certain  spectral  regions. 

B)  HITRAN  - This  Is  a computational  technique  (computer  code)  capable  of 
using  the  line  list  as  Input  data  together  with  an  appropriate  model  atmo- 
sphere and  creaelng  a synthetic  spectrum  through  repeated  monochromatic 
calculations.  It  Is  very  time-consuming  on  the  computer,  but  provides  the 
only  method  for  synthetically  generating  high  resolution  transmittance 
spectra  (spectral  resolution  between  0 and  5 cm'^). 

C)  BAND  MODELS  - Band  models  have  been  used  In  the  past  In  lieu  of  detailed 
molecular  spectroscopic  data  to  provide  low  spectral  resolution  trans- 
mittance models.  They  have  fundamental  limitations,  because  they  depend 

on  a variety  of  aceumptlons  regarding  the  spectral  distribution  and 
Intensity  distribution  of  molecular  absorption  lines.  Perhaps  more 
accurate  band  models  can  be  developed  using  the  line  list  as  Input  data. 

In  this  way,  models  could  be  built  based  on  the  real  spectral  and  Intensity 
statistics  of  a particular  absorption  band  belonging  to  a particular 
molecular  species. 


0)  LOWTRAN  - This  Is  a particular  low  spectral  resolution  (20  cm”^)  trans- 
mittance model  (see  McClatchey  et  al.,  1972).  It  Is  a one-parameter  model 
in  which  an  attempt  has  been  made  to  optimize  the  transmittance  accuracy 
In  terms  of  the  most  appropriate  single  parameter  for  a given  molecular 
species.  LOWTRAN  was  built  as  a computer  program  (Selby  and  McClatchey, 
1972)  In  order  to  address  a wide  variety  of  systems  applications  In  which 
atmospheric  transmittance  plays  an  Important  role  and  also  in  order  to 
provide  a very  rapid  and  efficient  computational  capability.  John  Selby 
will  describe  some  recent  modifications  to  the  LOWTRAN  code  Included  In  a 
new  version  of  the  code,  LOWTRAN  3.  The  modifications  are  primarily  In  the 
2.7-ym  water  vapor  absorption,  the  Introduction  of  a revised  aerosol  model 
and  more  flexible  Input  requirements  In  terms  of  meteorological  parameters 
and  atmospheric  models. 

The  second  slide  Is  a diagram  showing  the  flow  of  information  from  the 
generation  of  the  fundamental  spectroscopic  data  through  the  generation  of 
transmittance  and  background  codes  to  be  used  In  connection  with  system 
design  and  operation.  An  Important  point  Indicated  here  Is  the  feedback 
loop  Indicating  that  comparisons  of  transmittance  models  with  field 
measurements  Is  made  and  may  lead  to  appropriate  modifications  In  the 
developed  codes.  Scattering  data  are  also  provided  as  Input. 

The  third  slide  Indicates  the  molecules  Included  In  the  ARCRL  Atmospheric 
Absorption  Line  Parameters  Compilation  and  the  number  of  molecular  absorp- 
tion lines  currently  Included  In  the  compilation.  It  should  be  noted  that 
the  lines  Included  In  this  compilation  are  based  on  their  importance  on  the 
atmospheric  transmittance  problem  and  therefore,  are  not  necessarily 
sufficient  for  use  In  dealing  with  the  emission  of  hot  plumes  containing 
the  same  molecular  species.  Some  special  studies  have  been  carried  out 
which  Indicate  that  the  line  data  currently  available  (As  of  April  1975) 
can  be  used  with  some  confidence  up  to  temperatures  around  1000K,  but  that 
their  use  as  higher  temperatures  will  require  additional  Investigation  and 
probably  Inclusion  of  additional  hot  lines. 
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The  fourth  slide  shows  a series  of  solar  spectra  observed  from  different 
altitudes  (and  displaced  by  ZQ%  with  respect  to  each  other).  These  spectra 
result  from  a contractual  program  with  the  Denver  Research  Institute.  The 
point  to  be  made  Is  simply  the  complexity  of  the  atmospheric  transmittance 
and  therefore,  the  need  to  obtain  Information  of  the  many  absorption  lines 
Involved  In  accurate  transmittance  calculations  are  to  be  performed. 

The  fifth  slide  shows  the  equations  that  relate  the  line  parameters  to  the 
monochromatic  absorption  coefficient.  It  Indicates  the  need  for  the  follow- 
ing four  fundamental  parameters:  line  position  or  frequency,  v,  line  In- 
tensity, S,  line  half-width,  o,  and  the  energy  of  the  lower  state  of  the 
transition,  E . 

The  sixth  slide  provides  the  fundamental  equations  associated  with  the 
line-by-line  (HITRAN)  computational  technique.  Equation  1 defines  the 
monochromatic  transmittance;  Equation  2 Indicates  the  necessity  of  summing 
over  all  1 lines  belonging  to  all  j molecular  species  In  order  to  determine 
the  monochromatic  transmittance  at  the  frequency,  v;  Equation  3 Indicates 
the  necessity  of  Integrating  through  a nonhomogeneous  atmospheric  path  (In 
which  pressure,  temperature  and  molecular  abundances  are  not  uniformly 
distributed);  Equation  4 Indicates  the  usual  requirement  (except  for  laser; 
transmittance  calculations)  of  performing  a convolution  of  the  monochromatic 
transmlttances  with  an  appropriate  slit  function  (or  filter  function). 

The  seventh  slide  shows  the  results  of  a HITRAN  calculation  corresponding 
to  a 10-km  horizontal  path  at  an  altitude  of  12  km  In  the  4.8-vim  region. 
Similar  plots  covering  the  entire  spectral  region  from  0.76  ym  to  31.25  ym 
have  been  published  by  McClatchey  & Selby,  1974.  Several  additional  reports 
on  laser  transmittance  are  also  referenced. 

Slide  Number  8 demonstrates  the  potential  applicability  of  band  model 
techniques  through  use  of  the  molecular  line  parameters  by  comparing  the 
results  with  a line-by-line  calculation  In  the  15-ym  spectral  region. 


Slide  Humber  9 shows  a similar  band  model  comparison  In  the  9,6-ym  band  of 


ozone. 


Slide  Number  10  describes  the  general  applications  of  this  transmittance 
modeling  activity,  In  addition  to  the  obvious  laser  transmittance  and  low 
resolution  transmittance  modeling  capability,  the  same  tools  can  be  used 
with  slight  modification  to  address  the  problem  of  atmospheric  emission 
modeling.  With  the  reservations  and  limitations  expressed  above,  the  same 
data  and  techniques  can  also  be  used  for  exhaust  plume  modeling. 


Slide  Number  11  compares  an  emission  calculation  with  measurements  made  at 
about  0.5  cm"1  spectral  resolution.  It  Is  likely  that  any  discrepancies 
Indicated  here  are  more  likely  a measure  of  the  uncertainty  of  atmospheric 
temperature  and  molecular  abundance  than  an  Indication  of  uncertainties  In 
the  molecular  absorption  parameters. 


Slide  Number  12  shows  a similar  emission  comparison  In  the  16-30-ym  region 
where  the  measurement  was  made  from  a balloon  platform  at  an  altitude  of 
70.2  kft. 


Slide  Number  13  shows  the  results  of  a general  atmospheric  emission  survey 
computed  with  the  molecular  line  list  and  degraded  In  spectral  resolution 
to  5 cm"  . Additional  background  calculations  are  provided  In  the  report 
by  Garlng  et  al.,  1972  (Secret  Report). 

Slide  Number  U shows  the  results  of  applying  the  molecular  line  data  to 
the  emission  of  a hot  (6UCK)  source  as  viewed  through  a 2-km  horizontal 
atmospheric  path  at  sea  level.  The  overall  spectral  feature  Is  the  "blue 
spike"  ai.d  at  high  (Infinite)  resolution  the  Individual  emission  lines  can 
be  seen  on  the  short  wavelength  side  and  the  individual  atmospheric  absorp- 
tion lines  can  be  seen  on  tne  long  wavelength  side. 


TRANSMITTANCE  MODELLING 


. LINE  LIST 
HITRAN 

. BAND  MODELS 
. LOWTRAN 

SLIDE  No.  1. 


Laboratory  and  Theoretical 
Molecular  Spectroscopic  Parameters 

Data  Bank  of  Spectral  Line  Parameters 

Optlcal/IR  Atmospheric  Transmission 
Computer  Code 

Comparison  with  Field  Measurements  ■ 

♦ 

Atmospheric  Transmittance  and  Background 
System  Design  and  Operation 


— - Atmospheric  Model 
— Scattering  Data 

anrl  Thancu 


SLIDE  No.  2. 


55 


TRANSMITTANCE  TRANSMITTANCE 


Table  I 

MOLECULES  INCLUDED  IN  COMPILATION 


Molecule 

Abundance  (ppm) 

Variable  (3xl024  molecules/cm2) 

No.  Entries 

h2o 

38,145 

co2 

330 

32,839 

°3 

Variable  (IxlO20  molecules/cm2) 

19,328 

N20 

0,28 

14, 969 

CO 

0, 075 

354 

ch4 

1.6 

1,741 

•A 

o 

yc 

CVJ 

490 

SLIDE  NO.  3. 

.INI  KAMI)  A I MOSh’lll  HIC  I HANSMISSION 


RAVI  LIN  9 TH  <M  ) 


WAVILINflTH  ( h I 


SLIDE  No.  4. 


56 


SLIDE  No.  5. 
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I Introduction 


In  my  experience  the  users  of  atmospheric  transmittance  computations 
can  be  separated  into  two  quite  distinct  groups:  the  system  planners  and 
the  data  gatherers.  Those  of  us  who  are  expected  to  give  technical  advice 
concerning  atmospheric  transmittance  should  be  aware  of  the  different 
types  of  information  required  by  members  of  these  two  groups. 

The  system  planner  is  concerned  with  the  conditions  which  will  limit 
the  performance  of  his  system)  not  the  ability  to  know  the  transmittance  with 
great  accuracy  at  any  particular  time  or  place.  The  planner  needs  parametric 
studies  based  on  climatic  data  to  tell  him  the  limits  of  transmittance  his 
system  is  likely  to  encounter  under  operational  conditions.  If  inaccuracies 
in  transmittance  codes  do  not  unduly  bias  results,  then  some  uncertainty 
in  transmittance  computation  is  acceptable  in  view  of  the  uncertainty  in- 
herent in  the  use  of  climatological  information. 

The  data  gatherer,  on  the  other  hand,  wishes  to  extract  all  the  in- 
formation possible  from  observations  made  through  the  atmosphere.  The 
data  gatherer  cares  little  for  climatic  data,  but  requires  the  best  possible 
meteorological  data  for  the  particular  event  of  Interest  to  him.  He  also 
requires  the  best  possible  transmittance  codes  consistent  with  the  accuracy 
of  the  data  available  to  him. 

The  accuracy  of  computed  atmospheric  transmittance  is  limited  by  a 
number  of  factors,  Applicable  meteorological  data  may  not  be  available 
for  the  time  and  place  of  particular  interest.  This  is  particularly  true  of 
ocean  areas  viewed  by  space  based  sensors.  Even  if  meteorological  data 
are  available  at  the  time  and  place  of  interest^ the  data  itself  may  be  of 
limited  accuracy.  As  others  at  this  conference  have  pointed  out,  there  is 
little  information  in  the  routinely  recorded  meteorological  data  to  permit 
the  estimation  of  infrared  attenuation  due  to  the  particles  in  the  atmosphere. 

If,  in  the  spectral  region  of  interest,  water  is  an  important  absorber^ 
the  meteorological  data  may  still  be  inadequate  even  through  radiosonde 
data  including  humidity  observations  are  available.  For  example  the 
hygrometers  on  radiosondes  are  useless  when  the  ambient  temperature  is 
below  -40  C,  a condition  which  usually  occurs  within  10  km  above  the 

♦This  work  was  sponsored  by  the  Defense  Advanced  Research  Projects  Agency 
Order  No.  2843  and  SAMSO  Number  FO4701 -74-C-0075. 
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II  Source-Absorber  Line  Correlation  Effects  - Line  by  Line  Calculations 

The  problem  being  considered  is  summarized  in  figure  1.  Of  the 
several  emission  and  attenuation  mechanisms  implied  by  figure  l^only 
the  molecular  absorption  by  gases  in  the  atmosphere  of  radiation 
emitted  by  hot  gas  in  the  target  is  considered.  Figure  2 summarizes  the 
basically  simple  equations  of  radiative  transfer  applicable  for  a single 
frequency  of  light  when  the  assumption  of  local  themodynamic  equilibrium 
is  valid.  Quasi-monochromatic  target  radiance  spectra  and  transmittance 
spectra  of  slant  paths  through  the  atmosphere  have  been  computed  for  the 
2.7  ||tn  spectral  region  based  on  the  relations  shown  in  figure  2 and  the 
Air  Force  Cambridge  Research  Laboratories  atlas  of  spectroscopic  line 
parameters,  The  specific  conditions  for  which  spectra  discussed  here 
were  computed  are  shown  in  figure  3.  A detailed  discussion  of  the 
techniques  may  be  found  in  reference  1.  The  quasi-monochromatic 
emission  arid  transmission  spectra  can  then  be  combined  to  provide  tho 
apparent  radiance  at  the  end  of  an  atmospheric  path.  Figure  4 shows 
a short  segment  of  such  an  apparent  radiance  spectrum. 

Since  most  sensors  have  bandwidths  much  wider  than  the  width  of 
individual  lines,  it  is  convenient  to  consider  the  atmospheric  transmittance 
averaged  over  some  interval,  typically  a few  wavenumbers  wide  and  in- 
cluding a number  of  ?*nes.  However,  as  soon  as  the  bandwidth  exceeds 
that  of  an  individual  line,  the  numerical  value  of  the  "average"  trans- 
mittance depends  on  the  assumptions  made  about  the  spectral  character- 
istics of  the  source.  We  illustrate  this  by  considering  several  different 
definitions  of  average  transmittance  and  demonstrating  the  quantitative 
difference  in  transmittance  computed  by  these  various  definitions  from 
the  correct  quasi-monochromatic  spectra  computed  for  the  conditions 
listed  in  figure  3. 

The  effective  average  transmittance  T is  the  applicable  value  if 

© 

the  objective  of  sensor  measurements  is  the  inference  of  the  source 
radiance. 
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earthfc  surface.  Because  of  the  difficulty  of  measuring  extremely 
low  humidities  with  sensors  which  have  passed  through  a region  of 
high  humidity,  many  of  the  older  models  of  the  stratosphere  are  un- 
realistically wet.  Much  of  the  meteorological  data  routinely  provided 
by  the  USSR  also  indicates  impossibly  wet  stratospheric  conditions. 
Personnel  of  the  Air  Force  Air  Weather  service  inform  me  that  even  in 
the  troposphere  when  the  ambient  temperature  is  above  -40  C some  of 
the  hygrometer  systems  employed  in  the  past  could  be  significantly 
affected  by  solar  illumination  and  thus  produce  erroneous  results. 

There  are  a number  of  reasons  why  meteorological  data  of  the  desired 
accuracy  may  not  be  available  for  many  events  of  Interest, 

The  conversion  of  meteorological  measurements  to  the  required  in- 
put for  transmittance  modelling  codes  must  be  done  carefully  to  avoid 
introducing  errors.  For  example,  if  relative  humidities  are  reported 
for  temperatures  below  0 C,  are  they  referred  to  saturation  vapor 
pressures  over  ice  or  over  water? 

Finally,  the  transmittance  codes  themselves  have  limitations.  One 
limitation  is  in  the  area  of  aerosol  attenuation  computation  discussed  by 
others  at  this  conference.  Another  is  the  limitations  Imposed  by  the 
approximations  required  to  obtain  computationally  efficient  band  models. 

The  spectrum  of  the  source  observed  through  the  atmosphere  can 
affect  the  apparent  transmittance  of  the  atmosphere,  as  alluded  to 
previously  by  Dr.  Wolfhard,  The  rest  of  this  paper  is  devoted  to  a 
quantitative  demonstration  of  the  size  of  this  particular  effect  by  means 
of  quasi.-monochrom&tic  spectral  computations  and  subsequently  s 
discussion  of  a band  model  transmittance  approach  which  includes  these 
effects. 
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(1) 


J c ^ (v)  dv  j'1  Lt  (s.)  T (v)  dv 
V1  V1 


-V  5 NJ 

J h ^ ^ / 

Vj 


Lt  (v)  dv 


Is  the  ratio  of  the  integrated  apparent  radiance  (v)  to  the  integrated 

target  radiance  L.(v).  Spectra  of  T for  two  optical  paths  are  indicated  by 

solid  lines  in  Figsw  5 and  6.  The  value  Av  a vg  ” Vj  ■ 20  cm-1  was  used 

in  computing  these  spectra.  Three  values  of  T.  obtained  for  bandwidths 
- 1 ® 

greater  than  20  cm  in  the  two  wings  and  the  center  of  the  absorption 

region  shown  in  Figs.  5 and  6 are  listed  in  Table  1.  The  entry  in  Table  1 

-1  __ 
titled  "20  cm  Intermediate  Average"  is  the  average  of  the  T curve  from 

Fig..  5 or  6 for  the  same  radiometer  bandwidth.  This  indicates  that  the 

average  of  low  resolution  spectra  does  not  produce  the  same  results  as 

the  average  obtained  from  high  resolution  spectra. 

The  average  transmittance T is  the  most  frequontly  computed  or 

approximated  trahsmittance. 


T « ± f Z T dv 


(2) 


The  average  transmittance  for  the  two  sample  paths  was 'computed  from  the 

» 1 

high  resolution  transmittance  spectra  by  averaging  over  B 20  cm  . Re- 
sults are  Indicated  by  dashed  curves  in  Figs.  5 and  6.  Average  tr&nsmittances 
for  wider  bands  are  shown  in  Table  1.  For  sources  of  uniform  spectral 

radiance,  T should  bb  the  same  as?  , but  for  line  sources  such  as  hot  gases, 

6 

they  clearly  are  not  the  same. 

The  factor  T is  the  quantity  approximated  by  most  atmospheric  trans- 
mittance band  model  procedures.  Estimates  ofT  by  two  such  programs  for 
the  two  sample  paths  are  included  for  comparison  in  Table  1.  In  the  table, 
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LOWTRAN  refers  to  the  program  developed  by  R.A.  McClatchey  at  the 

2 3 

Air  Force  Cambridge  Research  Laboratories,  ’ and  ATLES  refers  to 
a band  model  developed  at  Aerospace  to  account  for  source  effects.* 

Spectra  computed  with  the  use  of  LOWTRAN  are  indicated  by  the  broken 
curves  in  Figs.  5 and  6.  The  comparisons  given  in  Table  1 indicate  that  the 
band  models  approximate  T reasonably  well.  The  discrepancy  between  the 
T spectrum  and  the  LOWTRAN  estimate  for  the  75-deg  slant  path  is  not 
understood  at  present. 

An  attempt  to  include  source  effects  can  be  made  by  using  low  resolu- 
tion source  and  transmittance  functions  to  compute  low  resolution  effective 
transmittance 

T ■ ^1 T (Avi)  Lt  W (3) 

L^ 

The  results  of  carrying  this  out  for  the  wings  and  center  of  the  2.7 
absorption  are  given  in  Table  1.  The  accuracy  of  estimating  T#  in  this  way 
is  only  a little  better  than  the  average  transmittance  T obtained  without 
considering  the  source  spectrum  in  any  way. 

The  techniques  of  band  modelling  can  be  applied  directly  to  the  definition 
of  the  effective  average  transmittance,  Eq.  (1),  to  obtain  These  pro- 

cedures are  discuseed  later  in  this  paper  and  in  detail  in  Ref.  4.  Included 
in  Table  1 is  the  radiometer  transmittance  predicted  by  this  special  band 
model,  which  is  the  only  band  model  procedure  that  explicitly  takes  into 
account  line -correlation  relations  between  the  source  and  atmosphere  in 
estimating  T , It  is  also  the  only  procedure  that  does  not  overestimate 
the  transmittance.  In  fact,  if  any  conclusions  may  possibly  be  drawn 
about  the  band  model  from  this  limited  set  of  data,  it  would  appear  that 
the  model  overestimates  the  effects  of  line  correlation. 

Based  on  the  limited  test  cases  presented  in  this  discussion,  one  may 
conclude  that  all  average  transmittance  procedures,  which  do  not  account 
for  possible  line  correlations,  tend  to  overestimate  the  effective  trans- 
mittance of  the  atmosphere  for  a line  source.  In  some  spectral  regions 
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and  over  lomt  atmospheric  paths,  this  can  lead  to  errors  of  up  to  79 
percent  in  the  estimated  effective  transmittance. 

Ill  Source -Absorber  Line  Correlation  Effects  - Band  Model  ATLES 

Presented  now  are  spectra  produced  by  the  band  model  approach  for 

estimating  the  effective  average  transmittance.  This  band  model* 

5 

is  based  on  the  work  of  Lindquist  and  Simmons  . In  this 
model,  the  hot  plume  is  assumed  to  be  part  of  the  optical  path  to  the 
observer.  The  resulting  path  is  therefore  highly  inhomogeneous  due 
primarily  to  the  high  temperatures  in  the  plume.  The  usual  Curtis  - 
Godson  approximation  for  accounting  for  inhomo genoities  along  the 
optical  path  is  no  longer  valid.  The  quantity  which  must  be  estimated 
accurately  in  the  equation  of  radiative  transport  is  the  derivative  of  the 
transmittance  along  the  optical  path.  The  ATLES  band  model  program  in 
based  on  making  an  estimation  of  this  derivative  using  the  Lindquist- 
Slmmons  approximation. 

Figures  7 and  8 compare  the  average  transmittance,  T,  spectra  com- 
puted with  the  ATLES  band  model  with  that  computed  by  quasl-monochromatic 
means  on  the  basis  of  Individual  lines  for  the  conditions  specified  in  figure 
3.  It  is  this  transmittance  that  is  applicable  for  determining  the  attenuation 
of  radiation  from  continuum  sources.  The  agreement  between  the  band 
model  calculations  and  the  quaai-monochromatic  computations  is  satisfactory. 

Figures  9 and  10  compare  the  average  effective  transmittance,  T , as 
obtained  with  the  ATLES  band  model  with  that  obtained  from  quasi-mono- 
chromatic  calculations  for  the  conditions  specified  in  figure  3.  For  both 
atmospheric  paths,the  band  model  estimate  of  T is  lower  than  that  obtained 
by  the  line  by  line  calculations,  even  though  they  are  both  based  on  the  same 
spectral  line  compilation. 

Up  to  this  point, only  a source  at  20  km  altitude  hae  been  considered. 
Figure  11  shows  the  average  transmittance,  T,  and  the  average  effective 
tranemlttance,  T , for  a spectral  interval  in  the  wing  of  the  2.  7p.m 
atmospheric  absorption  band.  Results  obtained  with  the  ATLES  model  and 
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with  the  AFCRL  program  LOWTRAN2  for  T ere  shown.  Shown  as 
isolated  points  are  results  obtained  from  quasi  •monochromatic  com- 
putations. The  agreement  between  the  line  by  line  results  and  ATLES 
is  adequate  except  for  the  tendency  of  ATLES  to  underestimate  T^.  As 
expected,the  LOWTRAN2  results  agree  well  with  the  T curve. 

The  surprising  feature  in  figure  11  is  the  large  absorption  at  high 
altitudes  evidenced  in  the  results.  This  large  apparent  absorption 
results  from  two  sources.  The  first  is  the  inapproprlateness  of  the 
Lorenta  pressure  broadened  line  shape  used  in  the  calculations  of 
figure  11  for  the  high  altitude  region.  The  second  effect  is  the  Improved 
correlation  of  source  and  absorber  line  positions  and  strengths  as  the 
plume  temperature  decreases  at  increasing  altitude. 

The  effect  of  the  line  shape  function  has  been  investigated  and  results 
shown  in  figure  12.  The  effectjas  expected,  is  to  Increase  the  high  altitude 
transmittance.  However,  even  with  the  more  appropriate  Voigt  line  shape 
there  is  still  an  appreciable  absorption  due  to  line  correlation  at  an 
altitude  of  60  km  as  indicated  by  both  the  ATLES  and  line  by  line  computations. 
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IV  Summary 


Significant  reductions  In  affective  atmospheric  transmittance  occur 
when  there  is  a correlation  between  the  emitting  and  absorbing  spectral 
features  and  when  the  sensor  has  a bandwidth  greater  than  a single 
spectral  line.  The  conditions  leading  to  this  effect  are  common  in  the 
observation  of  aircraft  and  missile  plumes  through  the  atmosphere  by 
many  types  of  sensors.  These  effects  can  be  demonstrated  by  quaal- 
monochromatic  spectral  computation  and  modelled  efficiently  by  im- 
proved band  model  techniques  described  briefly  here. 

The  reduction  in  transmittance  due  to  line  correlation  occurs  even 
at  very  high  altitudes  where  the  atmosphere  is  often  considered  completely 
transparent. 

At  Aerospace, we  are  continuing  to  examine  this  problem  by  using 
Improved  plume  models,  by  improving  the  band  model  formulation, and 
by  obtaining  a consistent  set  of  band  model  parameters  for  both  high  and 
low  temperatures. 
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Figure  Captions 


Atmospheric  Attenuation  problem  summary 

Basic  equations  describing  the  transfer  of  radiation  for  a single 
frequency  under  conditions  of  local  thermodynamic  equilibrium. 

Conditions  employed  to  compute  the  atmospheric  transmittance  curves 
presented  in  this  report. 

Apparent  spectral  radiance  of  a typical  target  after  traversing  a path 
through  the  atmosphere  from  20  km  to  space  at  a 75-deg  zenith  angle, 

Transmittance  computed  in  several  ways  for  a horizontal  path  100 

km  long  at  20  km  altitude  (the  solid  curve'  - is  the  effective  average 

transmittance  for  a typical  missile  plume.  The  dashed  curve  ------ 

is  the  average  transmittance.  These  two  curves  were  obtained  from 

high  resolution  calculations  averaged  over  20  cm"^,  The  

curve  is  the  approximation  to  the  average  transmittance  provided  by 
the  AFCRL  LOWTRAN  computer  program). 

Transmittance  computed  in  several  ways  for  a slant  path  from  20  km 

to  space  at  a 75-deg  zenith  angle  (the  solid  curve is  the  effective 

average  transmittance  for  a typical  missile  target.  The  dashed  curve 
------ ia  the  average  transmittance.  These  two  curves  were  obtained 

from  high  resolution  calculations.  The  curve  is  the  approxima- 

tion to  the  average  transmittance  provided  by  the  AFCRL  LOWTRAN 
code). 

Comparison  of  the  average  transmittance,  T,  computed  by  means  of 
the  ATLES  band  model  with  the  results  of  line  by  line  computations  for 
a 100  km  horizontal  path  at  20  km  altitude. 

Comparison  of  the  average  transmittance,  T,  computed  by  means  of 
the  ATLES  band  model  with  the  results  of  line  by  line  computations  for 
a path  from  20  km  to  space  at  a 75  deg  zenith  angle. 

Comparison  of  the  effective  average  transmittance,  T , computed  by 
meane  of  ATLES  with  the  results  of  line  by  line  computations  for  a 100 
km  horizontal  path  at  20  km  altitude. 

Comparison  of  the  effective  average  transmittance,  T , computed  by 
means  of  ATLES  with  the  results  of  line  by  line  comptuAtions  for  a 
slant  path  from  20  krn  to  space  at  a 75  deg  zenith  angle. 

The  variation  in  atmospheric  transmittance  with  source  altitude  for  a 
path  from  the  source  altitude  to  space  at  a 73,7  deg  zenith  angle. 

The  effects  of  line  shape  function  on  the  computed  transmittance  at 
high  altitudes. 
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Table  1 . Radiometer  Transmittance 
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FIGURE  2 


Monochromatic  Atmospheric  Effects  Program 

INITIAL  TEST  CASE 
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• AVERAGE  TRANSMITTANCES  IN  WINGS  AND  BAND  CENTER 


v{cm~  | 
FIGURE  9. 


ATMOSPHERIC  TRANSMISSION  AND  EMISSION  PROGRAM 


I 


David  C.  Andlng  S 

Science  Applications,  Inc.  ; 

l.a  Jolla,  California 

•■t 


1 


91 


i 

;f 


'1 

v 


i'j  DQDriBnTKn  oi/ur  ar  AMV*_un<Tt 


ATMOSPHERIC  TRANSMISSION  AND  EMISSION  PROGRAM 


David  C.  Anding 
Science  Applications,  Inc. 

La  Jolla,  California 

Described  herein  is  an  atmospheric  transmission  and  emission  com- 
puter code  which  is  an  outgrowth  of  a study  performed  at  the  Infrared  In- 
formation and  Analysis  Centen  (IRIA)  at  The  University  of  Michigan’s 
Willow  Run  Laboratories,  in  1967.  Since  its  conception  in  1967  the  code 
has  undergone  many  revisions  and  improvements  in  selected  spectral 
regions.  The  result  is  a code  which  may  be  considered  state-of-the-art 
for  wavelengths  greater  than  4.0 nm,  but  remains  in  its  conceived  form 
for  shorter  wavelengths.  The  code  has  been  published  in  an  Aerospace 
report.  * 

Slide  1 delineates  the  basic  capabilities  of  the  code.  The  effects 
of  molecular  absorption,  aerosol  extinction,  and  thermal  emission  from 
both  molecules  and  aerosols  are  included.  The  applicable  spectral  re- 
gion extends  from  1 to  30  p,m.  The  molecules  whose  effects  are  repre- 
sented include  HgO,  COg,  Og,  CH^,  NgO,  HNOg  and  Ng.  The  basic 
assumptions  inherent  in  the  calculations!  procedure  are  local  thermo- 
dynamic equilibrium  and  single  scattering  for  aerosols.  The  effects  of 
refraction  are  not  included.  The  inputs  required  for  a given  calculation 
are  the  altitude  distributions  of  pressure,  temperature,  aerosol  density, 
HgO,  Og  and  HNOg  densities,  and  the  mixing  ratios  of  the  uniformly 
mixed  gases  - COg,  NgO  and  CH^.  The  outputs  of  the  code  are  trans- 
mission and  path  radiance  versus  wavelength. 

* J.  Hamilton,  J.  Rowe  and  D.  Anding,  Atmospheric _ Transmission .and 
Emission  Program,  Aerospace  Report  tfo.  tOr-0073  (3o66-0S)-3, 
June  1978. 
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The  code  is  divided  into  two  sub-codes,  CATM  and  ATMRAD,  which 
operate  in  tantum.  The  next  two  slides  present  the  specific  calculations 
that  are  performed  by  each  of  these  codes.  CATM  converts  rudimentary 
atmospheric  data  to  a spherically  shelled  model  atmosphere  which  is  re- 
quired as  input  to  ATMRAD,  the  code  which  performs  the  transmission  and 
radiance  calculation.  The  inputs  to  CATM  are  pressure,  temperature  and 
HgO,  Og  and  HNOg  densities,  each  expressed  as  a function  of  altitude  (at 
whatever  altitudes  they  are  available  to  the  user),  and  the  mixing  ratios  of 
C02,  NgO  and  Ch^.  The  outputs  are  the  input  parameter  values  re- 
evaluated at  the  prescribed  spherical  shells.  In  essence,  CATM  constructs 
a spherical  shelled  model  atmosphere  from  rudimentary  atmospheric  data. 
Five  standard  atmospheres  have  been  compiled  for  general  usage.  These 
are  denoted  Arctic  Summer,  Arctic  Winter,  Temperate  Summer,  Temperate 
Winter,  and  Tropic  Mean. 

ATMRAD  (Atmospheric  Radiance)  calculates  spectral  path  trans- 
mission and  radiance  from  1 to  30  Mm  for  any  geometric  path  within  a 
spherically  shelled  model  atmosphere.  ATMRAD  also  includes  the  radi- 
ance contribution  of  a Planckian  radiator  within  the  optical  line-of-sight. 

The  inputs  required  are  the  model  atmosphere  for  CATM,  the  geometric 
path,  the  wavelength  interval,  and  the  spectral  emlssivity  of  the  Planckian 
source.  The  outputs  are  spectral  path  transmission  and  total  radiance  at  the 
point  of  observation  along  the  direction  of  the  line-of-sight.  A schematic  of 
the  geometry  is  shown  in  slide  4.  A Planckian  source  is  denoted  Target, 
located  at  the  nth  shell,  and  the  point  of  observation  is  denoted  Detector, 
located  at  the  earth's  surface.  The  radiance  at  the  detector  (in  the  direc- 
tion of  the  line-of-sight)  is  calculated,  including  the  atmospheric  modified 
target  radiance  and  the  atmospheric  radiance  originating  along  the  line-of- 
sight  (radiance  originating  at  points  off  the  line-of-sight  impinging  upon 
the  detector  are  not  included).  The  number  of  shells  (n)  into  which  the 
model  atmosphere  is  divided  for  any  given  altitude  regime  is  an  option  of 
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the  user.  Generally,  as  the  path  becomes  more  horizontal  the  shell 
density  is  increased  in  proportion  to  the  secant  of  the  zenith  angle. 

As  shown  in  slide  5 the  path  transmission  is  the  product  of  the  in- 
dividual transmissions.  This  approximation  introduces  little  error  if 
either  the  transmission  is  slowly  varying  over  the  wavelength  interval  of 
the  applicable  transmission  model,  or  there  is  minimal  correlation  be- 
tween the  spectral  absorption  lines.  Also  shown  in  slide  5 is  the  equa- 
tion evaluated  for  the  computation  of  radiance. 

Slide  6 displays  the  transmission  models  used  by  the  code.  In 
total  there  are  seven;  Goody  model,  strong-line  Goody,  Elsasser, 
strong-line  Elsasser,  the  continuum  models  for  HgO  and  N2,  and  the 
aerosol  extinction  model.  All  parameters  have  the  common  interpretation. 
Slide  7 presents  the  transmission  models  that  are  used  for  the  respective 
molecules,  listing  the  approximate  resolution  for  each  wavelength  interval, 
the  procedure  that  was  used  to  obtain  the  transmission  model  coefficients, 
and  the  source  of  the  data  used  in  the  coefficient  evaluation.  The  asterisk 
denotes  those  transmission  models  and  coefficients  which  have  been  modi- 
fied since  the  code  was  conceived  in  1987. 

Slide  8 is  a comparison  of  the  one  set  of  aerosol  extinction  coeffici- 
ents used  by  ATMRAD  with  those  published  by  McClatchey  in  1974.*  This 
comparison  is  given  to  demonstrate  the  sensitivity  of  the  extinction  co- 
efficient to  particle  size  distribution  and  aerosol  complex  index,  both  of 
which  are  considerably  different  for  the  two  cases  shown.  Wavelength 
dependent  extinction  and  scattering  coefficients  are  specified  as  input  and 
as  such,  ATMRAD  can  perform  calculations  for  any  aerosol.  Five  addi- 
tional sets  of  aerosol  coefficients  are  available  for  use  ranging  from  100% 

* R.  A.  McClatchey  and  J.  E.  A.  Selby,  Atmospheric  Attenuation  of 

Laser  Radiation  from  0.76  to  31.25  um,  A^CR1j-TR-74-O0o3^ 

January  WTT. 


95 


maritime  haze  to  100%  continental  haze.  These  coefficients  are  from 
the  work  of  B.  Finn  of  AFCRL. 

Possibly  the  single  most  important  coefficients  in  the  6 to  14  pm 
spectral  region  are  those  for  the  water  vapor  continuum.  Shown  in  slide  9 
is  the  self  broadening  absorption  coefficient  used  by  ATMRAD  (solid  line) 
which  is  a least-squares  fit  to  the  measurement  data.  In  slide  10  trans- 
mission at  10.59  pm  (as  calculated  by  ATMRAD)  is  compared  with  mea- 
surements of  McCoy*  using  a foreign  to  self  broadening  ratio  of  0.005. 
Observe  that  the  absorption  is  overestimated,  which  when  extrapolated 
to  long  paths  at  high  humidity  can  be  considerable.  Because  of  the  con- 
sistency and  repeatability  of  the  McCoy  data,  and  hence  the  likelihood  of 
its  correctness,  it  was  felt  a modification  to  either  the  self-broadening 
coefficients  or  the  Kf/Kfl  ratio  was  in  order  to  bring  ATMRAD  results 
in  agreement  with  the  McCoy  data.  This  was  achieved  by  reducing  the 
ratio  of  Kf  to  Kg  to  a value  of  0.001.  Subsequent  discussions  with  Burch 
and  Long  indicated  that  a least -squares  fit  to  the  data  of  slide  9 probably 
yielded  a self -broadening  coefficient  that  is  too  high  because  of  possible 
systematic  errors  for  the  larger  data  points.  By  reducing  the  self-broad- 
ening coefficients  to  coincide  with  the  lowest  values  shown  in  slide  9, 
agreement  between  ATMRAD  computations  and  McCoy’s  measurements 
can  be  achieved  using  a value  for  Kj/Ks  of  0.005.  It  is  recommended 
that  the  available  data  base  be  carefully  reviewed  and  the  best  values  for 
both  self  and  foreign  broadening  be  selected  for  use  in  future  computations. 
As  a further  refinement  to  HgO  continuum  absorption  a temperature  de- 
pendence has  been  adopted  which  is  shown  in  slide  11. 

Since  the  publication  of  the  Aerospace  report  certain  improvements 
have  been  made  to  the  code  as  part  of  SAI'b  internally  funded  R & D 

* J.  McCoy,  D.  Rensch,  R Long,  Applied  Optics,  Vol.  8,  No.  1,  1909. 
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activities.  These  are  shown  in  slide  12.  The  re-evaluation  of  the  COg 
and  HgO  coefficients  was  done  primarily  to  increase  the  resolution  capa- 
bility. The  data  base  and  procedure  used  are  noted  on  the  slide. 

Slides  13  through  20  present  example  results  fran  the  ATMRAJD, 
including  comparisons  with  measurement  data,  line-by-line  calculations, 
and  LOWTRAN’n.  Slide  13  displays  the  spectral  region  from  1,0  to 
4.0  pm,  simply  to  typify  the  resolution  capability  for  this  spectral  region. 
Slide  14  displays  a comparison  between  ATMRAD  results  and  open  air 
field  measurements  performed  by  Convair.  The  measurement  path  param- 
eters are  noted  on  the  slide.  The  results  presented  in  slides  15,  16,  17 
and  18  were  generated  subsequent  to  the  recent  code  modifications  to  demon- 
strate consistency  between  band  model  results  and  the  data  base  used  to 
generate  the  band  model  coefficients.  Slide  15  is  a comparison  of  band 
model  results  for  GC^  with  the  measurement  data  of  Burch.  Slides  16, 

17,  and  18  are  comparisons  between  band  model  results  and  line-by-line 
results  (degraded  to  the  band  model  resolution)  for  HgQ  in  three  different 
spectral  regions.  For  each  case  the  parameters  are  noted  on  the  slides. 

Slide  19  is  a comparison  of  radiance  values  computed  by  ATMRAD 
with  measurement  data  from  Nimbus.  The  measurement  data  were  ac- 
quired by  the  Infrared  Interferometer  Spectrometer  (IRIS)  over  a cloud- 
free  ocean  area  for  which  the  surface  temperature  and  atmospheric  temper- 
ature and  humidity  were  known.  These  known  parameters  were  used  in 
ATMRAD  to  perform  the  radiance  calculations.  It  is  felt  that  the  results 
are  very  satisfying,  particularly  in  the  window  region  centered  at  1150  cm  *. 

Slide  20  is  a comparison  of  transmission  calculations  made  by  ATMRAD 
and  IiOWTRAN  for  a zenith  path  through  a common  atmosphere.  The  pre- 
cise reason  for  the  difference  between  the  two  calculations  is  not  known. 
However,  it  is  felt  that  the  discrepancy  probably  arises  from  slight  dif- 
ferences in  the  treatment  of  all  the  contributing  mechanisms,  i.  e. , con- 
tinuum absorption,  local  line  absorption,  and  aerosol  extinction. 
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Of  considerable  importance  to  many  applications  is  accurate  pre- 
diction of  8-14  pm  window  transmission,  particularly  when  the  trans- 
mission is  less  than  10%.  To  achieve  this  both  an  accurate  scattering 
model  is  required  and  an  accurate  representation  of  HgO  continuum  ab- 
sorption. The  scattering  models  are  discussed  elsewhere  in  these  pro- 
ceedings and  will  not  be  considered  further  here.  Based  upon  a brief  yet 
reasonably  thorough,  study  of  the  literature,  a few  statements  ' seeming 
our  understanding  of  H2C  continuum  have  been  made.  These  are  given  in 
slide  21. 

First,  there  still  remains  a lack  in  the  understanding  of  the  mechanism 
of  continuum  absorption.  Is  it  HjO  dimer,  the  far  wings  of  neighboring 
lines,  or  caused  by  contributions  of  both?  Second,  based  upon  the  lab- 
oratory data  of  Long  and  Burch,  self-induced  absorption  is  known  to  an 
accuracy  of  approximately  20%  at  room  temperature  for  wavelengths 
longer  than  8.0  pm.  Because  of  the  sparsity  of  data,  the  accuracy  is  con- 
siderably worse  at  shorter  wavelengths  and  at  temperatures  cooler  than 
approximately  296  Kelvins.  In  general,  for  terrestrial  temperatures,  the 
temperature  dependence  of  the  continuum  is  poorly  known.  Third,  the 
foreign-induced  absorption  coefficient  is  known  to  within  about  a factor  of 
2.  This  uncertainty  arises  because  the  foreign-induced  coefficient  has  not 
been  measured  directly,  but  must  be  inferred  from  measurements  in  which 
both  self-  and  foreign -Induced  absorption  are  present,  using  the  self-induced 
coefficient  as  a basis.  Small  errors  in  the  self-induced  coefficient  can 
cause  large  errors  in  the  foreign -induced  coefficient.  La  -itly,  and  of  con- 
siderable relevance,  measurements  of  continuum  absorption  have  not  been 
made  for  H20  partial  pressures  greater  than  approximately  75  percent  of 
the  saturation  vapor  pressure.  Therefore,  the  application  of  the  data  base 
to  predictions  of  atmospheric  transmission  in  tropical  atmospheres  where 
the  humidities  are  near  00%  results  in  an  extrapolation  of  the  data  and  may 
result  in  considerable  error,  particularly  for  long  paths. 
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In  summary,  ATM  RAD  may  be  considered  a state-of-the-art  band 
model  transmission  and  emission  code  for  wavelengths  greater  than  4.0  pm 
in  regard  to  its  treatment  of  molecular  absorption  and  emission.  For 
shorter  wavelengths  ATMRAD  is  vintage  1967.  For  aerosol  extinction 
ATMRAD  yields  results  that  are  consistent  with  the  accuracy  of  the  ex- 
tinction and  scattering  coefficients  used  as  input.  As  new  aerosol  param- 
eters are  disseminated  from  AFCRL,  these  can  be  used  as  input  to 
ATMRAD. 
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ATMOSPHERIC  AEROSOLS:  MODELS  OF  THEIR  OPTICAL  PROPERTIES 
Eric  P.  Shettle,  Robert  W.  Fenn  arid  Frederic  E.  Volz 

AFCRL 

In  previous  AFCRL  reports  (Elterman  1964,  1968,  and  1970,  McClatchey 
a.o.  1 970) atmospheric  models  have  been  presented  which  can  be  used  to 
derive  the  transmission  of  visible  and  Infrared  radiation  along  a given 
path  through  the  atmosphere.  The  models  developed  by  Elterman  consider 
atmospheric  light  attenuation  due  to  scattering  by  aerosol  particles  and 
air  molecules  and  absorption  by  ozone  molecules.  IR  transmission  calcula- 
tions In  these  models,  therefore,  were  limited  to  a few  wavelength  values 
In  watervapor  windows  where  It  was  assumed  that  molecular  absorption  Is 
negligible.  The  models  developed  by  McClatchey  arid  others  Included  the 
absorption  by  a 1 1 major  molecular  constituents  In  the  atmosphere  In 
addition  to  aerosol  scattering  and  absorption. 

The  aerosol  component  In  these  models  was  based  on  experimental 
measurements  which  were  made  during  and  prior  to  the  mid  1960’s.  At 
this  time  there  was  sufficient  experimental  data  available  to  define 
an  average  stratospheric  and  upper  tropospheric  aerosol  profile  with 
some  different  haze  concentrations  In  the  lower  troposphere  (up  to  a 
few  km  altitude)  with  exponential  vertical  decrease  In  particle  concentra- 
tion. 

The  vertical  distribution  of  aerosol  attenuation  In  the  upper  tropo- 
sphere and  stratosphere  In  these  models  was  primarily  based  on  several 
years  of  searchlight  measurements  In  a fixed  location  (Elterman  1966, 

and  1 968;  also  Elterman  et  al.  1969).  Ivlev  (1967  and  1969)  has  made 
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a review  of  the  avallaai'i  "Aoerlmenta I datn  ud  rnrough  1967.  Based  on 
this  review  Ivlev  d re sirred  a mode!  of  tne  verttcal  distribution  of 
aerosol  particles  and  their  extinction  of  visible  light,  Ivlev's  model 
Is  similar  to  Elterman's  (1968,  1970),  which  also  forms  the  basts  for 
aerosol  models  used  by  McClatchey  at  at.  (1970),  up  to  a height  of  about 
30  km.  Above  this  height  Ivlev's  model  diverges  rapidly  becoming  several 
orders  of  magnitude  I a roe r than  the  • del  of  McClatchey  et  al.  (1970). 

During  the  past  decade  In  this  country  and  elsewhere  extensive 
additional  measurements  from  ground  as  well  as  airborne  and  space  plat- 
forms have  been  made  of  aeresel  concentrations,  their  size  distribution, 
and  optical  properties,  to  warrant  the  development  of  updated  aerosol 
models  which  also  describe  some  of  the  temporal  and  spatial  variations 
In  atmospheric  aerosol  distributions  and  properties. 

One  result  of  particular  significance  from  these  recent  measure- 
ments Is  that  the  stratospheric  aerosol  concentration  during  the  middle 
and  late  1960's  was  still  above  normal  background  levels  (see  Elterman 
et  al.,  1973;  Fox  et  al.,  1973;  Hofmann  et  al.,  1974;  and  Russell  et  al., 
1974),  due  to  a residual  of  the  aerosols  Injectod  Into  the  stratosphere 
by  the  eruption  of  Mt.  Agung  during  the  spring  of  1963.  It  was  measure- 
ments made  during  this  time  period  of  elevated  aerosol  concentration  which 
served  as  the  major  Input  to  Elterman’s  and  Ivlev’s  models. 

In  this  study  a number  of  different  aerosol  models  for  each  of  4 
different  altitude  regimes  has  been  developed.  The  vertical  distribution 
of  the  attenuation  coefficients  for  these  models  Is  shown  lit  Fig.  I. 

(I)  For  the  Boundary  Layer  (below  2 km)  10  models  have  been  defined 
which  describe  rural  end  urban  environments  as  well  as  the  maritime  sea 
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aerosol,  for  several  surface  visibilities  between  2 and  50  km. 

(II)  For  the  upper  troposphere  there  are  two  models  which  represent 
spring  and  summer  conditions  versus  fall  and  winter  conditions. 

(III)  In  the  stratosphere  (up  to  30  km)  models  are  presented  for 
background,  moderate  and  high  volcanic  conditions  for  each  of  the  two 
seasonal  models.  The  vertical  distribution  for  the  moderate  volcanic 
models  for  both  seasons  Is  essentially  that  of  Elterman's  model. 

(Iv)  For  the  upper  atmosphere  (above  30  km)  two  models  are  being  j 

I 

used:  One  of  these  corresponds  to  the  most  likely  background  conditions 

and  the  other  represents  the  high  aerosol  concentrations  which  may  be  / 
observed  at  times  In  shallow  layers  at  various  altitudes  and  which  can 
be  of  significance  for  long  horizontal  propagation  paths. 

The  vertical  aerosol  concentration  profiles  and  the  size  distribu- 
tion are  described  by  analytic  functions. 

The  data  wl II  be  presented  so  that  any  model  for  one  region  can  be 
used  with  any  of  the  models  In  a different  altitude  regime.  So  In  effect 
It  Is  possible  to  compose  100  different  combinations  of  aerosol  models, 
covering  the  altitudes  between  0 and  100  km. 

For  each  of  these  models  the  coefficients  of  extinction  and  scat- 
tering as  a function  of  altitude  will  be  presented  for  20  wavelength 
values  between  0.25  and  40/um,  Including  11  laser  wavelengths  (Fig.  2). 

For  each  of  the  aerosol  components,  composing  the  various  models,  the 
coefficients  for  ext t net  I on,  scattering,  absorption,  the  particle  albedo 
and  the  scattering  asymmetry  coefficient  will  be  presented  as  a function 
of  wavelength.  The  angular  scattering  function  will  be  shown  for  each 
model  for  e few  significant  wavelength  values. 
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These  aerosol  models  and  their  optical  parameters  are  being  published 
In  a format  similar  to  that  of  the  AFCRL  Report:  Optical  Properties  of 
the  Atmosphere,  by  McClatchey  at  at.  (1970)  to  allow  for  transmittance 
calculations  Including  molecular  absorption  and  aerosol  attenuation. 

The  optical  properties  of  these  Models  will  be  compiled  Into  a computer 
program  subroutine  suitable  for  use  alth  the  LOWTRAN  program  (Selby  & 
McClatchey,  1972)  for  more  detailed  calculations.  The  report  on  these 
new  models  will  also  present  a discussion  of  the  aerosol  models  and  of 
the  experimental  data  foundation  for  them,  and  It  will  give  guidelines 
for  the  selection  of  the  proper  model  for  s specific  environmental 
condition. 

These  atmospheric  optical  models  In  their  present  form  do  not  give 

any  Information  with  regard  to  the  probability  of  occurrence  of  an>r 

\ 

particular  condition.  Yet  euch  Information  la  frequently  needed  by  the 
users  of  such  models. 

For  leek  of  better  date,  one  te  presently  forced  to  resort  to 
statistical  date  on  eurface  visibility  and  humidity  distributions. 

However,  these  peremeters  ere  not  adequate  to  derive  or  predict  slant 
path  visibility,  spectral  contrast  reduction,  IR  transmission  or  other 
complex  quantities.  An  Improvement  of  th'r  situation  can  only  coma  from 
maasuremants  which  sra  directed  towards  obtaining  statistical  data  on 
soma  of  the  basic  atmospheric  optlcal/IR  properties,  and  which  can  be 
used  to  derlvo  correlatlone  between  these  specific  optfcel  quentltfes 
and  the  more  routinely  observed  atmospheric  peremeters. 
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A measurement  program  which  has  these  objectives  Is  presently 
being  Implemented  under  a NATO  Research  Study  Group  effort  In  western 
Europe.  Under  this  measurement  program  of  the  Optical  Atmospheric 
Quantities  In  Europe  (OPAQUE)  a network  of  at  least  6 stations  will  be 
set  up  to  measure  on  an  hourly  basis  24  hours  a day  a minimum  sot  of 
visible  and  IR  atmospheric  parameters  (Fig.  3).  Some  stations  will  also 
record  some  additional  properties  such  as  angular  llaht  scattering 
functions,  aerosol  content,  IR  sky  and  terrain  radiance.  These  stations 
will  be  Instrumented  Jointly  by  several  NATO  countries  (Fig.  4)  and 
located  in  various  regions  of  western  Europe.  The  measurements  are 
scheduled  to  start  In  fall  1975  and  continue  for  a 2-year  period. 
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ABSTRACT 


This  report  describee  computer  models  which  have  been  developed 
to  determine  the  effects  of  atmospheric  constituents  on  optical 
sensor  systems  and  used  primarily  to  model  the  effectB  of  aerosols , 
such  as  clouds,  fogs,  and  smokes.  Models  which  calculate  response 
to  only  single-order  scattering  are  described,  and  sample  results 
obtained  by  applying  these  models  are  presented.  A more  sophisti- 
cated model  which  uses  Monte  Carlo  techniques  is  also  described. 
This  latter  model  possesses  the  capability  of  modeling  a plane- 
parallel  atmosphere  of  up  to  100  layers  containing  both  aerosol 
and  molecular  constituents. 
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INTRODUCTION 


Several  computer  models  have  been  developed  at  NWC  and  have  been 
used  extensively  to  model  the  effectB  of  aerosols  (clouds,  fogs  and 
smokes)  on  various  active  optical  sensor  systems.  The  investigation 
has  three  major  goals: 

1.  To  develop  viable  methods  of  discriminating  between  aerosol 
and  target  reflections, 

2.  To  develop  computer  models  that  characterize  Bensor  system 
performance. 

3.  To  perform  design-optimization  studies  of  sensor  systems. 

All  of  thess  have  been  accomplished  to  some  extent;  however,  work  is 
continuing  on  the  development  of  better  techniques  for  discriminating 
between  aerosol  and  target  backscatter. 

The  purpose  of  this  report  is  to  explain  the  utility  and  applications 
of  these  modeling  techniques.  Illustrations  of  typical  results  are 
presented  and  suggestions  for  other  areas  of  application  are  put  forward. 

The  rationale  for  performing  theoretical  analysis  of  this  nature 
is  that  optical  system  designers  must  determine  the  effects  of  aerosols 
on  moat  optical  systems.  Experimental  data  collected  in  natural  environ- 
ments are  difficult  and  expensive  to  obtain.  On  the  other  hand,  theoreti- 
cal models  provide  information  rapidly  and  are  inexpensive  to  exercise. 

Use  of  modeling  capabilities  allows  the  optical  designer  first  to  perform 
a design  optimization  and  then  to  collect  experimental  data  as  needed 
at  a much  lower  cost. 

Figure  1 illustrates  the  aerosol  problem  for  a generic  optical 
sensor  system.  The  target  signal  1b  attenuated  by  the  aerosol,  and  back- 
scatter  from  the  aerosol  produces  clutter  and  false  alarms.  Figure  2 
illustrates  the  aerosol  problem  in  communications  systems.  Here,  the 
aerosol  attenuates  the  desired  signal  and  also  distorts  the  modulation 
content  of  the  signal.  Other  systems  investigations  where  aerosol 
modeling  techniques  might  be  applicable  include  active  and  semiactive 
seekers,  target  designators,  range  finderB,  landing  systems,  use  of 
smoke  Bcreens  as  countermeasures,  and  laser  detection  systems. 
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• TARGET  SIGNAL  ATTENUATED 

• AEROSOL  BACKSCATTER  PRODUCES 
CLUTTER  AND  FALSE  ALARMS 

OPTICAL  TRANSCEIVER 

FIGURE  1.  Aerosol  Problem  for  Optical  Radar. 


SINGLE  SCATTERING  MODEL 


The  first  type  of  model  discussed  1b  one  that  calculates  only 
first-order  scattering  from  homogeneous  aerosol  clouds.1  This  model  uses 
Mis  light  scattering  theory  to  calculate  the  intensity  of  light  scattered 
from  spherically  shaped  particles.  The  particle-size  distribution  and 
the  complex  index  of  refraction  of  the  scattering  medium  are  required 
for  use  of  this  model.  With  a particular  beam  geometry  defined,  this 
model  will  calculate  the  extinction  of  a target  signal  due  to  the  aerosol, 
and  the  intensity  of  the  light  scattered  from  the  aerosol  itself.  These 
calculations  can  be  performed  at  any  deelred  discrete  wavelength  and  for 
various  types  of  aerosols.  Figures  3 and  4 illustrate  the  capability  for 
calculating  the  extinction  coefficient  and  the  backscatter  phase  function 
at  several  wavelengths  in  the  visible  and  infrared  regions.  These 
calculations  were  performed  using  a Deirmendjian  Model  Cj_,  fair  weather 
cumulus  cloud  particle-size  distribution.^  The  extinction  coafficient 
is  fairly  constant  at  all  wavelengths,  but  the  backscatter  phase  function 
has  several  deep  minima  that  are  very  significant.  To  illustrate  the 
effect  of  these  minima,  the  beam  geometry  illustrated  in  Figure  5 was 
modeled  at  wavelengths  of  0.905,  3.0,  and  10.5  micrometers,  The  system 
is  immersed  in  the  fair  weather  cumulus  cloud,  and  a 100-nanoaecond 
square  pulse  is  emitted.  The  cloud  and  target  returns  are  Bhown  in 
Figures  6 through  8,  as  the  leading  edge  of  the  pulse  moves  out  in  range. 
This  corresponds  to  target  returns  from  all  ranges  from  zero  to  45  feet, 
and  cloud  return  from  the  depth  of  the  cloud  penetrated  by  the  leading 
edge  of  the  pulse.  A comparison  of  these  results  revealB  that  the  cloud 
backscatter  is  approximately  a factor  of  40  lower  at  3.0  micrometers 
than  at  0.905  micrometer  and  a factor  of  25  lower  at  10,5  micrometers 
than  at  0.905  micrometer.  However,  the  transmission  of  the  target 
signal  ia  approximately  a factor  of  2.5  larger  at  10.5  micrometers  than 
at  the  other  two  wavelengths,  which  results  in  an  additional  improvement 
in  the  target/cloud  backscatter  ratio  at  10.5  micrometers. 

The  effect  of  modeling  different  particle-size  distributions  1b 
shown  in  Figure  9.  Deirmendjian  size  distributions1  are  used  in  theae 
calculations  to  represent  a sea  fog  in  various  stages  of  development. 

The  mode  radius  is  the  radius  at  which  the  greatest  number  of  particles 
are  concentrated.  Five  size  distributions  with  five  different  mode 


’'"Naval  Weapons  Center.  Analytical  Models  for  the  Design  of  Lidav 
Systems  (U)  by  R.  E.  Bird.  China  Lake,  Calif.,  NWC,  October  1973.  (NWC 
TP  5576,  publication  CONFIDENTIAL. ) 

2 

Deirmendjian,  D.  Electromagnetic  Scattering  on  Spherical  Poly- 

dispersions,  New  York,  Elsevier,  1969.  290  pp, 
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radii  have  been  used  to  generate  these  data.  An  interesting  result  is 
that  the  extinction  coefficient  becomes  larger  at  a wavelength  of  10.6 
micrometers  than  that  for  a wavelength  of  0.905  micrometer  for  size 
distributions,  with  their  major  concentration  of  particles  with  radius 
greater  than  6 micrometers.  However,  the  backBcatter  coefficient 
continues  to  improve  at  10.6  micrometers  relative  to  0.905  micrometer 
as  the  particle  size  increases. 


MULTIPLE  SCATTERING  MODEL 


In  some  applications,  multiple-order  scattering  can  play  a significant 
role.  A computer  model  has  been  constructed  that  uses  Monte  Carlo  tech- 
niques to  calculate  all  orders  of  scattering  that  are  of  interest. 3 This 
technique  traces  photon  trajectories  through  random  scattering  events 
with  built-in  efficiency  techniques  which  greatly  reduce  the  computer 
time  required  to  produce  a good  statistical  sample.  An  Infinite  plane- 
parallel  atmosphere  with  up  to  100  layers  can  be  modeled.  The  aerosol 
and  molecular  content  can  be  varied  in  each  successive  layer,  and 
scattering  and  absorption  due  to  both  of  these  constituents  can  be 
calculated.  The  model  can  simulate  conical,  fan,  and  point-beam 
geometries  with  no  restrictions  on  the  beam  orientation.  The  polarization 
state  of  the  scattered  light  is  calculated  through  the  use  of  Stokes 
parameters.  An  illustration  of  these  capabilities  is  shown  in  Figure  10, 
Figure  11  illustrates  one  possible  profile  for  naturally  occurring 
aerosols  as  a function  of  altitude;  Figure  12  is  a possible  profile  for 
ozone  aa  a function  of  altitude.  A good  representation  of  these  profiles 
as  well  as  other  atmospheric  constituents  can  be  modeled  in  thiB  program. 

Experimental  data  were  taken  under  controlled  conditions  in  a fog 
simulation  facility  for  the  purpose  of  checking  the  accuracy  of  this 
model.  The  results  of  this  comparison  are  Bhown  in  Figures  13  and  14 
for  two  different  beam  geometries.  The  Monte  Carlo  results  agree  very 
well  with  the  experimental  data  in  most  situations. 

Comparisons  were  also  made  between  several  second-order  scattering 
theories  and  the  Monte  Carlo  results.  An  example  of  a comparison  with 
data  generated  at  the  University  of  Florida*^  is  shown  in  Figure  15. 


3Naval  Weapons  Center.  Calculations  of  Multiple- Scattering  Effects 
on  Active,  Optical  Sensors  in  Cloud  Environments  (U)  by  R.  E.  Bird,  China 
Lake,  Calif.,  NWC,  August  197A.  (NWC  TP  5667,  publication  UNCLASSIFIED.) 

^Anderson,  R.  C.  and  E.  V.  Browell,  "First-and-Second-Order  Back- 
scattering  from  Clouds  Illuminated  by  Finite  Beams,"  APPL  OPT,  Vol,  II 
(1972),  pp.  1345-51. 
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The  agreement  la  quite  satisfactory  for  this  particular  comparison,  but 
poorer  agreement  was  obtained  with  other  second-order  Bc&ttering  models. 

One  application  where  multiple  scattering  can  have  an  important 
effect  is  in  the  shape  of  the  return  for  a short-pulse  syBtem.  Figure 
16  illustrates  the  pulse  shape  calculated  using  the  single  scattering 
model  and  the  Monte  Carlo  model  for  exiting  a cloud.  The  returned  pulse 
as  a function  of  time  is  plotted  for  a 10-nanosecond  Bquare  pulse  for 
both  models,  and  for  a 15-nanosecond  Bine  pulse  for  the  multiple- 
scetterlng  case.  The  effect  of  multiple  scattering  on  both  the  leading 
and  trailing  edges  of  this  pulse  is  readily  apparent. 


CONCLUSION 

The  computer  models  described  here  have  proven  extremely  useful 
for  optical  design  optimization  studies  and  for  modeling  the  effects  of 
aaroaols  on  optical  systems  under  various  conditions.  These  modelB  are 
suitable,  for  application  to  optical  trackers,  seekers,  target  deaignatora, 
landing  systems,  communications  systems,  and  laser  detection  systems. 
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FIGURE  6.  Target  and  Cloud  Return  for  Pulsed 
0,905-Micrometer  Light  with  Option  Immerued. 


FIGURE  7.  Target  and  Cloud  Return  for  Pulsed 
3,0-Micromater  Light  with  Optics  Immersed. 
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FIGURE  8.  Target  and  Cloud  Return  for  PulBed 
10.5-Micrometer  Light  with  Optics  Immersed. 
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FIGURE  10.  Multiple  Scattering  Anriysis. 
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FIGURE  14.  Target  and  Fog  Return  Signal  Versus  Meteorological 
Range  Through  18  Feet  of  Fog  for  2-Degree  Source  and  4-Degree 
Recelvar. 
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FIGURE  15,  Comparieon  Batvaen  Second-Order  Scatter  and  Monte  Carlo. 
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While  very  impressive  progress  has  been  made  in  modeling  atmospheric  molecular 
absorption,  the  aerosol  portions  of  existing  atmospheric  models  apply  only  to 
long-term  average  conditions.  It  has  been  pointed  out  (reference  1)  that  optical 
system  designers  require  information  regarding  the  departure  of  local  meteorological 
optical  conditions  from  the  long-term  average  conditions  currently  represented  by 
our  atmospheric  models.  This  paper  compares  techniques  for  modeling  aerosol  trans- 
mission effects  and  compares  Filippov  and  Mirumyants  models  with  those  of  Elterman  (3) 
and  McClatchey,  at  al  (4). 

The  modeling  techniques  currently  in  use  can  be  grouped  according  to  the  fundamental 
experimental  data  on  which  the  model  is  based,  Both  Elterman  and  Filippov's  models 
rely  on  direct  attenuation  vs  wavelength  measurements,  In  contrast,  McClatchey,  at  al, 
and  the  techniques  used  by  Bernhardt  and  streets  (5),  and  Hodges  (6),  for  example,  rely 
on  experimental  measurements  of  particle  size  distributions  and  indices  of  refraction. 
In  the  latter  case,  the  attenuation  versus  wavelength  is  calculated  using  Mie  theory. 

The  development  of  useful  atmospheric  models  requires  a blend  of  both  experimental 
data  and  theory.  It  is  the  author's  opinion  that  the  current  aerosol  models  should 
be  extended  to  include  an  indication  of  expected  departures  from  long-term  average 
conditions.  Flllipov  has  attempted  to  improve  aerosol  models  using  direct  attenuation 
measurements, 

Filippov's  data  Indicates  that  a classification  of  "optical  weather"  at  a particular 
site  by  season  and  weather  type  is  possible,  One  can  obtain  an  indication  of  expected 
deviations  of  local  conditions  from  long-term  average  conditions  by  comparing 
Filippov's  data  with  the  McClatchey,  et  al,  model. 
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Filippov  uaed  an  empirically  derived  expression  for  the  attenuation  coefficient, 

a (X)  “ [Kq  + Kx  X-n] 

where  X is  the  wavelength  in  microns  and  V ie  the  meteorological  range.  The  values 
of  Kq,  K^,  and  n depend  on  the  weather  type  and  are  listed  in  Table  I,  The  trans- 
mission is  given  by 


T - e“a(X)R 


(2) 


where  R is  the  range.  Figure  1 shows  a comparison  of  equation  (2)  with  the  Elterman 
and  McCletchey  models  for  horizontal  propagation  where  V ■ 5 KM  and  R ■ 10  KM.  It 
is  apparent  that  the  transmission  for  a DF  laser  could  vary  from  80%  to  ^20% 
depending  on  the  weather  classification  even  though  the  meteorological  range  is  5 KM 
for  all  conditions.  The  McClatchey,  et  al  model  products  T « 40%  for  the  same  con- 
ditions, Figure  2 illustrates  an  extreme  example  which  illustrates  the  impact  of 
local  "optical  weather"  on  a practical  system.  The  case  considered  is  a DF  laser 
radar  propagating  over  a 10  KM  path  (R  ■ 20  KM  for  two-way  path)  with  a meteorological 
range  of  5 KM,  The  MeQlatchey  model  products  T • 15%,  whereas  Filippov's  data  covers 
the  range  of  T ■ 5%  to  60%.  These  comparisons,  while  strictly  applicable  only  to  a 
specific  geographical  location,  do  indicate  the  kind  of  deviation  from  long-term 
average  conditions  which  can  be  encountered. 


Two  specific  measures  should  be  pursued  to  improve  our  current  ir  aerosol  attenuation 
models.  First,  data  regarding  the  sensitivity  of  the  results  of  existing  aerosol 
models  to  changes  in  the  assumed  aerosol  characteristics  should  be  published. 
Secondly,  a direct  attenuation  measurement  program  should  be  initiated  to  begin  to 
provide  data  on  the  geographical  and  seasonal  variations  of  ir  attenuation.  This 
date  should  be  used  to  derive  empirical  models  for  local  conditions  and  to  check 
the  theoretical  predictions  which  are  based  on  assumed  particle  size  distributions. 
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TABLE  l 


r ■ 
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WEATHER  CONDITIONS 

K0 

K, 

n 

SPRING  AND  FALL  HAZE 

0.04 

mm 

1.02 

WINTER  HAZE 

0.0 

■H 

1.24 

STABLE  SUMMER  HAZE 

0.06 

1.88 

NEW  SUMMER  HAZE 

o 

• 

o 

Bfl 

1 .88 
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BACKGROUND  PAPER 
ON 

A DESIGN  STUDY  FOR  AN  AIRBORNE  INFRARED  TRANSMISSOMETER 

1.  This  paper  concerns  a design  study  f or  an  airborne  instrument  to 
measure  the  transmissivity  of  the  atmosphere  to  infrared  radiation 
between  the  wavelengths  of  2 and  13  micrometers.  The  design  study 
was  accomplished  by  Block  Engineering,  Incorporated,  of  Cambridge, 

MA  for  the  6585th  Test  Group  (AFSC)  at  Holloman  AFB,  NM.  This 
paper  addresses  the  potential  benefits  of  the  instrument,  its  specifica- 
tions, its  principles  of  operation,  its  possible  configurations,  and  its 
costs  of  procurement  and  operation.  This  paper  recommends  that  DOD 
designate  and  fund  an  agency  to  buy  and  fly  the  transmlBBometer  to  im- 
prove and  verify  the  accuracy  of  widely  used  transmissivity  models. 

2.  Benefits  of  an  Airborne  Infrared  Transmissometer: 

a.  The  immediate  benefit  of  a transmissometer  would  be  measure- 
ments of  atmospheric  transmissivity  during  tests  of  infrared  systems. 
Systems  sensitive  to  atmospheric  transmission  conditions  are  Forward 
Looking  Infrared  Imaging  Systems  (FLIRs),  imaging  infrared  guidance 
units  for  missiles,  air  to  air  missile  guidance  units,  infrared  trackers, 
and  remote  atmospheric  sensors.  For  example,  atmospheric  trans- 
missivity muBt  be  known  in  order  to  properly  analyze  resultB  of  com- 
parative tests  of  infrared  systems  that  are  not  accomplished  simultane- 
ously. Also,  knowledge  of  transmissivity  during  tests  of  IR  systems 
allows  analysts  to  define  their  performance  capabilities  under  various 
atmospheric  conditions.  Such  knowledge  is  important  for  determining 
proper  Inventories,  deployments,  and  employments  of  infrared  guided 
munitions.  * 

b.  The  most  important  benefit  of  a tranamisBometer  would  be  its 
capability  to  provide  data  needed  to  verify  models  of  atmospheric 
transmissivity.  Models  of  infrared  transmissivity  are  currently  used 
by  developers,  testers,  and  weapons  performance  analysts  to  assess 
potential  and  real  capabilities  of  infrared  systems.  The  accuracy  of 
these  models  &b  functions  of  various  atmospheric  constituents  has  not 
been  established.  Their  parameterization  of  the  effects  of  certain 
atmospheric  constituents,  such  as  aerosoLs,  needs  to  be  improved. 
Measurements  of  atmospheric  transmissivity  in  a variety  of  atmos- 
pheric conditions  could  provide  the  data  needed  to  verify  and  improve 

the  models.  Once  verified,  the  models  could  be  used  with  full  confidence. 
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A special  program  to  build  add  fly  the  airborne  infrared  transmit- 
icmaeter  to  gather  the  required  data  la  considered  appropriate. 
Simultaneous  measurements  of  atmospheric  parameters  with  trans- 
missivity would  be  required. 

tiu  A third  benefit  of  a transmit  some  ter  would  be  its  capability 
to  provide  basic  scientific  information  pertaining  to  spectral  locations 
of  absorption/emission  lines,  pressure  and  doppler  broadening  of 
theta  lines*  and  information  on  the  kind  and  concentrations  of  atmos- 
pheric constituents.  It  could  provide  the  first  spectral  Information 
available  on  Infrared  transmissivity  in  various  aerosols. 

d.  A particular  benefit  of  this  tranamlssometer  la  that  its  design 
study  Is  finished  and  paid  for. 

3.  Specifications  of  the  Airborne  Infrared  Transmissometer  are 
Illustrated  in  Figure  I.  They  are  as  follows: 

a.  Spectral  interval:  2>13/jm 

b.  Wavenumber  resolution:  ^0,5j:m“l 

c.  ’ Signal  to  noise  ratio:  > 10  to  1 (design  goal  =*  100  to  1) 

d.  Range: 

. (1)  Air  to  air:  500  to  25.  000  feet  when  extinction 
coefficient  <0.1  km*  * 

. (2)  Air  to  ground:  500  to  50.  000  feet  when  extinction 

coefficient  <0,  l km-1 

e.  Scanning  time:  1 to  60  seconds  (integration  time  for  an 

Interferometer) 

f.  Calibration: 

(1)  Internal  reference:  ±5% 

4.  Basic  Operating  Principles  of  Recommended  Design:  The  recom- 
mended design  by  Block  Engineering,  Incorporated  is  called  the 
Seurce-Interferometer-Retroreflector  (SIR)  configuration.  Four 
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additional  configurations  were  addressed  lathe  design  study  as  possi- 
bilities for  performing  transmissivity  measurements.  The  SIR  design 
is  illustrated  by  Figure  11.  A 2500  degree  Kelvin,  incandescent  source 
is  moduled  by  a moving  mirror  in  the  Block  type  Michelson  interfero- 
meter to  produce  an  interferogram  of  the  source  radiation,  which 
passes  through  transmitter  optics  into  the  atmosphere  for  transmission 
to  a retroreflector  aboard  a target  vehicle  or  on  the  ground.  The 
reiroreflector  returns  the  radiation  through  the  atmosphere  to  re- 
ceiver optics  ahead  of  & detector  which  senses  the  interferogram. 
Appropriate  electronics  are  employed  behind  the  detector  to  distin- 
guish the  interferogram  (AC  signal)  from  background  radiation  (DC 
signal)  within  the  field  of  view  of  the  detector  optics,  A Fourier 
transform  of  the  interferogram  and  proper  geometric  scaling  yields 
a high  resolution  transmissivity  spectrum. 

5.  Airborne  Infrared  Transmissometer  Test  Bed  Aircraft: 

a.  The  design  study  assumes  the  transmissometer  will  be  carried 

in  a space  stabilized  gimbal  platform  in  the  nose  of  a C-130  aircraft. 
This  gimbal  moves  in  two  dimensions  for  pointing  and  tracking  pur- 
poses, and  is  aimed  by  a visual  or~infrared  tracker  carried  with  and 
boresighted  to  the  transmissometer  optics.  This  arrangement  is 
illustrated  by  Figure  III.  • 

b.  An  alternate  approach,  illustrated  in  Figure  IV,  requiring 
additional  design  work,  would  be  to  use  pods  to  carry  the  transmis- 
someter as  well  as  the  retroreflector.  This  would  allow  for  more 
flexibility  in  the  choice  of  the  test  bed  aircraft.  Some  difficulty  may 
be  encountered  in  the  design  of  appropriate  protective  windows  ahead 
of  the  transmissometer  and  the  retroreflector.  A stabilized  gimbal 
may  or  may  not  be  required  in  the  pod  holding  the  transmissometer, 
depending  upon  the  abilities  of  the  pilot  to  properly  point  the  aircraft 
at  the  retroreflector  in  an  airborne  pod  or  at  a retroreflector  array 
on  the  ground. 

6,  Costs:  Figure  V indicates  likely  costs  to  be  incurred  by  the 
transmissometer.  The  basic  design  study  cost  $70K.  Design  of 
the  alternate  pod  configuration  is  estimated  to  cost  approximately 
$30K,  Costs  of  building  the  transmissometer  and  retroreflector  are 
estimated  to  be  $170K.  Various  configurations  of  data  processing 
and  graphics  display  equipment  are  estimated  to  cost  anywhere  be- 
tween $8QK  and  $300K,  the  Utter  for  a ruggedized,  airborne 

r sal -time  system. 
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C-I30  TEST  BED 


aircraft  pods 


FIGURE  IV 
6585TG 


7.  Other  Configurations:  Several  other  instrument  configurations 
were  investigated  by  Block  as  possibilities  to  accomplish  the  trans- 
missivity measurements.  These  are  discussed  and  illustrated  in 
the  design  study.  Their  primary  advantage  is  greater  range  capa- 
bility. Their  disadvantages  involve  greater  cost,  lower  accuracy, 
and  higher  pointing  and  tracking  requirements. 

8.  Summary:  The  design  of  an  instrument  capable  of  measuring 
the  transmissivity  of  the  atmosphere  in  an  airborne  environment 
has  been  completed.  Its  design  costs  have  been  paid.  This  instru- 
ment would  be  capable  of  providing  basic  spectral  transmissivity 
data  important  for  establishing  transmissivity  information  during 
tests  of  infrared  systems  or  to  verify  the  accuracy  of  various  com- 
puter models  of  transmissivity.  We  believe  the  most  important 

and  most  economical  use  of  the  transmissometer  would  be  to  improve 
the  accuracy  of  existing  transmissivity  models.  These  models  could 
then  be  used  with  full  confidence  by  systems  developers,  testers,  and 
performance  analysts.  Transmissivity  of  the  atmosphere  must  be 
considered  when  deciding  (1)  which  weapons  to  develop,  (2)  which 
weapons  perform  the  best,  and  (3)  what  proper  force  levels  to  pro- 
cure and  maintain  in  various  operational  environments.  In  the 
future,  streamlined  transmissivity  models  might  also  be  used  by 
weather  support  personnel  to  indicate  proper  employment  tactics 
of  available  weapons  to  operations  personnel.  Accurate  transmis- 
sivity models  are  essential  in  order  to  relate  standard  meteorolog- 
ical variables  observed  world-wide  to  the  performance  capability 
of  Infrared  guided  munitions.  For  these  reasons,  this  paper  recom- 
mends that  DOO  designate  and  fund  an  agency  like  the  Air  Force 
Cambridge  Research  Laboratories  to  build  and  fly  the  Airborne 
Transmissometer,  in  order  to  gather  data  required  to  improve  and 
verify  their  widely  used  transmissivity  models.  This  approach 
would  seem  to  be  more  economical  than  for  one  or  more  test 
agencies  to  buy  and  fly  a transmissometer  during  tests  of  infrared 
systems. 
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COMPARISONS  OF  TRANSMITTANCE  CALCULATIONS  BY  WO  METHODS 


A.  J.  LaRocca 

Environmental  Research  Institute  of  Michigan 


I was  invited  to  attend  this  meeting  after  hearing  from  Dave  Anding, 
who  suggested  that  I present  some  of  my  results  showing  comparison  of 
calculated  values  of  transmittances  using  the  LOWTRAN  2 model  and  the 
model  developed  by  him  and  others  while  he  was  at  the  Environmental 
Research  Institute  of  Michigan  (ERIM).  The  comparisons  are  part  of 
the  product  of  a State-of-the-Art  Report,  sponsored  under  a contract 
with  the  Infrared  Information  and  Analysis  (IRIA)  Center,  entitled 
"Atmospheric  Transmittance  and  Radiance:  Methods  of  Calculation,"  by 
LaRocca  and  Turner.  The  report  is  now  in  publication  and  will  be  dis- 
tributed to  IR1A-1RIS  subscribers  in  about  a month  or  so. 

Because  the  report  is  pertinent  to  the  subject  matter  of  this 
meeting,  a few  brief  statements  about  its  content  are  in  order.  This 
is  probably  best  done  by  quoting  the  abstract  which  is  reproduced  as 
follows: 


"The  effort  represented  by  this  report  was  a result  of 
the  need  to  bring  a description  of  the  state-of-the-art  of 
mathods  of  calculating  atmospheric  transmittance  and  radiance 
up-to-date.  The  report  is  broadly  dividsd  into  ths  categories 
of  scsttsrlng  and  absorption,  with  the  greater  striae  laid  on 
absorption.  The  essential  material  la  presented  in  Sections 
3,  6,  7,  end  8,  in  which  specific  methods  are  described. 

Section  3 is  devoted  to  scattering  calculation  techniques, 
while  Sections  6,  7,  and  8 cover  methods  of  calculating 
transmittance.  The  first  of  these  is  the  Bo-called  line-by- 
line direct  integration  method,  which  requires  a detailed 
compilation  of  the  characteristics  of  individual  molecular 
lines.  Some  familiar  numerical  integration  techniques  are 
used  to  effect  quadrature  in  the  most  convenient  and  economical 
way. 

"The  second  of  the  absorption  methods  of  calculation  pre- 
sented is  the  band-model  technique.  In  this  method,  the  line 
spectrum  is  approximated  by  some  mathematically  manlpulatable 
distribution  function  with  undetermined  band-model  parameters. 
By  comparison  of  calculated  results  with  laboratory  experi- 
mental data  the  parameters  are  defined,  and  the  band-model 
ia  used  for  calculating  transmittance  under  any  required 
meteorological  conditions. 
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"The  third  general  eet  of  techniques  is  given  the  heading 
’Multi-Parameter  '-ulytical  Procedures. ' These  techniques  are 
derived  from  the  band-model  concept,  incorporating  a larger 
number  of  parameters,  with  presumably  greater  accuracy  in  the 
resultant  calculations. 

"The  rest  of  the  report  iB  either  tutorial  or  supportive, 
presenting  details  of  information  which  is  required  as  input 
to  the  calculation  procedures.  The  major  input  is  the 
meteorology  required  to  describe  absorber  concentrations, 
pressures,  temperatures,  and  other  necessary  physical  entities. 

"An  assessment  is  made  of  the  various  techniques  of  cal- 
culation in  terms  of  accuracy,  computer  time  needed  to  perform 
calculations,  adaptability  to  specific  problems,  and  practi- 
cality." 

The  specific  purpose  of  this  presentation  is  to  show  a comparison 
between  the  results  alluded  to  above.  The  method  developed  by  Anding 
is  entitled  the  Aggregate  Method  in  the  State-of-the-Art  report 
because  of  the  nature  of  the  approach  to  the  calculation.  In  the 
Foreword  to  the  report  it  iB  stated  that  the  cut-off  period  for  material 
contained  in  it  is  around  mid-1974.  This  makes  the  results  about  to 
be  shown  a bit  out-of-date  because  Anding  has  since  changed  his  make- 
up of  the  Aggregate  method  to  some  extent,  as  he  related  to  me  shortly 
before  the  meeting.  And  we  have  heard  in  the  course  of  this  meeting 
that  LOWT RAN  2 is  about  to  be  changed  to  LOWTRAN  3. 

The  calculations  were  made  with  all  meteorological  Inputs  the 
same  so  that  the  comparisons  would  point  out  differences  only  in  the 
techniques  for  performing  the  calculations.  In  order  to  show  in  this 
presentation  the  range  of  whatever  differences  occur  I have  Included 
two  standardised  atmospheric  models,  one  relatively  dry  and  one 
relatively  wet.  The  dry  model  is  represented  by  what  is  called  an 
Arctic-Winter  atmosphere  and  the  wet  by  a Tropic  atmosphere. 

Figure  117  shown  a comparison  between  the  transmittance  calculated 
by  the  Aggregate  and  LOWTRAN  2 method.  The  figure  is  divided  into 
an  (a)  and  (b)  parts  representing  respectively  the  long  and  short 
wavelength  parts  of  the  infrared  spectrum.  Figure  117  is  for  an 
Arctic  Winter  atmosphere.  The  comparison  is  reasonable  good  showing 
some  compatibility  in  the  two  methods  for  a fairly  dry  atmosphere. 
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Figure  113,  (a)  and  (b),  shows  the  Bane  type  of  comparison  for 
a Tropic  atmosphere.  The  divergence  in  transmittance  values  here 
i.-i  somewhat  larger,  especially  in  certain  spectral  regions.  Some  of 
this  difference,  according  to  Anding,  is  attributed  to  differences 
in  the  method  of  handling  the  H^O  continuum,  as  well  as  differences 
in  the  coefficients  considered  in  the  two  methods.  It  is  most 
difficult  to  resolve  these  differences,  because  comparisons  with  an 
independent  method,  say  the  results  of  field  measurements,  is 
inadequate.  Reliable  field  measurements,  in  which  the  atmospheric 
conditions  can  be  reproduced  by  calculation,  are  hard  to  come  by. 
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WORKSHOP  ON  ATMOSPHERIC  TRANSMISSION  MODELING 

January  28,  1975 
Institute  for  Defease  Analysis 

Afternoon  Session:  Workshop  on  Physics  end  Engineering 

of  Modeling  the  Atmosphere 

Moderator:  J.  J.  Gallagher 

Engineering  Experiment  Station 
Georgia  Institute  of  Technology 

The  afternoon  session  was  an  open  session  on  the  physics  and  engineer- 
ing of  modeling  the  atmosphere.  The  session  vas  taped,  however,  one  tape 
did  not  record  satisfactorily,  resulting  in  considerable  difficulty  in 
transcribing  the  session.  The  difficulties  of  Rosemary  Wood  were  deeply 
appreciated  in  this  effort. 

The  session  was  started  with  a call  for  questions  and  comments  on  the 
papers  presented  in  the  morning.  The  first  paper  addressed  was  that  of 
Dr.  L.  Bibarman,  entitled  "The  User's  View  of  Atmospheric  Models:  1. 

Thermal  Imaging  Systems."  A question  was  asked  how  Biberman's  IR  measure- 
ments were  made. 

Bibarman:  In  1970,  a decision  was  made  to  gather  meteorological  data  for 

one  year.  I have  the  data  on  visibility,  relative  humidity,  dew  point, 
temperature,  all  the  usual  parameters.  I used  the  LOWTRAN  III  Modal. 

These  ere  calculated  reeulte  end  the  validity  hinges  quits  clearly  on 
validity  of  the  scattering  model,  and  1 pointed  out  that  when  things  get 
bad,  they  get  bad  because  of  scattering.  There's  not  much  water  in  January. 
When  stuff  goes  out  for  six  hours,  it's  snowing  or  there's  bad  ground  fog 
or  similar  atmospheric  conditions.  So  when  the  equipment  is  about  to  go 
out  of  operation  or  it's  only  operating  for  half  its  normal  range  or  less, 
ite  because  of  scattering,  and  therefore,  when  one  wants  to  teks  e look  at 
how  long  the  equipment  will  be  of  no  value  to  serve,  one  really  has  to 
consider  this  on  the  basis  of  the  scattering  model  that  dominates,  and 
therefore  it’s  necessary  that  the  scattering  model  be  reasonably  appropriate 
or  reasonably  accurate.  When  I go  back  to  the  Air  Force  end  sey  that  you 
better  not  build  that,  because  the  Air  Force  model  is  no  damnsd  good,  that's 
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about  where  we  stand.  Now,  it  turns  out  that  we  had  six-seven  years  of 
scatter  modeling  work,  and  we  still  can't  use  it.  I'm  going  to  put  all 
the  pressure  that  I can  into  getting  a scattering  model  that  will  plug  into 
a computable  program. 

Walsh:  Can  we  ask  you  to  give  a few  thoughts  on  what  the  content  of  a 
scatter  model  should  be?  One  difficulty  with  the  scatter  model  is  that  it's 
continuously  variable  from  100%  transmission  to  zero.  People  talk  about 
London  fogs  in  which  you  cannot  saa  your  hand  in  front  of  your  face. 

Blberman:  OK,  John,  that'e  critical. 

Walsh:  You  need  to  reach  some  understanding  on  what  kind  of  model  you  want 
to  hava.  It  haa  to  go  all  tha  way  from  zero  to  ona.  Do  you  want  a highly 
refined  model  that  correlatee  with  visibility?  Can  you  in  fact  corralata 
with  visibility  because  you  know  that  at  10  micromatera  it 'a  somewhat 
different  than  what  happens  in  the  visible.  What  should  be  try  to  do  to 
get  to  the  model?  What  should  tha  model  contain? 

Woodman:  Can  I ask  one  question  first?  When  you  eey  scatter  model,  are 

you  talking  about  probability  of  occurrence,  or  ere  you  talking  about 
scattering  by  particles  in  air? 

Bibarmaii:  Scattering  by  particles  in  sir.  I know  when  the  scattering  gets 

bad  because  any  weather  station  tells  me  that  the  visibility  is  bed. 

Woodman:  Tha  LOWTRAN  III  does  not  have  a hydrometeor  subroutine  in  it. 

Does  it? 

Selby:  No. 

Woodman:  So  when  you  say  that  it's  snowing  or  foggy... 

McClatchey:  That's  another  set  of  conditions  thst  one  might  address.  It 

could  be  handled  by  LOWTRAN  III.  I think  one  might... 

Fenn:  You  didn't  have  any  infinite  fog  situation  in  those  cases.  Did  you? 

Blberman:  All  I can  tall  you  is  that  whan  you  gat  meteorological  data 

recorded,  they  tell  you  whet  the  relative  humidity  is,  what  tha  temparature 
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is,  what  the  ground  temperature  is,  and  what  the  visibility  is.  Very  seldom 
there  is  a little  thing  that  says  it’s  raining. 

Selby;  The  visibilities  that  you  had  were  from  zero  or  2-10  km? 

Blberaan;  Correct.  Visibility  in  a couple  cases  went  to  zero.  As  far  as 
I’m  concerned,  that's  either  fog  or  enow,  and  I’m  not  worried  about  that. 

I'm  worried  when  the  visibility  is  2 km,  and  I would  like  to  see  a model 
that  says  this  model  is  typical  of  a high  sea  state  at  low  altitude  above 
the  ocean,  and  that  ain’t  distilled  water.  Someone  said  it's  boulibals, 
and  it  Is.  If  you  take  a look  at  all  models  of  what  aerosol  content  is, 
of  sea  and  why  bubbles  break  and  they  throw  stuff  up,  and  water  evaporates 
and  you  leave  salt  nuclei  all  over  the  place,  and  they  ere  not  small  in 
slza  and  they're  aerious,  then  you  need  a modal  like  that.  We  have  a model 
somewhere  over  White  Sends,  New  Mexico,  but  I don't  believe  thla  for  Hanover, 
Germany  or  coastal  regions  around  the  North  Sea,  so  I think  that  we  have  to 
have  something  like  coastal,  wt  havs  to  have  some thing  like  continental, 
they  have  to  say  that  we  arc  using  a scatter  that  is  typified  when  the 
visibility  runs  between  2 and  4',  and  it  is  in  suburban  or  rural  continental 
mld-lattitude,  and  we  hnvft  urban  under  some  conditions.  I think  the  scatter 
for  2-4  kilometers  is  different  from  one  that  has  10  kilometer  visibility, 
and  so  1 don't  think  that  we  have  one  continuous  model.  We  may  have  a 
series  of  step  functions,  but  your  program  must  be  able  to  correlate. 

Walsh;  It's  still  a problem  no  matter  how  much  you  know  that  people  will 
continue  to  ignore  it.  1 had  e call  last  weak  from  a Navy  contractor  way 
down  tha  lint  on  building  some  system  for  a CO^  laser,  and  he  is  suddenly 
concerned  about  the  high  attenuation  that  he's  going  to  encounter  in  rain 
end  fog.  I told  him  that  years  ago  they  should  never  have  started  such  a 
development. 

Blbermsa ; There  will  always  be  stupid  people.  Let’s  not  address  thoss. 

Let's  try  to  address  more  responsible  people  in  the  world  who  can  see  that 
thinge  happen  properly.  There's  no  point  in  taking  a look  at  the  dumb  heads. 
Thsre’s  enough  dumb  heads  around,  and  no  matter  what  we  do,  ve  will  not 
impress  the  dumb  heads.  What  I think  ve  have  to  do  is  supply  data  for  ths 


193 


thinking  individuals  who  know  how  to  use  real  data  and  understand  constraints, 
and  all  the  other  things  that  you  have  to  apply  to  the  data. 

Fetm:  The  problem  that  arises  though,  I think,  1b  to  evolve  too  many  models, 
a large  number  of  then,  and  1 fully  agree  that  I do  not  think  that  it's 
possible  to  ascribe  a large  variety  of  conditions  to  one  model  by  just 
varying  some  parameters.  There's  just  different  physical  conditions  that 
require  different  models  but  if  one  can  dascrlba  these  different  models, 
then  one  problem  that  we  run  into  ie  tha  difficulty  in  educating  the  ueera 
to  use  the  proper  model  for  this  situation. 

Blberman;  Listen,  there's  a lot  of  guys  that,  don't  know  what  to  do  with  a 
Bessel's  function,  but  that  doesn't  mean  we  shouldn't  have  Besaal's  functions, 

Fetm:  Mo,  but  1 think  that  It's  one  aspect  thst  ve  have  to  keep  in  mind  end 

take  care  of  at  the  same  time  that  one  brings  out  the  model. 

McCLatchey:  Well,  the  modal  will  probably  have  to  contain  that  kind  of 

information.  The  model  has  to  be  boiled  down  ultimately  to  a number,  a 
limited  number  of  input  meteorological  parameters  that  are  ordinarily 
measured  at  which  point  the  model  will  determine  that  if  you  have  a certain 
relative  humidity  and  a certain  visibility  and  other  pieces  of  information, 
it  will  pick  out. . • 

Blberman:  You  will  then  calculate  a not  improbable  result,  that  you  can 

expect. 

Perms  The  difficulty  Is  that  the  parameters  that  seem  to  be  measured  era  not 
sufficient  to  describe  the  environment  sufficiently. 

Blberman i You'll  never  get  that  result,  not  except  on  a research  program. 
Then  you  can  only  expect  that  by  correlation.  They  are  not  even  sufficient 
to  base  a decision  on  thought  and  to  dateraina  what  the  general  condition  vau. 

Welsh:  Let  me  pt)Be  s specific  question  along  those  lines  just  to  follow 

up  on  that  point.  Do  you  think  that  it  would  be  worthwhile  to  try  to  decide 
this  question?  A model  should  not  be  a thing  that  gives  a monotonic  or  a 
specific  functional  relationship  bstwean  visibility  and  whatever  others  you 
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deem  appropriate,  and  an  absorption  at  10  micrometers.  Rather,  It  should 
just  let  you  arrive  at  a limited  set  of  statistical  parameters  which  would 
describe  a probability  distribution  of  a certain  attenuation  in  the  10 
micrometer  band.  Could  we  reach  some  conclusion  on  that  issue? 

Blbennan:  I'm  not  sure  that  I understand  your  question. 

Walsh:  Well,  the  question  is  the  following:  Suppose  that  I device  a 

p theme,  presumably  a functional  form,  that  said  when  the  relative  humidity 
and  visibility  and  temperature  and  whatever  else  that  seemed  to  be  Important 
was  in,  there  is  a calculation  that  leads  to  a single  specific  number  for  the 
performance  of  my  FLIR  system.  An  atmospheric  attenuation,  a single 
specific  number.  That's  one  approach.  The  oth-r  approach  would  be  one  that 
gave  me  a probability  of  performance,  a range  of  performance;  that  range 
would  be  narrower  than  the  total  range  of  performance  that  I could  get 
over  all  possible  parameters  of  the  atmosphere,  but  never-the-lesB  X would 
not  try  to  force  a single  number  out  of  it.  Is  it  worthwhile  to  try  to 
consider  that  kind  of  an  approach? 

Fenn:  I think  that  it  would  be  worthwhile  to  consider.  The  only  difficulty 

is  that  to  develop  that  probability  distribution  you  have  to  either  solve 
the  first  problem  that  you  described  or  in  other  words  you  would  have  to 
know  what  the  correlation  between  the  visible  visibility  and  the  IR  trans- 
mission would  be  (Walshs  You  think  it's  too  hard?)  to  develop  the  IR 
transmission  probability  distribution  from  the  visible  probability  distri- 
bution, or,  if  you  cannot  do  that,  to  go  out  and  measure  for  "um-teen"  yearc 
the  IR  transmission  to  develop  the  IR  distribution  from  that. 

Walsh:  You  say  you  prefer  the  single  functional  relationship? 

Fenn:  1 think  that  it  would  be  the  faster  apptoach. 

McClatchey:  I think  that  the  thing  is  that  we  have  to  commit  ourselves. 

I think  that  there  is  a reluctance  to  commit  oneself  to  the  scattering  pro- 
blem and  to  say  that,  based  on  the  best  information  that  we  have,  this  is 
the  answer.  How  certainly,  there  is  some  uncertainty  in  that. 
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Blbeman;  Well,  it  would  be  very  nice  If  the  meteorological  people  reported 
not  only  the  meteorological  visibility  which  could  be  very  short  but  also 
indicate  that  it.  was  snowing  and  when  they  do,  then  you  put  a step  function 
into  your  model  to  take  account  of  the  fact  that  it's  just  not  a very,  very 
heavy  haze,  but  it’s  actually  this  kind  of  a problem,  and  there's  one  other 
model  that  I've  heard  no  one  talk  about  in  thlB  connection.  We're  all 
talking  about  normal  peace  timo  (Bo-called)  conditions.  If  these  models 
are  going  to  be  of  value  and  we  take  a look  over  a hostile  area  like  a 
battlefield,  I think  that  your  conditions  are  like  your  volcanic  stuff  which 
occurs  immediately  after  the  first  3 or  4 shots  go  off  and  you're  going  to 
have  all  sorts  of  crud,  and  somebody  really  ought  to  establish  a project 
crud  and  look  through  It.  Mo,  seriously. 

Benedict;  This  is  exactly  the  aort  of  problems... 

Woodman:  I'd  like  to  make  a comment  about  the  models  that  we've  been 

discuaaing.  How  do  you  start?  What  kind  of  project  Bhould  we  have?  One 
of  the  things,  1 think,  that  we  should  address  Is  - do  we  collectively  feel 
that  there  is  credibility  in  the  Russian  approach  to  this  problem?  They've 
added  another  variable,  another  parameter  that  is  measurable  but  somewhat 
arbitrary.  But  1 think  that  we  should  decide  whether  we  should  invent  the 
wheel  or  perhaps  find  out  what  the  consensus  Is  about  the  Soviet  approach. 

They  present  data  and  claim  that  they  can  characterize  and  distinguish  a 
summer  haze  from  e stable  spring  or  summer  haze.  Mow,  1 have  difficulty 
reading  the  Russian  publications,  and  I'm  sure  all  of  you  do.  But  that 
seams  to  have  some  merit.  Mow,  are  we  going  to  work  around  that  or  purhapa 
consider  that,  aa  a possible  way  we  can  go  with  a new  aerosol  model.  I 
guess  that  I would  like  to  see  if  Dr.  Fenn  has  a comment  on  that  approach 
because  it  certainly  is  an  approach. 

Fenn:  Well,  I certainly  would  agree  that  it  is  possible  to  find  a correlation 

between  the  optical  properties  end  the  general  meteorological  conditions. 
Simply,  we  know  that  there  is  a physical  relationship.  We  do  not  know  too 
veil  what  the  physical  relations  are,  but  they  do  exist.  Consequently, 
there  has  to  be  e relation  between  the  optical  properties  and  th  general 


meteorological  conditions.  Whether  the  choice  of  tying  the  optical  pro- 
perties to  the  seasons  is  the  best  one  or  not,  I would  question  that.  Well, 
certainly  the  correlation  exists,  evidence  is  in  the  data,  and  it's  not 
surprising  that  one  finds  that  either.  We  do  know  on  a day  with  a clear 
arctic  polar  air  mass  moving  in,  the  visibility  is  going  to  be  much  better 
than  on  a hazy  summer  day. 

Bibermaa : That  depends  on  where  you  live.  If  you  live  in  South  Chicago, 

and  they  burn  soft  coal,  and  it's  clear,  and  you've  got  a nice  clear  arctic 
air  mass  moving  in.  It  might  be  clear  at  some  altitude  but  it  sure  isn't 
on  the  edge  of  Michigan  Avenue. 

Yes. . . 

You've  got  it  on  your  collar  and  on  your  shirt,  and  everything  else. 

Penn;  But  we  have  to  look  at  different  regimes, 

One  regime  is  the  typo  of  air  mass  - the  larger  scale  circulation.  Another 
regime,  hidden  from  that,  is  the  purely  local  effects  or  local  pollution 
from  a city;  if  we  are  east  of  the  city,  west  windB  may  deteriorata  the 
conditions.  If  we're  west  of  the  city,  an  east  wind  which  may  mean  clear 
arctic  air,  however,  in  that  particular  situation,  it  may  make  things  worse 
because  we  just  now  are  getting  into  internal  pollution  or  its  local  pollu- 
tion. So  this  is  one  reason  why  a simple  relationship  like  visibility 
conditions  for  summer  or  winter  may  be  meaningful  in  a certain  type  of 
environment,  but  may  not  be  applicable  at  all  to  another  type  of  environ- 
ment. Then,  one  has  to  look  at  that.  For  instance,  on  the  first  of  January, 
it  may  be  that  the  large  scale  circulation  regime  is  the  dominating  factor 
whereas,  on  the  third  of  January,  it  may  be  a local  phenomenon  that  determines 
the  local  optical  conditions,  and  these  are  things  that  one  would  have  to  con- 
sider. 1 do  not  think  that  we  know  enough  about  these  things  at  the  present 
time  to  describe  them  for  all  general  conditions.  We  know  that  these 
correlations  exist,  and  we  can  even  today  describe  them  fairly  reasonably 
for  certain  conditions  but  not  all  of  them. 

Woodman:  But  this  would  seem  to  be  a logical  starting  point  to  agree  that 

we  need  more  than  Just  the  visibility,  and  that  perhaps  some  type  of  classi- 
fication of  weather  type  and  air  mass  could  be  of  value  to  us.  And  I look 
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ac  what's  being  done  and  one  of  the  practices  in  modeling  is  to  incorporate 
a combination  of  a continental  and  maritime  aerosol,  and  if  they  model  the 
growth  of  the  particles  with  the  increase  of  the  relative  humidity.  That's 
fine,  it  helps  us  understand  the  phyrics,  but  the  person  who  is  using  the 
model  cannot  go  out  and  measure  the  ratio  of  continental/marltime  particles. 

Fenn:  No,  I think  that,  if  one  breaks  it  down  into  rural,  urban  and  maritime 
particles,  this  would  be  a description  that  would  describe  things  on  the  local 
regime  scale.  You  can  come  up  with  another  regime  in  which  you  relate  the 
optical  properties  to  the  seasons  or  air  mass.  And  certainly  this  tie  in 
to  the  air  masses  is  reflected  in  the  data  the  Russians  have  because  in 
summer,  in  winter  you  have  different  air  masses.  That's  why  they  get  these 
consistent  differences.  But  I don't  think  that  you  can  say  well  one  is  a 
substitute  for  the  other.  1 think  that  this  is  an  example  of  two  different 
regimes  that  control  these  optical  structures.  One  is  the  location  thing, 
the  urban  vs  rural;  or  continental  vs  maritime  but  these  are  local  things. 

The  other  reggae  that  may  dominate  things  is  the  circulation  of  air  masses. 
Those  have  to  be  looked  at  but  separately.  In  my  opinion  what  one  should  do 
is  on  a given  day  the  model  that  should  be  picked  should  be  the  summer  or  say 
the  subtropical  air  mass  type  model  combined  with  some  parameters  which  make 
it  fit  in  an  urban  environment. 

Biberman:  Well,  lets  ask  a question.  Who's  going  to  use  the  models  and  for 

what  purpose?  It  is  clear  that  nobody  is  going  to  plan  a mission  for  the  17th 
of  April,  1975  and  sit  down  and  calculate  like  crazy  what  the  scattering  is 
going  to  be  and  what  the  absorption  is  goind  to  be.  What  he  needs  to  do, 
however,  is  to  have  a background  and  to  know  that  over  some  period  of  time 
statistically  the  probability  to  do  such  and  such  arises  from  two  factors  - 
the  absorption  and  the  scatter  which  typically  in  April  is  going  to  give  him 
something  like  a 90%  probability  of  success  and  a 10%  probability  for  a 
failure.  And  that's  all  that  he  really  needs  for  his  planning  purposes. 

In  the  same  way  in  the  design  of  equipment,  there's  no  sense  in  designing 
a piece  of  equipment  or  a system  that's  going  to  see  beyond  the  curvature 
of  the  earth  if  you're  going  to  use  it  at  altitudes  below  50  feet.  And  in 
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the  sane  manner,  there’s  do  point  in  doing  range  capability  based  on  clear 
atmosphere  when  the  probability  is  less  than  so  much  that  you’re  ever  going 

to  achieve  them.  I think  that  this  saves  us  from  over-designing  and 

statistically  over-designing  equipment  and  it  prevents  us  from  spending 
fortunes  on  things  that  really  do  us  no  real  good.  That  does  not  mean  that 
we  do  not  build  some  special  purpose  dedicated  equipment  which  might  be 
used  when  the  opportunity  might  be  required  for  some  very,  very  high  priority 
mission,  but,  in  general,  when  we  build  a hundred  something  or  other  aircraft 
we  don't  put  100  or  so  of  these  aboard  each  aircraft  with  capability  of 
seeing  47,000  miles  because  the  probability  of  seeing  47,000  miles  In  any 
real  situation  is  going  to  be  terrible.  In  fact,  we'll  probably  see  four 
miles,  and  so  I think  what  we  want  to  do  is  use  our  models  to  make  that  kind 
of  executive  decisions.  Administrative  decisions  on  how  you’re  going  to 
distribute  your  money  and  your  talent  and  your  resources.  So  we  only  have 

to  collect  data  for  some  typical  year  and  the  reason  that  I used  1970  is 

because  that's  where  the  data  is. 

Fermi  I think  that  there  are  two  principal  types  of  problems  and  applications 
in  this  area.  One  is  what  you  described  as  the  statistical  data  - what  is  the 
probability  of  occurrence  of  a certain  type  of  data  at  certain  time  of  year 
and  under  certain  type  of  conditions  and  at  the  present  time  and  for  the 
foreseeable  future.  There  is  no  question  in  my  mind  at  all  that  the  only  way 
that  one  can  approach  that  problem  is  by  looking  at  the  statistical  data 
that  we  have  available.  And  those  are  visibility  from  the  regular  meteorological 
observations,  and  surface  visibility  and  eyeball  spectral  range.  Daytime  only, 
no  night  time  with  some  exceptions.  The  standard  meteorological  data,  the 
temperature  and  humidity  which  allows  ua  to  make  an  evaluation  of  the  1R 
transmission,  and  that’s  all  that  we  have  at  the  present  time. 

Blberman:  Well,  I have  records  which  give  me  meteorological  visibility  24 

hours  a day. 

Fenn;  Yes,  there  are  some  that  give  it  in  day  and  night  time. 

Selby:  What  you  don't  have  though  is  that  subject  to  visibility,  what  is 

the  attenuation  at  10  micrometers... 


Fean;  You’re  right;  there  are  a lot  of  things  that  you  don't  have.  The 
first  thing  you  don't  have  is  what  are  the  illumination  levels.  For  a lot  of 
systems,  specifically  night  time  systems,  what  is  the  distribution  of  the 
natural  illumination  level  especially  at  night  time? 

Biberman:  That's  fairly  easy  to  go  back  and  reference. 

Fenn:  Yes,  for  astronomical  data  but  you  have  no  information  on  what  the 

atmospheric  factors  do  to  them.  And  the  other  principal  night  area  of  no 
data  at  all  is  tha  IR  transmission.  The  extrapolation  from  the  vislblle  to 
IR  transmissions  Is  highly  questionable  and  might  be  invalid  all  together. 

Biberman:  But  as  far  as  the  distribution  of  light  and  its  level,  if  I know 

the  hour  of  the  day  and  the  day  of  the  month  and  the  year  so  that  I can  get 

the  ephemeris  and  I know  the  cloud  cover  and  a few  things  like  that,  I think 

that  1 can  calculate  you  the  incident  illumination  and  the  spectral  distri- 
bution within  quite  a reasonable  factor. 

Fenn:  Yes,  but  another  point  of  view,  another  .’  tea  in  which  one  has  no  data 

at  all  is  the  contrast-  reduction  in  the  visible.  For  any  kind  of  system  that 
works  on  the  contrast,  it's  absolutely  inadequate  to  know  what  the  trans- 
mission is.  You  have  to  know  what  the  contrast  reduction  is  end  there's 
no  measurement  of  that  available.  But  for  1975  and  1976,  and  several  years 
after  that  we're  not  going  to  have  any  additional  data  of  that  nature.  All 

we're  going  to  have  is  the  surface  visibility,  so  the  only... 

Biberman:  What  about  the  Dutch  effort? 

Fenn:  Well,  that  was  a one  year  measurement  program. 

Biberman:  Well,  itB  pretty  good,  isn't  it? 

Fenn:  Yes. 

Biberman:  And  one  year  gives  you  a fair  amount  of  statistics. 

Fenn:  I wish  we  had  maybe  a dozen  of  these  measurements. 

Gallagher : Let's  go  on  to  the  next  paper  that  we  had  this  morning.  People 

keep  coming  back  to  this  problem,  and  it  looks  like  one  that  we're  not  going 


200 


to  get  rid  of  the  rest  cif  the  afternoon.  Does  anyone  have  son®  coaMatu 
Wolfhardt's  paper  on  Rocket? 

Benedict;  I'd  like  to  comment.  He  certainly  raised  a very  Important  point 


that  the  emission  from  plumes  is  different  from  the  continuum  and  he  aaid 
what  are  needed  are  better  data.  This  also  ties  In,  of  course,  to  the  various 
modeling  of  things  that  was  made  when  the  emission  and  absorption  were  baaed 
on  calculations  using  lines  from  the  AFCRL  line  list.  The  point  that  l want 
to  make  is  that  the  AFCRL  line  list  was  designed  solely  for  use  at  roam 
temperature. 

The  actual  knowledge  of  the  data  that  could  be  made  to  go  Into  a line 
list  that  would  give  both  the  emission  and  the  absorption  ot  water  vapor,  C02 
and  CO  at  much  higher  temperature  does  exist.  With  regard  to  CO,  ehere  already 
is  such  an  excellent  line  list  although  the  frequencies  are  not  the  meet 
up  to  date  but  those  could  be  easily  corrected.  That's  the  one  that's  due 
to  Birch.  As  far  as  CO^  goes,  here  the  data  do  not  exist  to  take  cere  of 
everything  up  to  quite  high  temperature  but  I would  like  to  make  the  point 
that  there  are  so  many  vibrational  states  that  are  excited  with  the  eery  low 
bending  frequency  of  C02  that  if  you're  talking  about  temperature  above  1000* 
the  density  of  lines  is  almost  equal  to  the  Doppler  width  so  that  ae  far  as 
C02  goes,  you  have  essentially  a continuum  source  with  a certain  distribution 

of  course  rather  than  a many  line  source. 

The  remaining  question  of  the  H20  is  one  that  the  data  ate  noc  ae 
complete  as  I would  like  them  to  be  but  they  certainly  do  p.ralc  a very 
considerable  extention  of  the  list  of  lines  rather  of  some  of  the  line* 
that  are  in  the  AFCRL  list.  As  far  as  all  the  levels,  up  to  6000  cm  of 
excitation  are  probably  known  to  an  accuracy  that  will  give  you  the  frequency 
of  Che  line  to  within  0.1  cm"1,  and  this  should  be  adequate  for  taking  care 
of  the  largest  amount  of  the  emission  from  plumes,  in  the  ..7  micrometer,  the 
1.9  micrometer,  6 micrometer  and  pure  rotation  regions  so  its  just  ■ witter 
of  getting  the  data  together  and  making  a very  much  expanded 
basic  data  that  can  be  used  either  in  actual  line-by-iine  calculatl 
U-  iei  calculations  will  be  there.  The  quantitative  accuracy  *a  far  a.»  line 
strength  of  H20  goes  may  not  be  terribly  good,  but,  from  all  the  examples 
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that  we  have  seen  this  morning,  the  data  in  the  AFRCL  line  list  are  not  • 
terribly  good  for  individual  lines  but  when  you  get  to  the  broadband  compari- 
sons of  models  and  data  the  overall  thing  is  quite  good  and  I believe  the  same 
thing  can  be  done  for  water  vapor  up  to  at  least  a 1000  degrees  without  too 
much  trouble.  To  go  beyond  a 1000  degrees,  one  needs  further  effort  on 
high  resolution  studies  on  flame  sources  in  the  longer  wavelength  region. 

Such  studies  do  exist  in  the  region  down  to  short  wavelength  and  up  to  4 
micrometers.  For  the  record,  my  name’s  Benedict  from  the  University  of 
Maryland. 

3iberman:  Don't  we  have  to  take  a look  at  the  criticality  of  the  problem 
from  two  different  points  of  view.  We  say  that  there  are  plumes  that  are 
very  hot.  And  if  we  take  a look  at  large  by-pass  engines  for  instance,  the 
effluence  is  not  extremely  hot  but  if  we  take  a look  at  LOX-hydrazine  or  some 
of  the  other  fuel  oxidizer  combinations,  that  can  be  used  in  some  of  our 
high  performance  things  like  some  of  our  larger  rockets,  then  the  exhaust 
temperatures  are  very  high  and  the  problems  that  we  raised  this  morning 
are  really  raised.  But  are  the  problems  really  so  bad  when  we  take  a look 
at  fairly  high  by-pass  engines? 

Long:  What's  the  temperature? 

Biberman:  A few  hundred  degrees  C. 

Long:  I don't  think  that  it's  really  much  of  a problem  in  handling  that 

problem. 

Benedict:  No.  There  are  two  types  of  input  data  that  are  needed  In  the 

analysis  of  hot  water  vapor  spectra.  What  I have  been  working  on  have  been 
the  hot  thin  flames  and  the  spectra  of  the  3un  spots  and  these  are  of  course 
sources  of  3000-3800  degrees  and  once  that  you  know  that  then  going  down  to 
low  temperatures  presents  no  problem  at  all. 

McClatchey : The  other  kind  of  traffic  problem  that  that  presents  is  why  you 

don't  want  to  add  stuff  unless  you  really  have  to  because  I believe  from  a 
computational  point  of  view  it  could  become  one  hell  of  a nightmare. 

Benedict:  Exactly. 


McClatchey : Just  trying  to  deal  with  an  order  of  magnitude  more  lines  or 

some  crazy  thing.  That  could  be  a real  problem  so  I think  that  we  have  to  be 
very  careful  that  we  establish  that  something  is  really  net  issary  or  we  have... 

Biberman:  The  adequacy  of  the  model  and  the  adequacy  of  the  data  bank 

concerning  the  atlas  of  lines  and  line  strengths.  I think  that  you  could 
use  a fairly  narrow  set  of  data  for  many  of  the  studies  that  we  want  to  do, 
on  aircraft  plumes,  on  soma  of  the  modern  engines  and  then  move  into  e 
completely  different  and  more  tough  problem  when  we  start  moving  up  drastically. 

Burch;  But  you  know  In  some  of  the  cases  I agree  with  McClatchey  in  that  the 
computation  gets  very  difficult,  but  can't  you  build  into  your  program  some- 
thing that  just  ignores  the  lower  strength  lines  for  applications  where  you 
don*  t need  them. 

Long:  Oh,  yeah,  you  could  put  them  in  and  not  worry  about  a horrendous 

computational  problem  in  those  cases  where  you  don't  need  them. 

Biberman:  Wait  a minute,  I don't  think  that  that's  really  right  because  the 

lover  strength  lines  when  they  predominate  in  number  may  be  the  most  important 
thing  that  you  have  to  consider. 

Benedict:  What  you  have  to  know  is  the  number  of  lines  in  a given  intensity 

range  and  a given  frequency  range. 

Biberman:  You  Just  can't  make  an  arbitrary  decision  by  any  means. 

Benedict:  There  are  two  different  types  of  pieces  of  information  that  should 

be  made  available;  one  is  the  overall  one  in  which  you  know  just  how  many 
lines  there  are  and  what  their  range  of  strength  is  and  approximately  where  they 
are.  If  you're  really  trying  to  go  into  a detailed  comparison  of  a specific 
region  then  you  ought  to  know  where  precisely  each  line  is,  each  strong  line 
is. 

Long:  The  point  that  I was  trying  to  make  was  that  it  may  be  a horrendous 

problem  in  assemblying  this  data  and  on  that  basis  you  take  care  in  how 
many  you  assemble.  Just  on  the  basis  of  the  computation,  you  can  always  get 
the  program  that  allows  for  that.  But  if  you  didn't  over  compute... 


Benedict : As  an  example  of  this,  there  is  a paper  in  the  Astrophysical 
literature  by  Amat  who  does  essentially  this  on  the  question  of  hot  water 
vapor  in  stars.  It's  seven  or  eight  years  old,  its  fairly  solid. 

Gallagher:  How  about  molecules  other  than  those  that  are  on  the  list?  Can 

they  contribute  in  this  area? 

Benedict:  In  the  hot  flame  region,  for  instance,  there  are  a number  of  strong 

lines  that  are  due  to  OH  for  example  and  thuBe  arv  very  valuable  in  that  they 
can  be  detected  at  these  temperatures.  But  ordinary  type  hydrocarbon 
flames  and  so  on,  we  are  left  with  the  Bimple  type  of  molecules  that  we  know 
pretty  well. 

Long:  There  are  some  other  things  that  could  be  added  to  this  particular 

data  set  such  as  UC1  which  has  shown  up  in  some  types  of  exhaust  products 
and  I am  sure  that  there  are  probably  some  otherB  too  that  are  still  in  this 
category  of  simple  molecules,  that  I am  sure  are  easy  to  put  together  and 
have  the  capability... 

Benedict:  They  don’t  clutter  up  the  list. 

Long:  Yes,  they  have  a small  number  of  lines. 

Benedict:  I chink  for  instance  adding  OH  if  there  is  a great  deal  of  importance 

and  interest  in  the  emission  problem  should  be  one  of  our  first  priorities. 

Long:  We  are  trying  to  find  out  with  respect  to  the  DF  laser  there  are  a 
couple  of  lines  that  we  see  absorption  in  the  most  pure  nitrogen  gas  that 
we  can  purchase  and  we're  trying  to  figure  out  what  it  is. 

Burch:  What  wavelength  is  it? 

Long:  You  mean  exactly  which  laser  line  is  it? 

• * a 

Long:  Yeah,  in  that  general  area.  Do  you  want  to  make  a guess? 

Burch:  Wall,  1 don’t  have  an  answer  but  I would  like  to  look  for  it  when  I 

get  home. 

Long:  I'll  call  you.  1 forgot  what  line  it  is. 
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Selby;  In  addition  to  the  individual  lines,  you've  got  a problem  with  the 
line  shapes  still  at  ordinary  temperatures  and  further  complications  at 
higher  temperatures.  You  may  need  to  put  into  our  program...  Well,  the 
higher  the  temperature,  the  higher  the  pressure  at  which  the  Doppler  shape 
where  you  have  to  go  to  the  right  regime  for  the  Lo**entz  width.  Also  there 
are,  of  course,  no  very  accurate  measurements  of  the  ...dependence  of  Lorentz 
width  on  temperature , but  I think  that  the  few  measurements  that  there  are 
Indicate  that  the  higher  the  temperature  the  less  the  effect  of  anything 
other  than  the  billiard  balls  where  the  collision  diameters  are  of  importance. 
So  it's  pretty  safe  to  take  a constant  collision  diameter  at  any  temperature 
above  a thousand  degrees. 

Gallagher;  How  about  the  paper  that  Corcoran  gave.  Have  you  comments  on 
that  as  far  as  the  laser? 

Rohde;  I have  a question  to  address  that  Dr.  Biberman  has  already  hinted  at, 
and  that  is  when  Rudy  Buser  comes  to  me  and  asks  me  whether  or  not  a laser 
system  is  going  to  work  under  certain  conditions,  we  also  have  to  consider  the 
battlefield  situations  and  even  if  I got  all  the  codeB  to  work  under  the 
statistics  and  the  weather  conditions  for  all  over  the  world.  In  a different 
battlefield  situation,  do  I know  if  the  laser  if  going  to  function  properly 
or  am  I going  to  wind  up  with  sc-ie  surprises  that  I didn't  expect?  So  before 
I have  to  go  back  to  him  with  x answer  I would  like  to  get  some  information 
as  to  whether  it  is  possible  to  make  some  transmission  measurements  in  a work- 
ing situation  out  at  Ft.  Carson  or  some  place  and  to  really  see  what  the 
conditions  are  before  they  start  this  game  what  the  visibilities  are,  what  the 
transmissions  are  and  now  make  the  same  measurements  while  the  game  Is  going 
on  so  that  1 could  have  a feeling  for  how  much  these  conditions  could  possibly 
change  and  whether  because  of  the  new  molecules  which  are  coming  up  or  be- 
cause of  the  shells  exploding  or  because  of  the  vehicles  that  are  in  the  area 
that  all  the  predictions  that  people  are  giving  me  are  going  right  out  the 
window.  That's  the  final  number.  You  know,  when  I go  back  to  my  boss  and 
tell  him  yes  this  system  will  work  in  a battlefield  situation  and  he  can 
believe  it.  And  that's  a problem  which  I cannot  answer. 
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?.-!  Berman:  Well,  "lve  seen  a number  of  Army  films  which  are  trying  to  demon- 

strate how  gocdsomething  is  andone  of  the  typical  things  that  they  usually  do 
is  show  an  infrared  set  looking  for  something  or  other,  and  it's  behind  a 
howitzer  or  some  other  kind  of  gun  and  some  one  pulls  a lanyard  and  the 
thing  goes  boom.  And  then  there's  a television  camera,  an  Infrared  camera  of 
some  sort  and  you  sit  and  watch  as  a function  of  time  how  Long  it  takes  for 
the  horizon  to  come  into  view,  or  the  thing  down  the  road  or  whatever  it 
happens  to  be.  Now  there's  a fair  amount  of  that  stuff  that  has  demonstrated 
the  problem  about  which  you're  talking,  but  which  has  not  been  quantified 
and  it  seems  to  me  that  it  is  not  a terribly  difficult  job  to  work  in  as  a 
companion  to  some  other  investigation  at  the  appropriate  time  and  I think 
that  you  certainly  need  to.  But  there  are  many  sources.  One  is  just  the 
muzzle  blast  and  all  that  kind  of  stuff  and  what  it  does,  it  blows  the  top  of 
the  earth  off,  hurls  it  into  the  air  and  settles  down  slowly.  And  then  the 
other  thing  is  that  when  it  gets  to  the  far  end,  it  digs  a hole  and  throws  a 
mine  up  into  the  air,  and  that  settles  down  slowly.  And  I think  if  there  is 
a serious  engagement  you  have  to  make  some  few  measurements  and  some  pre- 
dictions about  the  range  of  difficulties  that  you  get  into. 

Rohde;  Precisely. 

Gallagher;  How  about  the  questions  brought  up  about  linewdith  parameters? 

When  one  is  looking  at  narrow  laser  lines,  are  the  linowidth  parameters 
that  well  known  lu  say  that  vc  have  a particular  transmission  at  these 
wavelengths.  I think  the  people  up  a Cambridge  have  run  into  the  problems 
that  John  has  mentioned  just  looking  at  SOg  for  instance,  for  which  Clough  had 
trouble  with  the  Lincoln  Lab  people.  On  the  baBis  of  predicted  linewidths, 
the  question  c.ame  up  whether  they  were  getting  different  linewidths  in  their 
experiments . 

Benedict;  I don't  know  the  situation  on  S02  but  it  certainly  is  true  in  the 
experiences  that  1 have  had  that  it  is  never  possible  to  calculate  anything 
as  well  as  it  is  going  to  be  measured  when  it  can  be  measured  well.  When 
you're  doing  the  kind  of  measurements  like  Ron  Long  is  doing,  for  example, 
where  you  have  laser  lines  of  known  widths  and  known  frequencies  then  you 
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see  whether  the  position  as  calculated  for  a simple  molecule  like  N20,  the 
absorption  obeys,  he  tells  me  that  it  obeys  very  nicely.  When  it  comes  to  a 
very  difficult  molecule  like  1^0,  the  calculations  don't  agree  that  well. 

There  are  problems  remaining  in  linewidths,  particularly  in  the  line  '.apes 
on  the  wings  of  the  HgO  lines  which  of  course  leads  us  into  this  other  mess 
that  we  have  heard  so  much  about,  of  the  continuum  that  we  heard  so  much  about 
this  morning.  But  as  far  as  the  general  order  of  magnitude  thing,  20-3(Kr)0% 
or  so  of  the  present  calculations  are  gocd.  For  new  molecules  such  an  S0? 
where  there  have  not  been  as  far  as  1 know  any  accurate  linewidth  measure- 
ment, we  shouldn't  be  too  surprised  that  the  original  oneB  come  nut  wrong. 

McClatchey:  Well,  it  just  kind  of  occurred  regarding  the  data,  the  measure- 

ments of  the  Lincoln  Lab  people  on  water,  for  example,  have  not.  really  been 
introduced  empirically  into  the  list.  There  are,  of  course,  only  a limited 
nurdjer  of  lines  that  have  been  measured.  But  the  philosophy  of  having  these 
narrover  lines  in  high  J have  been  to  some  degree  taken  into  the  calculations. 

Benedict : Well,  this  is  a second  problem  of  course.  The  measurements  of  the. 

Lincoln  Lab  people  indicate  that  a limited  number  of  H^O  lines  are  very  much 
narrower  than  the  ones  in  our  original  list.  But  as  Bob  has  just  said,  It 
is  not  of  major  importance  except  to  those  lasers  which  happen  to  operate 
near  these  lines. 

Long:  How  do  you  take  that  into  account  now?  I»  it  in  the  data  tabulated  or 
In  the  program? 

McClatchey:  What  I'm  saying  is  that  we  hnven't  token  it  into  account  in 

detail. 

Long:  Oh,  you  said  to  some  limited  extent', 

McClatchey:  Well,  in  the  sense  of  the  calculated  linew.ldf.hn . 

Bcneiict : The  current  listing,  I don’t  know  if  its  on  your  lists  for  general 

distributions,  but  Clough  put  in  narrower  linen  for  the  pure  rotation  region, 
and  I heard  from  Doug  Woods  who  works  at  MAI  up  in  Ann  Arbor  thAt  they  had 
checked  some  of  those  measurements.  Their  measurements  were  indeed  narrower 
than  our  old  tape  but  they  were  wider  than  the  guesses  that  we  made  on  the 


basis  of  the  Lincoln  Lab  extrapolations  on  the  new  tape.s.  So  again  its  got 
to  be  aerJ ured  before  you  can  say  anything.  But  its  an  awful  lot  of  measur- 
ing if  you  want  to  know  everything.  My  point  is  as  far  as  specified 
transmission  for  specific  laser  lines,  that's  not  hard  to  measure.  For 
overall  statistical  transmission,  current  data  seem  to  do  a pretty  good  job. 

Gallagher:  I believe  I.  just  saw  Dr.  Long  flinch. 

Long:  Well,  you  said  it  wasn't  hard;  I think  that  it's  awful  hard, 

Benedict:  I didn’t  mean  to  detract  from  your  actual  work.  1 thought  that 
I was  giving  you  a compliment. 

Gallagher t Let's  move  on  to  the  McClatchey-Selby  paper  as  I think  that  there 
will  probably  be  several  comments  there. 

McClatchey.  Actually  there  are  two  comments  that  I would  like  to  make  that 
clear  up  the  situation.  (1)  Regarding  the  2.7  micron  region,  and  especially 
this  quest ionof  the  modification  of  the  water  vapor  data,  in  LOTRAN  and 
its  connection  with  line-by-line  calculation,  what  happens  there  is  that  in 
the  first  case,  the  LOTRAN  coefficients  depend  on  line-by-line  calculations 
degraded  to  the  appropriate  spectral  resolution.  And  what  we  found  was  that 
the  result  of  that  calculation  did  not  agree  as  well  as  we  would  like  with 
measurements.  So  in  the  iteration  which  amounts  to  the  input  to  LOTRAN  111, 
at  this  point,  it  was  decided  to  base  the  coefficients  on  measurements  rather 
than  on  the  line-by-line  calculations.  I say  this  because  Charles  Randall 
made  some  reference  to  discrepancies  in  the  line-by-line  calculations,  and 
some  of  your  other  curves.  A number  of  them  which...  1 guess  what  I'm  saying 
is  that  when  you  taka  all  of  the  lines  into  account,  as  they  are  down  on  the 
listing,  ns  individually  determined  by  various  measurements,  and  you  then  do 
the  calculations  for  low  resolution,  that  you  run  into  some  problems  whether 
these  problems  have  to  do  with  lineshapes  or  some  other  things.  I don't 
really  know  what  is  all  involved.  The  line-by-line  calculations  themselves 
do  not  seen  to  work  out  as  well  as  one  would  hope. 

Long:  Which  specific  region  are  you  talking  about? 

McClatchey : 2.7  microns. 
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Benedict : In  specific  narrow  regions,  the  overall  picture  seems  all  right. 


Long:  The  thrust  of  what  I was  showing  was  not  to  demonstrate  that  it  was 

bad,  but  that  it  was  good. 

McClutchey : But  I'm  trying  to  get  at  the  reason  for  this  modification  and 

what  we  did  to  modify  it.  At  least  in  the  context  of  the  LOTRAN  code,  it 
appears  that  we  needed  to  adjust  the  parameters,  the  coefficients,  more 
nearly  to  match  the  experimental  data  rather  than  to  match  coefficients 
determined  by  the  line-by-line  calculation.  There  is  a lot  of  interest  in 
this  region  that's  why  we  might  not  ordinarily  have  worried  about  it,  and 
of  course,  the  discrepancies  are  not  that  big  you  might  say,  but  from  the 
point  of  view  of  addressing  a specific  problem  they  asked  of  us  regarding  the 
2.7  micron  region,  we  made  these  adjustments  and  we  think  it  gives  a 
better  fit  for  the  range  of  atmospheric  conditions  that  we  used. 

Long:  On  a different  topic,  I would  comment  on  the  computation  of  laser 

lines  transmission,  that  one  of  the  problems  that  we  have  had  is  that  we  do 
not  know  where  the  laser  lines  are  well  enough  and  in  that  regard  we  are  now 
measuring  some  things.  Rao  of  Ohio  State  is  measuring  the  DF  lines,  he  has 
completed  the  measurement  of  the  DF  line  positions,  he  has  completed  measur- 
ing the  HF  lines  positions.  Of  course,  he  did  CO  sometime  ago  and  we  hope  to 
do  HCl  if  the  laser  works  so  we'll  try  that  in  a little  while.  And  his 
accuracy  is  now  about  0.002  cm  *.  Now  there  is  some  other  work  around,  I 
think  at  Lincoln  Laboratory,  frequency  heterodyning  or  what  have  you,  they 
talk  about  an  order  of  magnitude  better  accuracy  than  this.  I think  for 
most  of  our  purposes  that  this  will  be  good  enough.  And  in  some  cases  we 
found  that  this  was  quite  important  In  getting  a proper  comparison  with  the 
line-by-line  tape  deck.  When  we  had  the  right  lsser  frequency,  we  got  a 
much  better  comparison  than  we  did  before, 

McClatchey : OK,  I just  want  to  pursue  this  since  I heard  someone  mention 

this  to  me  last  week,  that  in  your  report  that  you  wrote,  you  showed  some 
substantial  discrepancy  between  calculations  based  on  the  tape,  especially 
in  the  case  of  HDO  lines  and  the  laser  measurements.  And  then  I was  told 
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that  because  of  this  some  shift  in  frequency  that  may  be  this  isn't  so  bad. 

Is  that  a correct  statement? 

Long:  I gave  the  information  to  Bill  here.  It  improved  it  slightly  with 

respect  to  EDO.  With  NjO,  it  brought  it  into  absolutely  remarkable  agree- 
ment between  the  measured  and  calculations.  But  with  HDO,  there  still 
seems  to  be  a problem  there. 

Benedict:  I suspect  that  the  HDO  intensities  that  we  put  in  are  too  low. 

There  are  measurements  going  on  at  SAI  that  Meredith  and  Woods  are  doing, 
which  I hope  will  clarify  this. 

Gallagher?  Do  the  results  of  Rao  on  CO  compare  with  the  heterodyne  mea- 
surements that  are  being  done  at  MIT? 

Long:  In  the  one  case  that  I know  of  in  which  a comparison  was  made, 
between  Javan's  measurement  and  his,  it  was  really  remarkable.  It  was  within 
the  Doppler  width  of  the  laser  line. 

Benedict:  As  far  as  CO  goes,  everything  I think  frequency  wise  is  in 

beautiful  shape  because  there  are  a number  of  laser  frequencies  that  have 

been  fitted.  There  are  Rao's  measurements,  there  are  the  Golash  and  Lee 

measurements.  They  used  the  interferometer  and  there  are  the  solar  spectrum 

measurements  up  to  very  high  J.  As  far  as  the  frequencies  of  the  laser  lines, 

-4 

they  are  certainly  known  to  10  wave  numbers.  I don't  think  that  we  can 
quibble  about  that. 

Selby:  Is  there  a tabulation  of  the  CO  laser,  Bill?  The  line  positions  and 
intensities  that  enter  this? 

Benedict:  The  new  constants  for  calculating  the  CO  line  positions  that  have 

been  published  are  the  ones  of  Kildalling  and  Ross.  These  aren't  quite  as 
good  as  they  should  be  if  you  go  up  to  very  high  J's.  If  you're  only  inter- 
ested in  the  absorption  problem  and  not  going  to  the  emission  problem  in 
vary  hot  flames,  then  this  is  good  enough.  The  frequencies  can  then  be 
improved  by  using  the  constants  of  Kildalling  and  Ross.  The  intensities  that 
are  in  Conte's  6 or  7 year  old  paper  1 think  are  just  about  as  good  as  they 
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can  be.  Nothing  has  come  up  since  then  to  cause  a particular  improvement, 
and  this  la  a very  nice  table  because  it  has  the  values  at  a number  of 
temperatures. 

Gallagher:  Are  there  any  other  comments  or  questions  on  the  McClatchey- 

Selby  paper?  Do  they  have  anything  to  say  about  their  own  work  that  might 
concern  the  problems  that  bear  more  looking  at? 

Harris : I have  one  question  on  this  thing  here...  John  has  mentioned  that 

he  Inserted  spectral  aerosol  attenuation  data.  How  much  detail  are  you 
going  to  put  in  on  that? 

Selby:  I showed  one  slide  that  showed  the  refractive  index  of  aerosol  as  a 

function  of  wavelength. 

Harris:  The  attenuation  coefficient  as  a function  of  wavelength? 

Selby:  Yes,  right.  So  that  was  digitized  and  put  into  LOTRAN  111. 

Harris:  That  was  for  HjO? 

Se lby : That  was  for  a composite  of  dust  and  water  soluable  materials. 

Harris:  Oh,  now  that  was  the  crux  of  the  whole  matter. 

Would  you  like  to  say  what  you  use  in  your  program.  We've  had  a dialogue  on 
this  thing  with  a couple  of  the  people  who  are  here  now.  That  the  atmosphere 
is  so  enormously  complex  that  we're  hard  pressed  to  specify  what  the  relative 
parameters  are  and  we  need  to  know  and  we  do  not  know  this  yet  to  be  able 
to  predict  the  things  that  We  want  to  know  for  calculations.  Now  we  can 
measure  the  thing  in  a given  case,  measuring  all  the  things  we  can  probably  do 
it;  measure  the  Mle  scattering,  measure  the  things  as  a function  of  wave- 
length, measuring  all  sorts  of  things,  we  can  probably  do  it.  But  we  don't 
know  because  the  atmosphere  is  so  complex  and  it  is  so  bad  that  there  is 
hardly  anything  left  on  which,  by  means  of  a few  things  we  can  predict  almost 
anything,  because  of  the  inherent  variability.  You  mentioned  two  or  three 
times,  and  about  the  maritime  aerosol  and  this  becomes  a variable,  the 
continental  is  a variable  and  the  earth  is  a variable  and  then  you  get 
mixtures.  Not  only  the  seasonal  types  of  things  and  the  air  mass  but  a lot 
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of  other  things  that  we  haven't  even  considered.  Now  the  simplest  thing  to 
do  would  be  to  get  the  information  from  using  visual  range  or  from  using  an 
Integrating  nephelometer  which  gives  you  pretty  good  visual  range  at  a point. 
But  that  doesn't  do  it  all  because  you  get  various  kinds  of  Mie  scattering 
from  all  sorts  of  combinations  of  things.  So  the  important  thing  for  your 
model  is  to  take  a reasonably  good  combination  of  these  things  which  would 
work  in  a rather  wide  variety  of  circumstances  and  some  sort  of  mixture 
of  things.  But  it  makes  a great  deal  of  difference  whether  you  are  going  to 
take  horizontal  measurements  and  whether  you're  over  ocean  or  you're  over 
land.  Or  If  Its  going  to  be  over  a vertical  path  and  “.hen  when  you  get 
above  a certain  altitude  It  is  the  aame  as  a continental  atmosphere  anyway. 
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The  computing  session  was  attended  by  only  two  members  of 
of  the  committee,  thus  the  exchange  of  technical  information 
was  limited  and  the  session  was  ended  early  to  adjourn  and  go 
into  the  experimental  session. 

Two  points  were  covered,  however;  one  was  the  experimental 
work  on  infra-red  laser  fusing  and  the  other  was  the  theoretical 
work  on  IR  transmission  calculations  being  carried  out  at  Avco. 
The  fusing  work  is  being  done  at  NWC,  and  the  program  calculates 
the  return  of  laser  energy  from  clouds  and  target  (which  in  this 
case  is  an  airplane).  This  return  energy  is  used  as  a designator 
for  a C02  homing  missile.  The  program  1b  up  and  running  and 
some  of  the  results  were  pointed  out.  One  of  the  most  signif- 
icant is  the  problem  of  backscattered  energy  off  the  clouds, 
the  missile  at  times  sees  more  energy  from  this  source  than  from 
the  target. 

Work  on  the  code  was  funded  by  the  NWC  and  the  finished 
program  does  nut  make  use  of  LOWTRAN  or  the  Aerospace  code. 

The  second  subject  discussed  was  a program  written  at  Avco 
on  atmospheric  IR  transmission.  This  program,  uses  a Goody 
model  for  the  line  strengths  and  half-widths  for  a homogeneous 
path.  A procedure  to  obtain  an  equivalent  homogeneous  path 
in  an  Inhomogeneous  atmosphere  uses  the  Curtis-Godson  approxi- 
mation. To  use  this,  tables  of  line  strengths  over  lines  as 


30  Km.  to  57  Km.  doppler  effects  were  taken  into  account.  Re- 
fractive effects  were  not  accounted  for,  but  earth  curvature 
effects  and  atmospheric  emission  were.  Wavelength  resolution 
cells  are  given  by  X/AX  ■ 100  (0.1  y at  10  y)  from  8 to  25 
microns  from  H2,  C02,  0^,  N20,  CH^,  CO  and  HNO^.  The  program 
was  made  to  run  on  an  IBM  360-44  of  20K  word  capacity.  Test 
runs  were  made  and  compared  favorably  with  experimental  measure- 
ments . 

The  remainder  of  the  discussion  session  was  concerned 
heavily  v/ith  the  effects  of  aerosol.  It  was  pointed  out  that 
it  is  very  difficult  to  make  aerosol  measurements  in  which  the 
experimental  conditions  are  well  controlled.  Combination  of 
aerosols  and  variations  along  the  path  of  aerosol  concentrations 
present  problems.  There  is  a great  need  for  high  precision 
transmission  data  to  make  range  predictions.  Calculations  are 
needed  on  different  types  of  aerosols.  The  models  are  useful 
for  averages,  but  a need  has  been  expressed  for  experimental 
measurements.  The  sensitivity  of  the  data  is  extremely  impor- 
tant. Thus,  it  is  Indicated  that  a 2-1  pread  of  data  presents 
a factor  of  10-1  in  utility  or  design  of  device. 

From  the  discussion  of  the  meeting,  it  is  evident  that 
there  is  an  important  need  for  experimental  data  in  all  phases 
of  propagation  as  inputs  to  the  calculations.  Simultaneous 
precise  meterological  data  is  needed  for  all  propagation  meas- 
urements. 

The  last  part  of  the  recordings  of  the  open  discussions 
was  not  clear  enough  to  present  detailed  information. 
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